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АВ5ТКАСТ 


A transient three-dimensional heat conduction code was devel- 
oped using finite differences. A stability restriction on the time step 
was avoided using a technique proposed by Brian. Computations from 
the code were validated using both the explicit technique and an 
available closed-form solution for small times. The maximum error was 
found to be within 0.019 percent for an 11 x 11 x 11 grid and time 
step of 17.117 seconds. The total CPU time to carry out the computa- 
tions up to 3,600 seconds using Brian’s technique was six times that 
required for the explicit technique with the same time step of 17.117 
seconds. However, as the time step was increased without altering the 
geometry, the CPU time using Brian’s technique decreased and was 
less than that used in the explicit technique for time steps larger than 
110 seconds. The validated code was also used in the analysis of the 


transient thermal response of a component on an orbiting spacecraft. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in 
this research may not have been exercised for all cases of interest. 
While every effort has been made, within the time available, to ensure 
that the programs are free of computational and logic errors, they 
cannot be considered validated. Any application of these programs 


without additional verification is at the risk of the user. 


1V 


П. 


ПІ. 


TABLE OF CONTENTS 


cesse clo DOS sem e TEES 1 
С ЕИ ENITSOROEPRODBIEENCT S eere ниво Рени 0). 1 
ЕЕ ЕЕЕ ЕЕ НОО ннен 2 
Е 32 £107 5250 93 1 00010) QD (0) D) 2 Ps Cc насоките 3 
СТО ИЕ ПУТЕР о иеа ненае ене ие З 
И ЕЕ Е НОО ОЕ аан Ынды н дыы ы... ынна 4 
1. Derivation of Nodal Temperature Equations.......................... 4 
2. Restrictions on Time Step in the Explicit Technique...... 7 
ЕС] ао ее scia eire eo enne 9 
DE И У ЕЕ C SENI CIE u s ea sam ARR 9 
l. Derivation of Nodal Temperature Equations....................... 11 
2. Solution of the Nodal Temperature Equations................... 17 
EB DEC REC CGC3lCnlStionsS нне 18 
ОТОО а ЗОО анала Ыбы азаннан... 20 
И СЕБЕШЕЮ СЕОМЕТКА АМО CONDITIONS................................ 20 
БОЛЕЕ АЕСОМРОТАТОМ на на 22 
Gu Е ION  COMPUTATIONS................ uu uuu cunei eu opone sue tnn 24 
ое УУТТАН 35 
K Up ненада наннан 35 


В ORBITING SATELLITE COMPONENT m eee 35 


1. Geometry and шина Соте ои co TM 35 

2. Numerical Caleülations mna I ass s ss 41 
У. ВЕЗОГТЬ... ао nn o EE 45 
VI. СООМСЕОЗОМ®.....:...... 5 --се-с-ест ООШ 46 
VII. КЕСОММЕМПАТІГОМӘ............ ИИИ 47 
APPENDIX A EXPLICIT TECHNIQUE NODE EQUATONS a ma 48 
APPENDIX B BRIAN'S TECHNIQUE NODE EQUATIONS 28 57 
APPENDIX C РЕОСЕАМ ЕХРИС nn 0-0 85 
APPENDIX D PROGRAM BRIAN. И 107 
APPENDIX E PROGRAM УАБПЫ 222 оН 187 
LIST OF REFERENCES ............... R P а UT НИЯ 191 
INITIAL DISTRIBUTION LIST ..... ev n PP CE 192 


vi 


1 
2 


LIST OF TABLES 


Шлгогпсегеге Це ит М аНнДайоп Кипз..... n ........................ 
Parameters Used in Orbiting Satellite Application ............................ 


» © N к 


сл 


10 


11 


1:2 


13 


14 


LIST OF FIGURES 


Left Surface of Rectangular Block Used їп е моае 7777777 5 
Flow Chart for Program EXPLICIT mnnn ОО 10 
Fiow Chart юг Ргоргат ВЕІАМ.. 2... 1.77 19 
Rectangular Model Showing Boundary Conditions Used in 

Validation Küns... eden ee RUE IT 21 
Flow Chart for Program VALID ИИ шшш 23 


Temperature Difference Versus Time for Time Increment of 
17.117 Seconds with Distance of 1 Meter Into Бјоек. ник 25 


Temperature Difference Versus Time for Time Increment of 
4.279 Seconds with Distance of 1 Meter Into Block......................... 26 


Temperature Difference Versus Time for Time Increment of 
17.117 Seconds with Distance of 0.4 Meters Into Block ................ 2m 


Temperature Difference Versus Time for Time Increment of 
4.279 Seconds with Distance of 0.4 Meters Into Block................... 28 


Temperature Difference Versus Time for Time Increment of 
17.117 Seconds on Lett Surlace of BloCle ЕИН 29 


Temperature Difference Versus Time for Time Increment of 
4.279 Seconds on Left Surface of Block c s 30 


Temperature Difference Versus Time for Time Increments of 
17.117 and 120 Seconds with Distance of 1 Meter Into Block....31 


Temperature Difference Versus Time for Time Increments 
of 17.117 and 120 Seconds with Distance of 0.4 Meters 


Temperature Difference Versus Time for Time Increments 
of 17.117 and 120 Seconds on Left Surface of Block........................ 39 


15 
16 


ЕТ 


18 
19 
20 
21 


Ии ропен оп а 5аге!епп Сїгсшаг ОгБЬї........................................... 37 


Geometry Used in Determining Geometric Factor for Earth 


атоо. О Д ео ИО 38 
Geometry Used in Determining Geometric Factor for Solar 

ее 39 
Temperature Variation of the Top Surface Center Node................. 42 
Temperature Variation of the Top Surface Center Node................. 43 
ERO UR FOIE есананлам а 22 2222... ннн 49 
саара атое Иннаа. 58 


іх 


LIST OF SYMBOLS 


Description 


albedo coefficient 


coefficient associated with 
coefficient associated with 
coefficient associated with 
coefficient associated with 
coefficient associated with 
coefficient associated with 


Biot number 


specific heat of material 


coefficient associated with 


coefficient associated with 
coefficient associated with 


solar aspect coefficient 


ж 5 Е , š 
T |] in Brian's technique 


[mm 


мах in Brian’s technique 


n+] 


Т нь. in Brians teclhnrque 
Ти in Brian’s technique 


Tù in Brian's technique 


n+l 


Tik № Ва КЕБІН ШЕШЕ 


ж . . 9 . 
Так in Brian's technique 


Тык in Brian s technique 


n+l 


Taa in Briair осше 


sum of satellite altitude and radius of earth 


coefficient associated with T* simultaneous equations in 


Brian's technique 


coefficient associated with T** simultaneous equations in 


Brian’s technique 


a | n+l, 
coefficient associated with T ^ simultaneous equations in 


Brian’s technique 


ALBEDO 


COND 


CONV 


SOLAR 


ST 
FA 
FE 
Fo 


rate of energy transfer into a control volume from albedo 
flux 

rate of energy transfer into a control volume from 
conduction 

rate of energy transfer into a control volume from 
convection 

rate of energy transfer into a control volume from earth 
thermal radiation 

rate of energy transfer into a control volume from imposed 
flux 


rate of energy generation 

rate of energy generation per unit volume 

rate of energy transfer into a control volume 

rate of energy transfer out of a control volume 

rate of energy transfer into a control volume from radiation 
heat transfer 

rate of energy transfer into a control volume from solar 
flux 

rate of increase of stored within a control volume 
geometric factor for albedo flux 

geometric factor for earth thermal radiation 

Fourier number 

heat transfer coefficient 

x-coordinate of node 

y-coordinate of node 


z-coordinate of node thermal conductivity of material 


Tari 


period of revolution for earth satellite 


imposed heat flux 
Ax 
С. 


(%) 


ambient temperature 


temperature at previous time level n 


temperature at half-time step interval 
temperature at n + = time level in the x-direction 
temperature at n + 2 time level in the y-direction 
temperature at current time level 

КјЕО 

Коко 

2RjFo 

2R9Fo 

2Fo 

FoAx?/k 

solar absorptivity 

time increment 

distance increment in x-direction 

distance increment in y-direction 

distance increment in z-direction 

emissivity; absorptivity of earth flux 


gravitational constant 


Xii 


G 


OS 
ФЕ 


mass density 
Stefan-Boltzmann constant 
solar flux 


earth thermal radiation flux 


The subscripts below have the following meaning: 


У] 
Vo 
V3 


V4 
V5 
V6 


oor WN >= 


left face 
right face 
bottom face 
top face 
front face 
back face 


-2Fooe&1R1Ay/k 
-2Fooe2RiAy/k 
-2Fooe£3R2Az/k 
-2Foo£4R2Az/k 
-2ЕосеБАх/К 
-2FooegAx/k 
FoBi, R? 

FoBi, R 


FoBi 3 R 


le = nje = nje 


FoBi , Rj 
FoBis 
FoBig 


Ay „ 
2FoR 1 (4, + 6 Еда + 0 Свл1 + a aF 0s) 


xiii 


2FoR (a БЕК Фе CO, Cane Ome merci) 
2FoR , AZ (q', > Е,.Е,ф,%0,Свл,ф.% а,аЕ, 9.) 
2FoR 552 (9 + £ Е,ф. + 0 Сза, фа авифо) 
2Fo (q' ев фо СО В 


Ах ^" 
2Fo- (q + & F; 0, а Сва еф. + сав, Ф5) 


ACKNOWLEDGMENTS 


The author would like to express his deep gratitude to Assistant 
Professor Y. Joshi, his thesis advisor, for his guidance and advice in 


the endeavor. 












та ты Ж фи. РАЋ 

i к | | Pu Ша n n 
12 i ie imt и 
"вући wa му му 


P У си 





| We СӨР | d ' 
mar yt 
a T 








Vx 


I. INTRODUCTION 


A. STATEMENT OF PROBLEM 

Three-dimensional transient heat conduction computations 
involving finite difference techniques can be costly and time consum- 
ing. In the explicit technique, there is a maximum time increment 
that can be used before the solution becomes unstable [Ref. 1]. This 
time increment results from an upper limit on the allowable values of 
the Fourier number [Ref. 2]. The implicit technique, on the other 
hand, requires more complicated computations than the explicit 
technique due to matrix inversion operations [Ref. 2]. Even the stor- 
age capability of a modern computer can be taxed when dealing with 
storage and manipulation of very large matrices [Ref. 3]. Both the 
explicit and implicit techniques can become expensive when a large 
number of nodes is required [Ref. 4]. 

For two-dimensional transients, an alternative technique is the 
alternating direction implicit (ADI) method. Calculations at each time 
step require only the handling of tridiagonal matrices [Ref. 5]. A 
direct extension of the ADI technique to three-dimensional problems, 
however, has been found to have only conditional stability [Ref. 5]. A 
different algorithm by Brian is unconditionally stable and still requires 
the handling of only tridiagonal system of equations, thus avoiding the 


more complex matrix manipulation of the implicit approach [Ref. 3]. 


B. OBJECTIVES 
The objectives of this study were: 


1. To develop a FORTRAN computer program employing Brian’s 
technique for three-dimensional transient heat conduction. 


2. To validate the FORTRAN computer program using the explicit 
technique with identical geometric configuration, boundary, and 
initial conditions. Comparisons with closed-form solutions avail- 
able for a semi-infinite geometry were also made. 


3. To employ the above technique in the application of temperature 
computations involving an orbiting satellite. 


This study was also intended to be a basis for future satellite heat 
transfer analysis and computations involved with Project Orion, the 
Naval Postgraduate School’s multi-disciplinary space research 
program. 


All objectives were achieved. 


II. DEVELOPMENT OF MODEL 


A. MODEL CONDITIONS 

The model is a three-dimensional rectangular block with geomet- 
ric dimensions expressed in Cartesian coordinates. A numerical tem- 
perature solution only permits the determination of the temperature 
distribution at discrete points or nodes [Ref. 2]. The method of finite 
differences requires that a nodal equation be written for each node in 
the block. The number of nodes in any coordinate direction is based 
on the distance increment in the same direction. 

The nodal temperature equation for each node was obtained using 
a conservation of energy equation within a control volume. Incopeia 
and DeWitt [Ref. 2] present the general form of the conservation of 


energy equation as 


ы” Еол Қ Еовһ = Egy (1) 


The nodal temperature equations were developed using equation 1. In 
the following sections, an equation for a node on the left surface of the 
rectangular block away from corners and edges is derived for both the 
explicit and Brian techniques. Equations for all other nodes are 
derived in an identical fashion. The explicit technique equations are 
tabulated in Appendix A and the Brian technique equations are tabu- 


lated in Appendix B. 


B. EXPLICIT TECHNIQUE 
1. Derivation of Nodal Temperature Equations 
Figure 1 shows the left surface of the rectangular block with 
associated boundary conditions. Equation 1 is rewritten assuming that 


all the energy components are into the control volume as 
E T E = E (2) 


where the energy generation term is assumed to be zero. 
Using finite difference notation, the energy storage term is 
n+l n 
: Ау To Б 
Е = рО EE EE 


2 At (3) 


The energy inflow term, E,,, is, in general, composed of convection, 
conduction, radiation, and imposed heat flux terms. For the develop- 
ment of the explicit technique, the radiation term was neglected. In 
the explicit technique, all involved nodal temperatures on the left side 
of equation 2 are evaluated at the prior time step. The convection 


term is thus 


pese ^x Az (T. — Tuk) 


(4) 





x 


Figure 1. Left Surface of Rectangular Block Used in the Model 


The net conduction term will have contributions from all three coor- 


dinate directions and is 


| Ау s; Ер, i ik Лл m к в. 
ЖЩ = К о == kAxAz Ay 
- s я (5) 
ae Ay i ЕТІ n" EM 
di 2 А2 
The flux term will be 
Е них = q” Ax Az (6) 


Substituting equations 3 through 6 into equation 2 produces 


та п п п п 


кТ..,-2Т -T 
Ау ша 1.j.k i*1.j.k ijk јак ТВ 
k 5 Az де + kAx AZ Ay 
Ау 1 а ЕН] F T sm = о 
kax p 7 e kAxz(T. - T5.) (7) 


n+] n 
„ Ау До T 
+q Axz = РСАх 5 42 — 


The temperature at the next time level, n+l, is next 
expressed in terms of temperatures at the previous time level, n. To 


accomplish this, terms in equation 7 are еи 





Multiplying each term in equation 7 by yields 


m 


2 


(Пра ШЕТ 7 2 9 J ЄЭ (Тула 2 Тым) 


2 


2 
+125) (т а 2T,,,) a 
| 2 С n 
+ EAXAY(T- - 159 +9 E E T LONE Та) 


Equation 8 is simplified using Biot and Fourier numbers, 





which are written as Bax and respectively. The terms 
k 


kAt 
pCAx^ ' 
2 2 Д 
Ax Ax А Ах А ‚г. 2 
ta (=) , and Е are identified as Ri, Ro, and R,, respec 


tively. Equation 8 now becomes 


S(T tT, POER (T ar Tik) 


i-1.j.k Blk € í k ijel.k i,j,k 
г ік (тыр Т.а DIS (9) 


; n » А n+ n 
+ BiR (т. zum *q R я э (Тык -T^u) 


Rearranging terms, the temperature at the current time level is 


expressed as 


n+l n n 
Tit F ека -1,j.k ат л) + 2Е „К Ginny = ши 
ОВ (Ток FT. a 27 ци) (10) 

; A 
ШЕВ Ет) 2FR g => 


In a similar manner, the temperature equations for each node 
are derived. 
2. Restrictions on Time Step in the Explicit Technique 
In the derivation of the explicit technique equations, the 
Fourier number was expressed as 


kAt 


Fo = OC Ax? (11) 


The stability of the temperature solution is dependent on the time 
increment [Ref. 2]. Rearranging equation 10 to isolate the node of 


interest at both the current and previous time levels yields 


n+] п 
Тос Боа Agen а es ое с 1 DE) 





„АУ 
k (12) 


Е 
à 1 - 2Fo - 2FoR , - 2FoR, - 2FoBiR d 


P 


The coefficient associated with the previous time level at the 
node of interest must be greater than or equal to zero [Ref. 2]. This 


implies 


1 
1-2Fo(1+R,+R,+BiR*)20 


(13) 
Rearranging terms, it follow that, for stability, 
Fo < EL 
1 «R, «R, e Bin*) (1щ 


Limitations on Fo were derived for each different node. The maximum 
allowable Fo was chosen as the least of all these values. Since the 
Fourier number is dependent on At, Ax, and material properties, it 
follows from equation 11 that the maximum time increment will be 
dependent on the distance increment. As Ax decreases, the maximum 
time increment must also decrease to have equation 14 hold true 


[Ref. 1]. 


3. Numerical Calculations 
A computer program (program EXPLICIT) was written to 
evaluate the transient temperature distribution within the block using 
the explicit technique. Figure 2 is a flow-chart description of the 
program. A complete listing of program EXPLICIT is provided in 
Appendix C. 


C. BRIAN’S TECHNIQUE 

The Brian technique provides unconditional stability, allowing 
large time increments for transient three-dimensional computations 
[Ref. 3]. As mentioned earlier, the explicit and fully implicit tech- 
niques may become unsuitable for computations involving a large num- 
ber of nodes or long transient times. The explicit technique relies on 
small time increments for stability, thus increasing the amount of 
computer time required for calculations [Ref. 5]. The implicit tech- 
nique requires inverting large matrices that may create problems with 
storage capability and computer time [Ref. 4]. 

Another possible technique, the alternating direction implicit 
(ADI), is applicable to both two- and three-dimensional problems. 
While the ADI method is unconditionally stable for two-dimensional 
cases, instability occurs for large time increments in three-dimen- 
sional problems [Refs. 4, 5]. 

The Brian technique modifies the ADI method to provide uncon- 


ditional stability and minimize computer time [Ref. 3]. 
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Figure 2. Flow Chart for Program EXPLICIT 


1. Derivation of Nodal Temperature Equations 
Figure 1 shows the left surface of the rectangular block with 
associated boundary conditions. Equation 2 is again the starting point 
for nodal temperature equation derivation and is 
E + E = E (2) 


IN GEN ST 


where the energy generation term is again assumed to be zero. 

The nodal temperature equation derivation for Brian's tech- 
nique uses half-time steps for each coordinate direction. This differs 
from the explicit technique, in which a whole-time increment is used 
and spatial coordinate direction is not considered. 

In Brian's technique, the first half-time level, n+, treats 
temperatures in the x direction as unknown with temperatures in the 
y and z directions known from the previous time level, n. Solving a 
tridiagonal system of simultaneous equations, the temperatures at the 


п + > level are determined. 


Temperatures in the y direction are then determined using 
x-direction temperatures at the n+ time level and z-direction tem- 
peratures at the M time level. Similarly, temperatures in the 
z-direction at the current n+1 time level are determined using n + 5 
temperatures from the x and y directions. 

The above procedure is demonstrated for a node on the left 


surface away from corners and edges. 
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The energy storage term is 


Ay Ж 
E, = pC Ax Az a (15) 


The energy inflow term, E,,, is composed of convection, conduction, 


radiation, and any imposed external heat fluxes. As a first step, only 
the nodal temperatures associated with heat conduction and energy 


storage in the x-direction are evaluated at the n + > time level. The 


convection term is thus 


ЕШТЕ hAxAz( T. -Ті.) 


(16) 


The net conduction term will have contributions from all three coor- 


dinate directions and is 


Ay к + TE = 2M T тч ds 
(17) 
И с АУ Tua pu ee. Ж 2T jk 
2 AZ 
The radiation term, assuming extensive surroundings, will be 
: 4 n* 
E др оАхАа (Те - ТО) (18) 


For applications involving an orbiting spacecraft, additional input flux 
terms need to be considered. These include the solar flux, the 


infrared emission from the earth, and the reflected component of the 


12 


solar flux from the earth’s atmosphere, the albedo. The solar flux term 


will be 


Е РЕ СО. (19) 


5 
The reflected solar flux (albedo) term will be 


E = АхА2Ғ, о аф. (2 0) 


ALBEDO 
The earth flux term will be 


E Кени Си ЕТЕ (21) 


Е 


The imposed heat flux term will be 
Е ых = АХА (22) 


Substituting equations (15) through (22) into equation 2 yields 


п+; па n+ 


Ay ШЕТ ш D 7 2m eat m To 
MI пина К Ay 
ду ЕЕ en 8 2 5 14 
+ КАх-2- АБ + oeAxAzÍ TŠ = Eus 


+ Ax AzaC „Ф„+ АхАФЕ, aad, + AXAZF, EQ, -- AXAZQ' 
(23) 
Ay d d 
) Б oe a ae 
2 


п 


as (P то, 
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The temperatures at time level n +> are next expressed in 
terms of temperatures at the previous time level, n. To accomplish 


this, terms in equation 23 are rearranged. 
Multiplying each term in equation 23 by [m yields 





1 1 1 
1 nt; ntz HÀ 
fp? yp? -от тї) (E) (Tha k пе) 





Ax Y n n ‚ бєАх° 
T iU, Ше meum 2 jk "Ұлу С -Т у) 





А2 
ОС ОЗЕ oF, aAx? Е. РАХ Ах? , 
—kAy S Ae C XB Т (24) 
hAx? _ 2pC Ax* | т J 

ш Клу f E - (T.- a ОКАТ bre c Dag 


Equation 24 is simplified using Biot and Fourier numbers, 


: ! h kAt . 
which are written as к Ах апа 2pCAx? ' respectively. Equation 24 
now becomes 


1 1 


: n+> ny Ч 
т ОГ ото P E 


2 i-1.j.k "e 
+IR,(T' tT eyes 2T) * SERAy(T5- T7,,) 

aC ., R, Ay aF, aR áy F eR Ay (25) 
t= ре Оз = 


R, Ay ; 1 | ntz n ) 
k а +В, В (т. = TG ee 2 Fo Tor Tg 








+ 
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Rearranging terms, the temperatures at the first half-time step level 


in the x direction are 


- вот? ШЕР 2Fo)T,, TT TESTE aos 


1-1,}.К i,j,k i.J+1,k 


+2FoR APR ) * FoR, (и Due 21 0), 


КЕ та 


(26) 








«С, ШЕ 
+2Fo Rok 9,7 2Fo R Ау k 








k 


Hn 


+2Fo R Ay + T- + 2FoBi Rae Eu) 


The next step is to express the temperatures in the y-direc- 
tion at time level n к? using x-direction temperatures determined at 
time level n + from equation 26 and z-direction temperatures at 


time level n. For greater clarity, the x-direction temperatures at time 


+— 
2 


level n + >, T ^ ,arerewritten as T*. Equation 26 now becomes 


n 


—FoT*^ 


ER t (1 + шо – Еот =i 


i*1.j.k i.j+1.k 


ын о xm 


+ 2Еок (Г ) +FoR, шаси ША = з) 


(27) 








aCe 
СЕ 4 
+2FoR Ay (т: - T "m = 





/? 


БЕЛЕ q я с n 
ОВ, ЛУ, 0, + 2FoR,Ay__ +2FoBi Ri(T.-T,,,) 


The derivation for the y-direction nodal temperature equation is simi- 
lar to the x-direction derivation. The y-direction equation for the same 


node is 
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+2 


(1+ 2FoR JT. - 2 ҒоК as. " =T jk 


n 


+T 


+FO(T t, pa "Та 7 2T f ) FOR (Tha Тоа 


i-1.].k 1+1,}.К Е 


4 





+ 2Fo R,Ay S(T* сте |$ 2Fo R,A 





-- 5 











aF, a E 
+2Fo R,Ay— Ф, k 


и 


+ 2Fo В, Аут + 2FoBi R\(T.-T,,;) 


КЕ" 


1 
п += 


2 ° 
, are rewritten 


The y-direction temperatures at time level n t, T 


as T**. Equation 28 now becomes 


(1+ 2FoR ,)T** -2FoR, Tt*,,=T! + Fo(T#,,, +T* 


+ Бо ТИ ое k (ri у 
oC oF, a МЕНЕ 
*2Fo R, Ay—- -6,* 2Fo R, Ay —— 6, * 2Fo R;Ay -— $c 


+ 2Ео К ays + 2FoBi R, ic ть) 


i-1 j.k ак -2 


-2T x) 


(28) 


Tisk ) 


(29) 


The expression for the z-direction temperatures at time level 


n+l is dependent on the x-direction and y-direction temperatures at 


time level n+ >, T* and T**, respectively. The z-direction nodal tem- 


perature equation is 
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EEOR,TO e or(I-2FoRj)T;; —FoR;T,,- 





íj. K=] јок 
Tijk *tFo(T?, i alis E 2T )*2FoR, ut bugs ) 
OE «С, 
Шкок лу (Т. – Та )*2FoR,Ay— 6, 
aF, a ЕЕ 
+ 2Е0о В Ay - ф. + 2Ео R,Ay— 6; е) 


и 


+2Fo R Ay + 2РОБ | (Ta -THA 


In a similar manner, the temperature equations for each node 
are derived. We note that during traverse along each coordinate direc- 
tion in equations 26, 28, and 30, a tridiagonal matrix results for the 
solution of temperature. 

2. Solution of the Nodal Temperature Equations 
The x-direction nodal temperature equations for the rectan- 


gular block are of the form 


* * * = 
АЛТ. Ш А uei eee ШЕП, (31) 
where the coefficients Aij» -e Dije are known from the previous 


time step. 


Each row of T* equations parallel to the x-axis is a tridiagonal 


matrix of coefficients Ai,,,, Bi,,;,, and Ci,,,. It was solved using the 


standard Tridiagonal Matrix Algorithm (TDMA) [Ref. 5]. After each row 
of T* equations is solved, the y-direction equations will be solved [Ref. 


3]. These are of the type 
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= DD 


As Tet, + Bo COMME M (3 


ж 
ЈАК i,j,k ij+1.k 


As in the x-direction case, each row of T** equations parallel to the 
y-axis is solved using a tridiagonal matrix solver. The z-direction equa- 


tions are of the form 


n+l n+l n+l 
Азд зук + Е Зак пак + C3, Тадич = Ds 


dk (33) 


The solution of the T?” equations is accomplished in a similar man- 
ner as used for both the x- and y-directions [Ref. 3]. 
3. Numerical Calculations 
A computer program (program (BRIAN) was written to evalu- 
ate the transient temperature distribution within the block using the 
Brian technique outlined in the previous sections. Figure 3 is a flow 
chart description of the program. A complete listing of program 


BRIAN is provided in Appendix D. 
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Enter Parameters: 

e Material Properties 

*« Geometric Dimensions 
* Other Constants 


пицайге: 

* Time 

* Printout Counter 

* Nodal Temperatures 





Increment: 
e Time and Counter 


Space Appltcation? 





Determine Solar and 
Albedo Fluxes Based 


on Sunlight and 
Darkness Times 


Set A, B, C, D Coeffi- 
cients in x direction 
Solve for Тук Using 
Tridiagonal Matrix 


Set A, B, C, D Coeffi- 
cients in y direction 


Solve for Tii Using 
Tridiagonal Matrix 


Set A, B, C, D Coeff- 
cients in z direction 


Solve for To Using 
Tridiagonal Matrix 


No Yes 
Validation Desired? 
Printout Interval No 
Reached? 
Print Temperature Differences and 
Reset Counter to Zero 














Print Temperatures and 
Reset Counter to Zero 


No 
Problem Time Reached? 


Yes 


CPU Time Calculated 


Figure 3. Flow Chart for Program BRIAN 
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III. VALIDATION 


A. SELECTED GEOMETRY AND CONDITIONS 

Both the Brian's technique and explicit technique programs 
required validation to ensure accuracy in the temperature solutions. 
Identical geometric and boundary conditions were used for these 
techniques in the validations. The left face had a constant surface heat 
flux with all other surfaces treated as adiabatic. Figure 4 shows the 
geometric and boundary conditions for the validation model. Table 1 
lists the parameters used in the validation runs. 

The Brian's technique and explicit technique programs were 
written to allow the user to conduct a validation test. The test involved 
taking the temperature differences at identical nodes in each program 
and then comparing each node's temperature difference on a graph. 

The finite difference solutions for the rectangular block were also 
compared with an analytical solution for a semi-infinite solid exposed 
to a uniform flux at its surface. The closed-form solution would be a 
reasonable approximation to the actual solution for small times when 
interior nodes are not influenced by the lateral boundary surfaces of 
the rectangular block [Ref. 2]. 

Incropera and DeWitt [Ref. 2] tabulated the closed-form solution 
for the case of constant surface heat flux in terms of complementary 


error functions. These were evaluated in the present study based on 
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x 
Figure 4. Rectangular Model Showing Boundary Conditions 
Used in Validation Runs 


TABLE 1 
PARAMETERS USED IN VALIDATION RUNS 





875 J/kg-K 
275K 









the polynomial approximations provided by Abramowitz and Stegun 


[Ref. 5]. 
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B. NUMERICAL COMPUTATIONS 

A computer program was written to evaluate the closed-form 
solution for a semi-infinite solid with a constant surface heat flux. The 
fifth-order polynomial approximations to the complementary error 
function have an error less than 1.5 x 1077 [Ref. 5]. 

Figure 5 is a flow-chart description of the closed-form solution 
program VALID. The parameters entered are identical to those used in 
programs EXPLICIT and BRIAN. A complete listing of program VALID 
is provided in Appendix E. Table 1 lists the parameters used for the 
validation runs. 

The validation period for programs VALID, EXPLICIT, and BRIAN 
was 3,600 seconds (one hour). The distance inside the semi-infinite 
solid was based on pre-selected nodes in the rectangular model. Three 
distances were used in the semi-infinite solid, corresponding to three 
selected nodes in the rectangular block. 

In the first validation run, each coordinate direction had 11 
nodes. This corresponded to an explicit technique time increment of 
17.117 seconds, which was also selected as the time increment for 
programs BRIAN and VALID. In the second validation run, the number 
of nodes was increased to 21 in each coordinate direction. The time 
increment used in this run was 4.279 seconds. A third validation run 
was conducted using 11 nodes in each coordinate direction. The 
explicit technique and program VALID used a time increment of 


17.117 seconds. The time increment for program BRIAN was 
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Enter Parameters: 

e Flux 

¢ Thermal Conductivity 
e Material Density 

e Specific Heat 

e Distance 

e Distance Increment 
e Printout Interval 


Present Constants: 
ел 
e Polynomial Coefficients 


Initialize: 
° Time 
e Printout Counter 


Increment: 
e Time 
e Printout Counter 


Calculate: 

« Temp. Difference 
Using Closed-Form 
Solution 


Printout Interval 
Reached? 


Yes 





Print Temp. Difference 


Reset Printout Counter 
to Zero 





Validation Period we 
Reached? 
Yes 


End 


Figure 5. Flow Chart for Program VALID 
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increased to 120 seconds. The increased time increment was used to 
determine whether the program BRIAN temperature differences 
would differ significantly from those previously obtained with the 


17.117-second time increment. 


C. VALIDATION COMPUTATIONS 

Graphical representations of the temperature differences experi- 
enced over the validation period are presented in Figures 6 through 
14. 

At a distance of one meter into the block (Figures 6 and 7), the 
temperature differences calculated in EXPLICIT and BRIAN diverge 
significantly from VALID after a period of 616 seconds. This is to be 
expected as the effects of the constant surface heat flux will provide 
higher temperatures in the rectangular block than in the semi-infinite 
solid. The rectangular block can thus be modeled as a semi-infinite 
solid for only about 10 minutes at a distance of one meter. 

At a distance of 0.4 meters into the solid and on the surface 
(Figures 8 through 11), both EXPLICIT AND BRIAN temperature dif- 
ferences were in close agreement with the closed-form temperature 
differences for almost the entire validation period of one hour. 

Increasing the number of nodes and, consequently, decreasing 
the time increment had virtually no effect on temperature differences, 
as can be seen in Figures 6 through 11. 

Increasing the program BRIAN time increment to 120 seconds 


did not degrade the accuracy of the solution. In the three cases 
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3600 
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TIME (SECONDS) 


600 
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Figure 6. Temperature Difference Versus Time for Time Increment 
of 17.117 Seconds with Distance of 1 Meter Into Block 
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Figure 7. Temperature Difference Versus Time for Time Increment 


of 4.279 Seconds with Distance of 1 Meter Into Block 
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Figure 8. Temperature Difference Versus Time for Time Increment 
of 17.117 Seconds with Distance of 0.4 Meters Into Block 
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Figure 9. Temperature Difference Versus Time for Time Increment 
of 4.279 Seconds with Distance of 0.4 Meters Into Block 
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Figure 10. Temperature Difference Versus Time for Time 
Increment of 17.117 Seconds on Left Surface of Block 
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Figure 11. Temperature Difference Versus Time for Time 
Increment of 4.279 Seconds on Left Surface of Block 
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Figure 12. Temperature Difference Versus Time for Time Increments 
of 17.117 and 120 Seconds with Distance of 1 Meter Into Block 
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Figure 13. Temperature Difference Versus Time for Time Increments 
of 17.117 and 120 Seconds with Distance of 0.4 Meters Into Block 
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Figure 14. Temperature Difference Versus Time for Time 
Increments of 17.117 and 120 Seconds on Left Surface of Block 


33 


involving distances of 1.0 and 0.4 meters into the surface and on the 
surface, the temperature differences did not vary by more than 0.045 
percent from those obtained using a 17.117-second time increment. 
Figures 12 through 14 show the 120-second time increment plot of 
temperature difference versus time being in close agreement to the 
corresponding 17.117-second time increment plots. 

Using a time increment of 17.117 seconds, the maximum error 
was found to be within 0.019 percent for an 11 x 11 x 11 grid. The 


maximum error was calculated using 


9% ERROR = —coMeuren — Texact 100% 

T EXACT (34) 
where TCOMPUTED is the Brian’s technique computed temperature and 
ТЕХАСТ is the closed-form temperature. Sampling points were taken 
at 600-second intervals. 

In a similar manner, the maximum error was calculated for a time 
increment of 4.279 seconds. The maximum error was found to be 
0.021 percent for a 21 x 21 x 21 grid. 

A time increment of 120 seconds was above the maximum allow- 
able time step from stability considerations for the explicit technique. 
Brian’s technique thus allowed time steps above the limitation of the 
explicit technique with very little compromise in solution accuracy. 

Based on the validation test results, program BRIAN can be con- 


sidered correct and validated for further use. 
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ГУ. APPLICATIONS 


A. GENERAL 
The Brian technique was used to solve for the temperature distri- 
bution in an application example. The example involves the transient 


thermal response of a component on board an orbiting satellite. 


B. ORBITING SATELLITE COMPONENT 
1. Geometry and Initial Conditions 
The transient thermal response of an orbiting spacecraft 
component will be investigated. The spacecraft component is in a non- 
geosynchronous circular orbit at an altitude of 1,609.364 kilometers 
(1,000 miles). The equation for the orbital period is listed in Agrawal 
[Ref. 7] as 


D (35) 


where u is the gravitational constant of the earth and d is the sum of 
the satellite altitude and radius of the earth. Eisele and Nichols [Ref. 8] 
list the values of u. and earth radius as 398,603.2 km?/s? and 
6,378.165 kilometers, respectively. 

The period of revolution in this case is 7,104.4 seconds or 
118.40 minutes. For a satellite in a circular orbit at an altitude of 100 
miles, Stevenson and Grafton [Ref. 9] listed the sun exposure time as 


72 percent of the orbital time. During one revolution around the earth, 
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the satellite will have 85.25 minutes of sun exposure and spend 33.15 
minutes in the earth’s shadow. For the example application, the satel- 
lite is assumed to start the orbital period at the instant it becomes 
exposed to the sun. 

Single spin stabilization is assumed, thus allowing the 
spacecraft component to spin about a principal moment of inertia axis 
[Ref. 7]. The spacecraft component spins about the z-axis, as depicted 
in Figure 15. With single spin stabilization, it is assumed that the 
component analyzed here will have all four sides continuously exposed 
to the sun with the top and bottom surfaces not receiving solar flux. 
The fraction of solar flux that is incident on the satellite surface is 
usually defined as the solar aspect coefficient [Ref. 7]. Thus, for this 
example, the solar aspect coefficients are assumed to be 1.0 for the 
side surfaces and 0.0 for the top and bottom surfaces. 

Stevenson and Grafton [Ref. 9] compiled the geometric fac- 
tors for both earth radiation and reflected solar radiation (albedo) on a 
flat plate. The geometric factor for earth radiation depends on the sat- 
ellite’s altitude and the altitude angle between the normal to the rect- 
angular component and the vertical to the satellite from the earth [Ref. 
9]. Figure 16 depicts the geometry involved. The geometric factor for a 
rectangular plate at an altitude of 1,000 miles (1,609.364 kilometers) 
and altitude angle of 60 degrees is listed as 0.3174 [Ref. 9]. 
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COMPONENT 


CIRCULAR ORBIT 


Figure 15. Component on a Satellite in Circular Orbit 
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ALTITUDE 





EARTH 


Figure 16. Geometry Used in Determining Geometric Factor 
for Earth Thermal Radiation 


The geometric factor for reflected solar radiation is depen- 
dent on the satellite’s altitude, the altitude angle, the sun angle, and 
the angle of axis rotation [Ref. 7]. Figure 17 depicts the geometry 
involved. The geometric factor is 0.0456 for a satellite at 1,000 miles 


altitude, altitude angle of 60 degrees, sun angle of 90 degrees, and 


rotation angle of O degrees [Ref. 9]. 
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RECTANGULAR PLATE 


ALTITUDE 





1. ALITTUDE ANGLE 
2. SUN ANGLE 
3. ROTATION ANGLE 


SUN 


EARTH 


Figure 17. Geometry Used in Determining Geometric Factor 
for Solar Reflected Radiation 
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ТАВГЕ 2 
PARAMETERS USED IN ORBITING SATELLITE APPLICATION 


Dimensions 0.15 m x 0.10 m x 0.015 m 
Material aluminum alloy 2024-T6 
* density 2,770 kg/m? 

e thermal conductivity 186 W/m-K 
e specific heat 1,000 J/kg-K 
Exterior Coating black paint 
e solar absorptivity ШТ 

° emissivity 092 
Imposed Flux 

° bottom surface 50,000 W/m2 
Orbital Data 

e altitude 1,609.364 km (1,000 mi) 
° period of rotation 118.40 min 
° exposure to sun 85.25 min 

° exposure to darkness 33.15 тїп 
Solar Aspect Coefficients 

* side surfaces 1.0 

e top and bottom surfaces 0.0 
External Flux 

e solar flux 1,353 W/m2 

e earth radiation 237 W/m2 
e albedo coefficient 0.3 
Geometric Factors 

e earth radiation 0.3174 

e albedo flux 0.0456 
Initial Data 

е heat transfer coefficients 0.0 W/m2-K 
e ambient temperature ок 

e initial nodal temperature 500 K 
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Agrawal [Ref. 7] listed the mean annual solar flux, mean 
annual earth radiation, and mean annual value of albedo coefficient as 
1,353 W/m?, 237 W/m?, and 0.30, respectively. These value are used 
in the example. 

The material coating on the component's exterior determines 
the solar absorptivity and emissivity. The component is assumed to be 
coated with black paint, thus providing a solar absorptivity of 0.97 and 
an emissivity of 0.92 [Ref. 10]. 

The component is rectangular with dimensions of 0.25 m x 
0.10 m x 0.015 m and is made of aluminum alloy 2024-T6. Incropeia 
and DeWitt [Ref. 2] list the values for density, thermal conductivity, 
and specific heat. 

The component is mounted on the exterior of the satellite. 
The bottom surface is attached to the satellite with all other surfaces 
exposed to the space environment. The flux on the bottom surface is 
assumed to be 50,000 W/m2. 

Table 2 lists the properties and other parameters used in the 
example. 

2. Numerical Calculations 

A time increment of 120 seconds (2 minutes) was used in 
program BRIAN. The program covered a period of 120,000 seconds 
(33.33 hours), during which the spacecraft component conducted 
16.9 orbits. The temperatures of top surface node were plotted versus 


time, as shown in Figures 18 and 19. The temperature versus time 
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Figure 18. Temperature Variation of the Top Surface Center Node 
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Figure 19. Temperature Variation of the Top Surface Center Node 
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plots in both figures demonstrate an imperfect sinusoid. Steady peri- 
odic response occurs following orbit. The peak temperatures during 
the third orbit vary from the second orbit peak temperatures by 0.27 
percent. The largest temperature difference in successive peak tem- 
peratures is 0.27 percent and occurs between the second and third 
orbits. 

The oscillations in temperature distribution may be attributed 
to overlaps in the sunlight and eclipse periods. The sunlight exposure 
time for the orbit is 85.25 minutes. A time increment of 120 seconds 
(2 minutes) would allow solar flux to be seen by the satellite compo- 
nent for an additional 1.25 minutes on the first orbit. Thus, 1.25 min- 
utes of the first eclipse period would be missed at the beginning of the 
eclipse period. At the end of the eclipse period, the component is in 
darkness an extra 1.60 minutes. These overlaps occur throughout the 


running of the program with the maximum overlap being 2 minutes. 


a4 


V. RESULTS 


Brian’s technique was validated for four of the six rectangular 
surfaces. For cases involving flux or convective boundary conditions on 
the front and back surfaces of the model, temperature distributions 
did not follow the pattern obtained in previous validation cases. This 


anomaly requires further investigation. 
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VI. CONCLUSIONS 


Brian’s technique was validated using the explicit technique and 
the closed form solution for a semi-infinite solid. For the same time 
increment, the explicit technique uses six times less CPU time. How- 
ever, as the time increment was increased, the CPU time used in 
Brian's technique became less than that used in the explicit technique 
for time steps larger than 110 seconds. Brian’s technique used time 
increments that would lead to instability if used in the explicit 


technique. 
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VII. RECOMMENDATIONS 


While developing the computer code for Brian’s technique, mate- 
rial properties were assumed to be independent of temperature. A 
more precise code would make thermal conductivity a function of 
temperature. 

Another improvement for the computer code when used for 
spacecraft thermal analysis computations involves the application of 
orbital mechanics. This would permit the use of non-circular orbits in 
the code. 

This study was restricted to the use of Cartesian coordinates in 
Brian’s technique. Further study is required in applying Brian’s tech- 


nique in cylindrical and spherical coordinates. 
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Figure 29. Geometry of Rectangular Model 
3. Bottom Face 


n+l 
T ра ОФ А рани I ил) 


i,j,k 


n 


кк (ш КО, 21 ) 


11-1. Ил „ЈК 


+2FoR (Щи pe: ) 


i.j.k+1 1.j.k 


4. Top Face 


n+l 


Гв ЕЕ ЕШ Ете м 2553) 


+ FoR er k pal k E 9 


+2FoR («е | =) 


49 


5. Front Face 


n+l 
T ШТ 2Fo(T, =T. 


ik ТРЕ ijk 
+ FoR (Tiat T арии 


+ FoR E TEE к) 


6. Back Face 


Thuc Dp 2B MN 


i.j.k Најк i.j.k 
n 
* FoR ВЕ а да к) 


+ FoR КЕ ae BT) 


C. PERIMETER NODE EQUATIONS 


1. Perimeter 1 


п+1 


и ЕО nie 


n 


* FoR M CU: а За) 


+ 2 ЕоВ. дате zum) 


50 


2. Perimeter 2 


n+l n 
T, = Е т т 2Fo(T, -1.jk pb) 
* FoR nec k E 52 Due) 
+ 2 ЕоВ ава SE) 


З. Perimeter 3 
n+l n 
M HORE) 


1,3. К Најк 1.j.k 


* FoR (Тұма = 2T. e) 


+2FoR = 126) 
4. Perimeter 4 
п+1 п 
1 = jue + 2Fo(T, -1 j.k i) 


* FoR (Disa D ол.) 


+2FoR e c =) 
5. Perimeter 5 
Тук“ Ту + 2Е0(Тълук Tiga) 
+2FoR (т к s) 


* FoR ГЕ Тее 2T, x) 


51 


6. Perimeter 6 


n+l n 
1. k Та 2Е0(Т j.k 25: 


+2FoR Це 4 -T jk 


+ FoR a ен этү) 


7. Perimeter 7 


п+1 п 
JE TOR 2Fo(T;,_ -l1.jk YI 


+2FoR or р ) 


ijtl.k ijk 


+ FoR T Тума” 27 а) 


* 2FoBiR (T. - T7,,) * 2Foq’R = 
8. Perimeter 8 
n+l n 
dS ver ur 2Fo(T,, jk X 
+2FoR pure Р Tm) 
+ FoR ое Tor 2T, x) 
ы. AM 
В БОБИ ү 1.) “ 2Еод“В, 


52 


9. Perimeter 9 
n+l 
Dee ОИЕ гт eu 2Т к) 


+2FoR, (T ax- Tia) 


1,3-1. i,j,k 


+2 FoR (@ а жй 


10. Perimeter 10 
ue cun ROS а 2T а) 
on (1 и) 
+2FoR “ШЕ ы 1 а) 


11. Perimeter 11 


n+] 
Ти l s Fo; -1.j.k TIUS 2e 


и оса пе) 
+ 2FoR I реа 1 =) 
+ 2FoBiR i. >т, 5) Е 2ғо к, 5У 


53 


D. 


12. Perimeter 12 


п+1 
T =a BO TT Ton ee) 


i,j,k 


*2FoR((T,,, ,- T, 


ijt] .k 1,3. К 


pa 


+ 2 БОК ч jk пак) 


) + 2Foq’R, EM 


+ 2 КОВК iei -T. k 


i,j,k 


CORNER NODE EQUATIONS 


1. Corner 1 


TOL ARS ZEN VN 


ij. k i+] j,k ijk 


+2FoR (во T ) 


i,j+1.k i.j.k 


+ 2FoR ти Ера ) 


ij.k+] ij. K 


A 
+ 2Foq’R, = 


2. Corner 2 


n+l 
i,j, k 


Te ge cho S T 


i-l.j.k i.j.k 


+2FoR (T, 2 ра ) 


i.j+l,k i,j.k 
+2FoR ns - 1. ) 


ij. k+] i,j,k 


A 
+ 2Foq’R,=— 


54 


3. Corner 3 


= = T k + 2Fo(T,, m xi) 
+ 2ЕоК Ten k E) 
+ 2FoR м WS, 
+ оа 
4. Corner 4 
Tv = E 2Fo(T, IE СТ за) 
+ 2ЕоК eee k ze 
4 2FoR ТА 1 21 а) 
+ 2Fod’R, 
5. Corner 5 
D = j. k + 2Fo(T,, m Uo) 
+ 2FoR (Ж: 1k E) 
* 2FoR, (T7, 44 7 T; 4) 


+ 2FoBiR ха ВВ к) 


55 


6. Corner 6 


То. = Е -T sx) 


ЛЕ 


ИР. 


јак i.j.k 


+2FoR (Cnm 


+2FoR aa ен -Тіш) 


+ 2FoBiR ien zum к) 


7. Corner 7 


T T. HOM MS T) 
ОРОВ а) 


432 FOR MM и еа 


+ 2FoBiR та – Тил) 


8. Согпег 8 


qd = TM ‚+ 2Fo(T,. NU = lee) 


і, ј. К 


„ОКОВ (ТЕ и = Ty 


i,j-1,k i,j.k 


+2FoR Е Тук) 


+ 2FoBiR tens Tx) 


56 





А. 


APPENDIX В 
BRIAN'S TECHNIQUE NODE EQUATIONS 


X-DIRECTION EQUATIONS 


1. Corner 1 


ЕО ОЕ От, нь т, к) 


мект Тарк 
+U нес ІЗ БП» (T. - NS EVO (Е non) 
У (То - T4,,) € W, € W, W, +X (TR, -TÀ) 

+X АШ Да 3) +X, СЕ SuSE Y SC 


2. Perimeter 12 


~FoT* + (1 + 2Fo)T * 


i-1.j.k 


ijk КОТ үк = 


l U (T ae ili ANE Гу) 


-VQ(T., - T ux) * Vs Ge = Тик) +, +W, 


wait. a |) +х, pros По СВОИ Ee 


57/7 


LEFT RIGHT 





Figure 30. Geometry of Rectangular Model 
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=> НОШЊА + (1 +2FoR Ns = Ор (ds 1 PI m 2080) 


i,j-l.k 


ЕП (Тыл ана в iu) aus (ies T) +W. 
HEN dat и To.) +X (Th. то Ва 
25. Согпег 4 
(1+ 2FoR,)T** - 2FoR Ts = = tye Uh ai 
ыша алм E уы bar T 9) 77 a (T=; EM 
£VL(T., - T.) &W, e WW, X (йет) 
xX (TA, - TL)*X (T4, - пен) а Е 


26. Регітетег 2 


n 


-Коқ + (Е РОВ КЕШЕ ЕСЕН ҮГІТ 
Шы к ды Tie oe (Тик p G В) 


LV (T. По ss: (т, - ща 4) 


С 


+ X (т - т) eI 
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27. Corner 8 


р вов т = Е ТИ (Тоа Ти) 


i,j-1,k i.j.k i,j. k-1 i,j,k 


UT T^, км, (Т-, (ТТ) 


+V,(T.,-Tij.) +W,+W,+W, +X (тест) 
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ex ти, тъ) ек (тт) ие 


C. Z-DIRECTION EQUATIONS 


1. Corner 1l 


(1+ 2FoR, )T tix 7 2FoR T ENDO ш CO. go i RES ) 
+U EI ique si eo П ЕИ у>У,(Т-; Так ) 


+V,(T-,-T*, )+W,+W,+W,+X,(T*" та) 
meta, та учи (та ставо, 


i.j.k 


2. Perimeter 8 


-FoR Ti p +(1+2FoR,)Ti,,-FoR,T,),=T** 
U (Tuna 5 Т оа ЛЕ ВТ ) 


-У.(Т..-Т%Қ ум, ич (тии ста) 


КА: 1.j.k 


+X Cee: - T) +У Боја 


го 





3. Corner 3 


n+] n+l 
-2FoR,T, 4 t T 2F6FR ЈЕ у Ра E 


tU (T. и 


no )+V (T= T ) + V, (T=, g ) 


ИМ (и )+W,+W,+Ww,+X,(T*% -TÀ) 


Rx (aus E .) +X (Шей - Tš.) US | 


1.j. k 


4. Perimeter 12 


(I+ 2FOR,)T Si, -2POR Tia =THE * UST - TA 


1,3. К 


+Fo(T#,,, +T? 


i+1,j,k - 2 


4 
+W,+W, +X,(T** 


5. Left Face 


T5, )*V(T--T*S )*V4(T., 


- Т. ‚) + Х (mess -т*,) ТУ, Ев, x 


n+l n+l 
-2FoR Тула + (1+ 260R) ia Doe ee ee 


ео № 477. +T 


+W,+W,+X,( 


ж 


wijk 2 


жж 
тен 


T ) +V (T.-, Мы ) UE (Т... 


-Т. ) (OM -T-,)+Y,Ec,,, 
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k 


) 


) 


ж ж 
= ie 


р 


ж 
1j.k 


) 


) 


6. Perimeter 11 


n+l n+l 
-2 FoR НЕ и. met (1 * 2 ВОК 5) Тук = Í 
Та mu -2 Í )+V,(T-, 


+W,+W, +X, (T -TLn)eX((TAA 


7. Corner 2 


PE FOR) T FoR, T „„с Та“ 


ijk ij kl ЈК 


+U Sere m Er. )+V (Te, = 


i,j,k 


urs plo ya Wwe BW. N +X er 


B (Теа E Tx ) 
к )+V,(T-, Ei уь 
-Ti,)*Y;Eo,,, 

+ U (Та 7 Тук ) 
) * Va(T4 SB ) 


+X,(Tr™ -Tl,)+X(T#"% т) +Ү,Е 


8. Perimeter 7 


- T) 


-FoR JM ,* * 2FoR;)T7, - FoR AW NS M s 
+ 9 (T£ x Hi аа D. T = 1 

+V (T-, EL. RSS -T** DEM 
+X (ынты - T4) То ата - TŠ.) ШҮ Еш 
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) 


9. Corner 4 


-2FOR „Тура (DS 2FORTE TO NM M NE Wb 


ij. k-1 ijk ijk јак ij.k 
ып (тына YS TD EAD LE 


i,j,k i,j. k 


+V (Te -TIA )+W, +W, +W, +X, (T: сте) 
+X nen -Ti,)«X, (Ghee - TŠ.) „У Ес a 


10. Perimeter 3 


(1: 2ЭБок, УК = 2FoRet. 3. = т 


ГТК i,j, k+1 ijk 
ж ж ва жж _ TT * 
+U (T + ШТ 27,5% ука (е p 


+V,(T-, zd у+У (Т. = )+W,+W, 


+ X (T - Рен CM - Тие) i. INE 


i.j.k 
11. Back Face 


-FoR Ti" , * (1 * 2FoR,)T7^, – БОВ „Таја 27: 


i,j. k-1 ij.k ЈК ij.k 
жж m ж ж ж _ тж 
+ U (To S Це и 2Т ) © Ud k Е ) 


ЈК 


+V (Tes -TA )+W, +X (TH -Те)+ У Есе 
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12. Perimeter 1 


п+1 n+l 
-2 FoR д ш (1 + 2FoR 2g - Ши 


Е U, (и. В = Ти UE E T ) 


1+1,}, К 
EV (Th = ) а (Т. IM ) +W,+W, 
+X,(T vat -Ti,) +X (TH 2S) +Y Бели 


13. Bottom Face 


(1+ 2FoR,)T?,-2FoR,T?,,,=T** 


КА: 1,3. к+1 i,j,k 
ж ж ж ж = ж ж 
+U ЕЕЕ T П 2T k ) 


Е к Í Tk Е 27 yk ) У (Теа Тъй ) 


4 4 
+. +X,(T зи =) +Y Ес ју 


14. Interior 


Soke 7 +(1+2FoR,)0,.,—-FoR,1,,.,= 0/7" 


i,j. k ij k+l i,j,k 
ж ж ж ж 2 ж ж 
+ ОТ Шы + T ma 2174, ) 


+Fo(T ži jkt Кш 2Т У е ) 


4 4 
+W,+X,(T** -T.,)+Y,Ec,,, 
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15. Тор Ғасе 


-2FoR T^^ "(Ре 2 овој ти ERA 


ке i.j.k ijk 
ж ж жж и жж > и ж 
= U, (T 1.j-1.k T je k 214 ) + Fo(T * -l j.k Е oT ) 


+V,(T-,-T**, )+W, +X, (T: Taa EO 


ij k 
16. Perimeter 4 


(L+ 2FoR,)T u 2FoR T a TA 


i,j. k+1 i,j,k 
ж ж - ж ж - * 
+U Ж ы +T N 214 укны (Тн de) 


+V (T=; = Г ) РМ Да тен DEL VE 


ijk 


+Х (Т. Б TA.) +У Ес, к 


17. Front Face 


-воВ ТР, + +2FoR,)©,,— FoR ioe — tee 


ij. k-1 1,}. К jj. kl 
+U neuen ij+1,k a TT ELE M 215% ) 
* U,(T ue c T. ) + V, (Ше. ze ) 


4 
+W. +X (Кон: -TORY Eca 
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18. Perimeter 2 


-2FoR JD a +(1 +2FoR ) TI? = Tee 


t p Eo АЕ sU (I. те шы) 


ij-1.k ij*l.k i,j,k i-1,j,k 


И, эм Шаа (Ше = ку +W., 


i,j. k 


EST 6^ dac T^. -Т-,) * Y Eo,,, 


ij. к ЈК 


19. Согпег 5 


(1+ 2FoR ,)T7/, - 2FoR T7, 2 T?* 


i,j,k 1,3. К+] i,j. k 


EU (Lii E Ті ) + қыт, ТІНІ ) 


11-1. 


У (Т-,- Тил J)! Ve. s DEVE ca а ) 


i.j. k 


4 4 
+, + + И +Х (те -Т-,) 


rx (r> =) вх (т -T2,) +Y,Ec,,, 


i,j. k ij. k 


20. Perimeter 6 


– БОВ „Тила + (1 + 2К0В8 ТА, – КоВ „Теа =T** 


i.j.k-1 i,j. кт 


О а = n. DUE = ds ) 


АРМ (T=, ER ) А (T = sd ) +W,+W, 


D -т.,)+х ти -T-,) +Y,Ec,,, 


i,j, k 
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21. Corner 7 


то Бобища * (1 - 2FoR j)TT, =т (ЧИНИ ВИ и 


i,j. k-1 i,j,k 


+U,(T! B T$ ) + V, (T=, -Tik ) + У, (T=, Тик ) 


Најк 


У (Тео Ти )eWQ We WeX((TAE -Т-,) 


1. }. К 


+Х (Там –Те,)+Х (Ту Ста) ЧУ Бо 


i,j,k 1,}, К 
22. Perimeter 9 


(1+ 2FoR,)T, , — ОБН ween ле 


i.j.k i.j.k+1 


+I (I = E T oer xo ENDS Eo. ) 


јак 


+V, (Te, -TA )+V,(T-,-T** )+W,+W, 


(Gne -Тт.,)-х. (ты Та) ЧУ Во 


ЕТЕ i.j.k 


23. Right Face 


- FoR,T;,,-(1-*2FoR,)T7, FoR, T a TA 


i,j.k-1 ij.k ЈК i,j.k 
ж ж = жж ж ж ^ ж 
Е 0, (Т i.j-1.k T ) * Fo (T i-1,j.k p m jk 21 ijk 


+V,(T-,-T**, )+W,+ Pelee -T,)+Y,Ec,,, 


Lj.k i,j,k 
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) 


24. Perimeter 10 

n+l n+l 
-2ЕоК,Т уа + (1 РОК ЈТ ји Тук + О (Тули - Tik ) 
нот," -1,j.k v E 1 1 DEVE (T=, = ете ) 


+V, (Te, -TIA tW, +W, +X, (THY -Th ) 


ijk 


X Иран - Ti.) +У Еби 


ijk 
25. Corner 6 
ШЕ Онова) оров т um TA 
Щи a (T ij- ne T Tik ) 4, (Ту uc Тін 


жж 


+V, (Te, zi ) +V, (T-, = ааа ) 
+W,+W,+W, +X (тү Xe (Uns - Ti) 


4 4 
+X (Tt -TL)*YEso,, 


26. Perimeter 5 
-Еов Тоа + (1 + 2208 „Ти – БОВ Та Е ТЕН 
+U „(Ги ШЕ 7 Heer П та по ) 
+V, (T=, -Tii МЕЦ EE Jade ue 


1,3. К 


+X (АА К т. - Ti.) +Ү,Ео, 1 
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27. Corner 8 


Бов Тоа t(l + 2FoR yn = Fa EE ) 
+U,(T nik E TE ) um (Т.., plos ) A (T=, zm ) 
+V,(T.,-T**, )+W, +W, + W,+X,(T*"% -TŠ,) 


soc (ous зара aren ае -Ti,) * Y;Eo,,, 


ijk 
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APPENDIX C 
PROGRAM EXPLICIT 


Meck Soko ck desk al al kok hob aR x a de c a Ж Жоп ЭК ж А xc te xc t ЖК гу Р se уга ҖЕ ЛЕ ЖК Де л аи K 
Ж 

* PROGRAM EXPLICIT 

ж 
ИИА # ЖО Ж Ж Ж Ж Ж Ж: Ж: Жк ЖК Ж Е ЛЕ Ж ЖЕ ЖК К Ж: ЖЕ ЖК л: Ж Л: ЖО: Жы 
ЕЕ. 


* 


Geer ICI] METHOD FOR TRANSIENT HEAT CONDUCTION IN THREE 
ВНЕ NSTONAL CARTESIAN COORDINATES 


а 

Ak dk XR RR МИКА em RR KE oR oh ce ch oh oe yc S IE YE МОИ Алам лл МЛ £ f Xe X tk Sk Ez SE e Y Жол Ж 
Me Meh oh wR sh sh sk hk sk tek oe ee oe Se ТИ МА ти ар th eR dox ТИ Moke eM ch oh ete oe se ch eo А оао se a ae Sh Sk йсй И гд сои су сй 
ж 

* DESCRIPTION OF VARIABLES USED IN PROGRAM EXPLICIT 

ж 

ae OTE: SI UNITS ARE USED IN THIS PROGRAM UNLESS 

3 SNHERWISE NOTED 


BEEN SER SPECIFIED VARIABLES 

BM GEOMETRIC DIMENSIONS: 

° ЖОП DISTANCE INCREMENT IX THE X DIRECTION 
^ DELS DISTANCE INCREMENT IN THE Y DIRECTION 
» МЕН ОА СЕ [INCREMENT IN THE Z DIRECTION 
5 [ОО ОО OPD ISATE OF NODE 


E J----Y COORDINATE OF NODE 
c K----Z COORDINATE OF NODE 
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LX---LENGTH OF ОВЕ ТА FHE ХЕ ИКЕ 
LY---LENGTH OF MODEL IN THE и Ле relic. 
LZ---LENGTH OF MODEL IN ПРИ ИЕ БЕО 


M---- NUMBER OF NODES IN THE X DIRECTION 
N---- NUMBER OF NODES IN THE Y DIRECTION 
P----\U MBER OF NODES IN THE Z DIRECTION 


Bo VIA TERIAL SP ROVER Ties: 


CP---SPECIFIC HEAT 
K1---THERMAL CONDUCTIVITY 
RHO--DENSITY 


CONTE 0 le: 


FLEA=SURFACETEL Y 

ЕКЕО--\ С МВЕК ОЕ ТЕМЕ ТЕРЬ ВЕТ ЕЕ БЕ | е 
TENWDEISSSTUETSESPERIDSNOPSIES 

Е-е ЕИ АХЫ СОР ИЕ 1 

TAMB--AMBIENT TEMPERATURE 

TINIT-INITIAL TEMPERATURE OF Aili ODES 

VAL---LOGIC NUMBER TO DETERMINE TTE РЗ ОУ 
ЕЕ PE D 


2. NON-USER VARIABLES 


A. GENERAL: 


BIO----BIOT NUMBER 
COUNT--COLNTER FOR DETERMINING THE NU MBEKOGE Ele 
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WME STEPS 
DELI---TIME INCREMENT 
FO-----FOURIER NUMBER 
TIME---DURATION OF CALCULATIONS 
TNEW---MAIRIX OF TEMPERATURES AFTER TIME STEP DELT 


B. EXTERNAL PROGRAMS: 


CTIME-----USED IN EXTERNAL IMSL ROUTINE FOR DETERMINING 
ГОЛОЕ УО ОВИМЕ Е 

ХОСТ------ Е ЕВА Е lol ROL fl NE FOR DETERMINING 
AwWOUN OL CPL TIME USED 

ШИ ЕО=- CPU TIME AT BEGINNING OF CALCULATIONS 

TIMEI----- ЕЕ РОВ ЕСТ РОО 

СМАСН----- Cob Dawe teh Ne | vSt ROWTINE FOR DETERMINING 
pRrOU NT OF CPL TIME LSED 


Е ЛЕУЗОК 1252525), ТХЕМ (25.25.25) 

КЕСЕГЕ НЫ. RHO. ALFA, FO. FLUX, R. RFE R2, KI 
Peele Ody. DELZ, DELT, TIME, TINIT, TNEW, BIO 
EAC TDIFEI TDIFE2; 1DIFE3 

EESC LA, LY, LZ 

REAL CTIME, TIMEO, ТІМЕІ 

D ECGERTI E M,A, P, COUNT, FREQ, VAL, NOUT 


DA TERNAL CIME, UMACH 


G ОО TIME., DELT, M, A, P, TDIFEIL TDIFE?2, TDIFE3, TINIT 
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COMMON COCA TIENEN 


She she езу hh sie, age NS M A SA A SA SA NA SA IR SA DAA SASA A A A ATA SA а зе Зее Зе ASA D Aa И A AAA 


Фр ale ate ale als «б, е + +, 42 42, чо sv), < OEV EOE OE OE E EVEN E ENV ао аб обо а ENVE TEE E ENE E E ETE E EE E E EVE obs ale 0's < sts У: ats ate ale 0% xÚ < x < we we uo elo uo «o, «йо +?, ob 
% 4% 57 745 795 797 745 547 745 9,7 927 79% 74% 247 зе 797 24% 74% am mae тоз mam oma туз a, tam mu o] Iam tu vam 247 түз 947 547 287 74" 2,5 545 5,5 245 25 5;% 5,5 595 5% 44% 245 SaS San 24% 2;% 4% “3% 24% за 29% 542% 99% 94% 79% 5% 4% 2% 2,5 24 9% 54% 29% 29% 


` ЧОИ G C I SDD IE 


TIMEO = CTIME() 


NOx Xx xoxo x deo ae oh ake сй ту изи зу sa ХА ҒАС АЛСА АҚСА ИЛ ХИ ХАЛ ok ho ok oe oe oe ok ж 


CR I NES a ee ic 1, 49, 42 „ «Po «P 2 EC) 9, è, 2, 6, ate +), 8, 4, 4, О te we to 4. 42 -+ 
sk z oh ofp oh ok s зи зи Зи Зи уа зе Ум Ум oe е Ма he oe She oe Se ok Мо Ме Ма мо мо sk Ме За ам зи е Не Rz OO Oe HK OK HE Oe Se MK sk Жен Уз fc y fc s< xx yas Ya das ve Ye sese sess 


VARIOUS PARAMETERS TO BE USED THROUGEHOUT НЕ РКО ЕН 
WILL NOW BE READ TN (Os Hie RO Greer 


УЕТТЕСУЕХТЕКТАМВОТИЕАУ ИЕ ФЕ РЕ ОРГА 


МЕШЕС ЗУ ЕХТЕКТТХІ Т ТИЕ ГХИЯЕЛЕ Р ИРЕ 
КЕЗ I] II 


READ(*.*) I 

WRITE(*.*) ‘ENTER K1, THE THERMAL CONDUCTIVITY.’ 
READ(*,*) KI 

WRITE(*.*) 'ENTER CP, THE SPECIFIC HEAT OF THE MATERIAL.’ 
READ(*,*) CP 

WRITE(*,*) (ENTER RHO, THE MATERIAL DENSITY.’ 

READ(*,*) RHO 

WRITE(*,*) ‘ENTER FLUX.’ 


WRITE@,*) ENTER. LX, THE PENG? ae 2 ЕЕ ОЕ 
КЕ ВИ ОВ 


МЕПЕС,) ЕХТЕВ ЕУ, ЪНЕ ЕХО wie Se ele. 


WRITE(™,”) EXTER LZ. НЕ БЕО ieee Giles 


8§ 


ENDO LZ 

ПИТ ENE ММ ТИЕ NUMBER OF SODES IN THE “DIRECTION” 
ЕАО <, =} М 

ИИК О ЕЕК К THE NUMBER OF NODES IN THE Y DIRECTION.’ 
EIEAD(*,*) № 

Sees. ) ENTER P, THE NUMBER OF NODES IN THE Z DIRECTION” 
ШЕАГХЗЗУР 

E S ESTER FREQ, THE NUMBER OF TIME STEPS BETWEEN 
ШЕРГЕ”) SUCCESSIVE PRINTINGS.” 

READ(*,*) FREQ 


Pee) CALLED USE 1: IF THE VALIDATION SUBROUTINE IS NOT’ 
КИ TEC“) DESIRED, USE 0. 
READ(*.*) VAL 


. .....,..,.....л. Ба ............,.7 3°. КА x ate e'e ge sta she sha ste ote 2% 590 ale ale ake ate ale ats als ole ahs oh ... 
5 255 55 55 2,5 2,5 2,4 245 Tam am am cat at за“ Gye See 245 2,% 24% 2% 4)% 7,% 7,% 2% 74% 7% 2% 5% 2ҙ% зе 24% за Резо за“ == Myr NR 


meee ly ONS OF CONSTANTS TO BE LSED THROLGHOULT THE 
PROGRAM 


COLNT 2 0 
ШЕ = 0.0 


DEIN = LX (M-1) 
DELY-LY (N-1) 
DELZ - LZ (P-1) 


R= DELX (DELY*DELZ) 
R1=(DELX**2) (DELY**2) 


R2- (DELX**2) (DELZ**2) 


ALFA- KI (RHO*CP) 
BIO-1I KI*DELX 
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FO-0.5 (1- RI R24 BIO*SQRT(R2) 
DELT = FO*(DELX**2) ALFA 


ә. 
4 


NOTE: THE FOLLOWING EXPRESSION FOR DELT WAS CSEDT OR 

т VALIDATING THE EXPLICIT TECHNIQUE WITH THE CEOSEDOBEOEST 
SOUTIEN: 

ж 

^- “DEPT=0025 ИРЕК а 

ж 

ЖИ ЕЖИК ЖАК ya per ye yc sk sk Ж 


s s zz šE e tente test tete tee mone tote t zh ДИ СЕ ДЕ И ЙЕ г ДЕ ДЕ ТЕ tede td dert dete tet td et tette td dett d tata fk ху ЕЕ г ЕЕ № 


INITIALIZE AEL TEMPERATURES ТЕП ИЕ ИЕЛЕ 
TEMPEK STURE: 


DO 101=1. M 
DO 15J21. * 
ООО Цена 


Jj b S] | 


20 ОЕ 
Јн АЕ РОУ DT 
Ioco SISSE 


МИЗ ВИ ОС р ЕИ $ he hh ok he oh oe hs ok ft He ois eo Ee HA dd gk d Rd A dk ОК Me ok ok ЖЖЖЖ ok ok Ns ok OE k 
ko XN XR oe He И ы Лл си Ж Зи ги 2и г э ОК oe oh ke ok ek ck He Re Hem Hh Hak RK жож ж Ж Xk k sk MRK IH RK KK HK А 

22 

* THE OLD TEMPERATURES WILE NOW BE SAVES FOR ПЕЕ СЕИ И 
= EN IHE TEMPERATURE EQUATIONS 


DO 30121, M 
DO 55 J21, N 


90 





РО ЕР 


T(J.K)5 TNEW(LJ.K) 


40 CON TENCE 
See CONTINUE 
CONTINUE 

ж 


oe Me ТАЗИ ри За зи tok Хр зи дима ти Зи зи he oe ae ae oe ee ee ae ee ee ae Mee oe a ЖЖ РМ sk Ne ahh oh a Ne ok oe Ne ee ae eh oe ok oe ok ok se ok Зи ЗА ake oh op ok oh SE 


oh we oh Нему ме ЖУ Не Ус cfs ss ote ols af sks ak stk o's ofc oh oh ots oe oe oh oh fx ok oh x deb x ale ЗМ Зе за е И: Уж ЭУ У Мо X p Ye y YZ Ж: yz yk 3 2 yE rh oe ok oe 


~ 


fee SCRE WENTING TIME AND COUNTER 


I C OU NT = COUNT + I 
MNE = TIME + DELT 


ж 


Yo af 2” 
NO. 


уе за уж зи Фл Фе Ур nfs ah ste sta ote ole ate ote ate ate afa afa 00, fo s*a a'a ala o'z op ato o'a afm e'e o'e sha аба фи Ма ум «е Ме зе Ме Ум е ale pha s'a ale ate st oh ale ole o'n Ул ул Ул в, ote x sie oh ehe > t» ух ~ ее, 4“... 
+“ =e" 18% 70% 54% зае %,% Mee Bee г." "а? за“ .. .4% 98% 64% 24% n 99 оа 5,“ 5,7 Bee од =," .9" .4. n, ." =e" „а. .. aa? ee" ~" а" “` це 2. ~." ө“ sen =e оре 4.. Sgt See See Fee + во оо" чо ое оф чо Pen У," . е? 74. ҒЫСЫ 2. 7. “~. 


06 ato аб обр ado ohe so Doe do o^g oe Jo wo alo de sto o*r ado nta Po Jo Pa Чо ato afo sta ata ВУ ate Ja Jo Fe ato da to da Jo do do to Je nto Jo Po Po aho ato ato 
е еде офи оде өүө Фу» Фф? өрө өш» өү? чог ефе едг еде еде оде еде еде ефе еде оце дуо еде вде еце рое еде ода гоч еде ғор еде 240 91 040 042 040 9,6 94% 94% tj tq tgo ғор еце оде еде еде еце 


5% 
+ 
` 
+ 
` 
` 
` 
"4, 
.. 
> 
52 
`. 
52 
5 
26 
26 
4 
& 
se 


feo Or EOLATIONS FOR THE THREE DIMENSIONAL BLOCK ARE 
BENED "U^ PRESENTED. 


ale ato oto ја ч x обр ohe oo do s'o a'a a% ata ee 402 «02 «42 272 абр абр абр ae «4, als ab ato ale oh ah abe ale af або als ale we ў x Po we ate ole 2р 3% аб, бл б обр бр абл абр Se ats We У“ уг Jo y: > «fo чб о 0% 2% aha a% „бр 49, обр af 
Per Mer 8,6 Ре" dpa Feu Fae врп 5,1 7,9 2,0 2,9 5,9 549 4,9 2% 44% 7% Ppr Pga Sgr PoP 504 20% 2% 24% 2% 49% 7/4 20% 2,4 29% 2; 29% 2,9 5,5 2. 2; 2/4 4,5 2,6 5,» 2% 5,5 2/9 2,0 2,2 2,5 2,5 4,4 2)» 2,7 290 24% РА 2,5 5,6 2,5 2,5 2,5 ғ,% 7,8 Pa? cao tao cao 242 одр рар рдо ТР 


Jo Jo Jo ato do ato ato Jo nto Jo 4» -» 4)» Jo Jo Po Jo Po oo Jo ate ote Po Jo <” Lo Pa oo Jo Lo to Jo «до Po oan 4)» ale 4)» Jo Jo Jo Jo Jo 4% мо Jo Po ata Po sto Jo ate we we Jo we ae Ја Jo Jo Jo o Pa Ja oo o". Jo aha do Jo Va Jo 
497 wo" Sor Po? oP PoP до“ рар роп де РР оре РОР PAP Pa" POP 6,0 949 06,0 9,6 64% 9,09 9,6 680 гор еде чоч едр рае pP Phn Pan Pg? едр еде едр гоч едр едр 6,- 9,6 94% 2,6 9,6 640 оф? 68% 9,4% едр еце ео? рог едр еде рое еде еар 0,9 940 «40 одо пр одо ој? 2,6 6,0 946 ее ера ед едх едр 


A TERIOR NODES 


РО 50 1=2, М-1 
DO 55 J2 2, N-1 
DO 60 К 2, Р-1 


ЕЕ) РОО + T(IJ+ ЕК - 2“ 
ШЕ” [OS T(IJK+I)+ HJ KI) - 2° 1011.K)) 


91 


& + FO о р 


60 сом Кр 
Sy =a CON FINE 
50 -CONTINGE 


ah sh oh ЕАР ЕЕЗ СС ЧЕК Ял: Me 2 9% 9% ok 9% 9% aie fe co xxx ek ok of ok ok ck ok ok ok ok ok of ok oh ok sh У ok ok oe SK of ok За за За У 2% ok ck oe 


+ 5 
HR she xk x ve sh she Se sh sh oh ate sb Ske oh os Se a che ak oe Мо е ake oh ce Se НЕ we see Ма СКМ СЕКСЕК КЕТ КК sx y 


CORNERED NOBRES 


7 СОКЪЕКО БИ 


ZA (24 === 
l | dg 
— — 
— 


TNEW(LJ.K}= T(LJ.K) — 2*R2*FO*( T(L.K - D) - T(LJ,K) ) + 
& 2*RI*FO*( T(LJ-- LK) - T(L,K) ) + 2*F0*( T(1+1,J.K) - 
& T(JLK))- 2?FO?RI*DELY*FLUX KI 


С СОКЕ 


C. — 
|| || 
— ‘77 

— 


TNEW(LJ,K)2 T(LJ,K) + 2*R2*FO*( TK 1) - ТОВ Е 
& 2ЗБКІЗЕОЗИТІЕІЕКІ- ПЛ TOTE 
& T(LJ,K)) * 2*FO*RI*DELY*FLUX/KI 


E ORNER NR. 3 


TNEW(1,J,K) = T(1,J,K) + 2*FO*( T(1+1,J,K) - T(1J,K)) + 
& 2*R1*FO*( T(1,J+1,K) - T(1J,K)) + 2*R2*FO*( T(1J,K-1) - 
& T(LJK) ) * 2*FO*RI^DELY*FLUX K1 


BOR NER NR. 4 


nod 
ng 
pep 


een WK) = IG J.R) + 2°FO*( I(I-1H,K) - T(IJ.K) ) + 


ПР СО “( ((ј. ЈК) - TUSK) + 22R2*FOSC (IJ, K-10) - 
ГК) 2*FO*RITDELY*FLUX МІ 


EMUORNER NR. 5 


= 


= 
~ 
K=] 


= 
| 
TNEW(1N,1) = TUN.1) + 2*FO*( TQ2 „Хл) -Та,Хи)) + 
& 2*R1*FO*( T(i,N-1.1) - T(1,N,1) ) + 2#R2*FO*( TUN, 2) - 
& T(1,N,1)) + 2*FO*BIO*(R1**0.5)*( TAMB - T(LN,1) ) 





| 9; 


“CORNER SRI 


E S а 
yi 
А 


TNEW(LJ,K) = T(1J.K) + 2*FO*( T(1-1,J,K) - ТАЈК) ) + 
& 2*RI*FO*( T(1J-1,K) - TUJ.K)) + 2*R2*FO*( T(J,K+1) - 
& T(1J,K)) + 2*FO*BIO*(R1**0.5)*( TAMB - T(I,J,K) ) 


* CORNER NR. 7 


I 
J=N 
ВЕ 


TNEW(UJ.K) = T(IJ.K) + 2#FO*%( T(I + LJ,K) - T(IJ.K) ) + 
& 2°RI*FO°S(T(IJ-IKO S НЫ O o OS np 
& ТОЛКУ) ОРО ВТО O ae ІШПЕЙ 


~ 
t 


СОК БЕР Каз 


i 
JEN 
K=P 


ТМЕК) = ТОЛК) 2 ВОТ ае = 


& 2*RI*FO*( T(J-LK) - T(LI,K) ) * 2*R2*FO*( T(LJ,K-1) - 
& T(LLK)) * 2*FO*BIO*(RI**0.5)*( TAMB - T(I,J,K) ) 


Shoe Se oe ok oie of oe of me ok se ЖК oe sh ok oe o doo X xo Me ke ok xxx xoc x oe Ж: Жж Ж. ЕЕ 


9% 4“ 72242 ж sie ste sh 2% Ма > sts ste у» 2% y 22.5... у ,.. 
9% 54% 25595 2,5 үз age ru rm AS uS NS. NS AS ASA... GR Iu tye mye syd ду 
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ADIABATIC SURFACES 


ж 


M EOST YZ SURFACE 
DO 70 J=2, N-I 
BOSK = 2, P-] 
ЕМ 


ПА ЕМУПЈК) = ТаЈК) *- 2*FO*( T(I-LK) - T(I4JK) ) * 
& RI*FO*( T(IJ-I,K) -- T(IJ 4 LK) - 27T(LJ,K) ) + 
& R2*FO*( TU.J.K-1) + TUJS.K 41) - 2°T(1J.K) ) 


KF CONTINUE 
ШОСПХІРЕ 


we 
9° 


* BACK YZ ADIABATIC SURFACE 


DO 80 J= 2, X-1 
DO $5 k 2 2. P-1 
[= | 


ВЕЛЕ KR) — T(1J.K) t 2*FO*( T(I- LJ.K) - T(,J.R) ) + 
EE ISEOCCI(LJELK) * T(EJ- LR) - 2*1(DLJ,K) ) + 
& R2*FO*( T(LJK-1) + T(1J,K +1) - 2*T(1J,K) ) 


fo) CONTINUE 
Bue CONTINUE 


ж 
sk z de he he ts Me Me Т de 46 Xe Ae de Se e ds АЛСЫН sie ak ap de sk s аи ПМ ЖК КЗК Ke С СКО kok dek k oi k ОИК Ж Фа Па ДА Я: 


> ate ale sh ake ale x: ate ate atp nis sts ба str to ato ahs Yo <°, x, ole ots Ми 279 3% Мо зе Мо фи а slo aly ale ale әб a's ~, ж aly ste ale ale So ole ate ats obs we ole эбе зб» э» ab * ва зо зб о0о збе збе або у и чо ове ма мо Ла 492 абе 48 
5 54% 49 205 зе“ 20 дуз Sy ора ро 209 949 949 909 949 9098 909 49% «,% «9% Суз у“ зе“ 5,9 =“ зе“ зе" 25551 See Spt Tye Hye Seu See See Age MG 5,5 55 51 25 2,5 5,5 545 595 295 59% 24% 79% 29% 29% 2% $.9 24% 29% 20% 2,5 29% 2,5 2,1 20% cat Sa 29% 2% 29% 2,5 29% 24: 46% 





* BOTTOM XY ADIABATIC SURFACE 


DO 90 T= 2. M-I 
DO 95 J=2, N-I 
КЕ 


TNEW(LLK) 2 T(LLK) * FO*( T LK) * T(I-1J,K) - 
& 2*T(LLK)) 4 RI*FO*( TUJ-LK) + TUJ+ 1K) - 2 ДК) ) + 
& 2*R2*FO*( T(1J.K+1) - TUJ,K)) 


Оз CONTENT 
ЭО ЕРЕ 


* TOP XY ADIABATIC SURFACE 


DO 100 I2 2. M-1 
DO 105 J2 2. N-1 
Nae 


TNEW(IJ K) = T(LJ.K) + FO*( TWd+1J5,K) + TU-1,J,K) - 
& 2°T(1J.FE)) + RISFOS( (IJ LK) t T) JP ИК p 
ОЕ вю 


io. ОБ ШЫР 
100 CONTINUE 
ж 
s: 
& з o f tt f d te te fe te ЕЭК ss teo dt te c dk sk ad d К ЖК ЖКК ЖЖЖ ЖЛ Ж Ж Ж ЖЖ у 
sk s te te the fe de siz sls sls te de oh de ecd de efe d o ob dete ode of te def oe o de te odo o oe dioe ooo d o dote dede ek ЖОЕ oot dece sk sk sk fo КОКО ЖОК 


ж 


"o  RIGH L.XNZ-GONVECTINE St REAGE 
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ЮО TIO I7 2. M-] 
DO 115 К =2, Р-1 
А 


ШЕ -ЕУСІЛК) - T,J,K) + FO*( Ti1+1,J,K) + T(I-1,J,K) - 
& 2*T(LJ,K) ) * 2^RI*FO*( T(LJ-L,K) - T(LJ,K) ) * R2*FO* 
5 (ТОЛЛК-І) - T(OJK-10) - 2?T(OLJ,K) ) + 
& 2*FO*BIO*(R1**0.5)*( TAMB - T(I,J,K) ) 

ШІ CONTINUE 

ШІ СОХТІХГЕ 


ы, 46 25 ats whe абл оби ate „4 ate ahs ats ah ole ate ale ah cfs 46 26 a's Ма +. at ate af ahs ale whe обл абр whe ale ats «Ко чо oh ale ak als whe wh зи 7: n REPE E TE TEE TE E PERSE E E E MN PE M M ME PE 0% 
wit oe oh rs te fe she фа А 20 35 130 25 2% 14 5%. 542% 2% р М си м X ner АА ЕЕ ЧЕ АЕК ЗМ за Ме Y. ЧУЕ ЕЕ ЕЕ ЧЕ ЗА СЕЧЕ УО Ма зе (9 e rte te Ps y o n tx: 


ПОКЕТ АКЛ КОСК +f ale LJ. ate ale 2% ate mts ots ата w. +ç <. <Ó ate ate ate 494 ate ate ate ate ale во Г% 2% Г% 4% 4% -% СВ ate fo ate ate ale als whe ats ats fs або э? ate ats x. ale ale xf. про ale ate we 47% ... ... аб 
24% 54% 949 989 949 949 949 9490949 948 949 949 SeS Sat бо 79 949 505 94% 29% 909 949 46% 94% 949 949 909 949 949 949 949 зе“ чу tg t) То“ 545 Syd 595 505 505 44% 94% 94% Tye HOS 54% tam taS Sot 249 94% 96% 94% 54% 949 949 944 tgo tuo tah tat 949 tat %4% 4% 949 04% 44% tam tat 24% 


SUUSEFT NZ FLUX SURFACE 


156-120 [2 2. М-1 
WO 125 K=2, P-l 
Ј= ] 
ШЕ ОЮ) = ЦЕЛК) + РОСТОК) + ТО 1J.K) - 
EENLUIUEPROS 4 2*RI*FO*CT(IJ- LR) - К) } + 
EN TUO CIAR: D — T(EJ.K 9 D) - 2“Т4Ј„К)) + 
ENS LL O-RI*DELY*FLUX KI 


ШӘ CONTINUE 
ПРИТОМ ПА СЕ 


ж 
= 
к 


ЕО e X ai x Y eae Yr ee o o o e елле Жж А ЈЕ КСЈ ЈОКО O K x x 


Be «Je a^ ats le We als ats ale ale als no no se n^ no o So M. uo a. o co n. n. ne s sk sh als eo We >. ~. Wh Де 7 we ~. ~. ме ме ме е ме Дем «do u^ co ato co mo. n. co uo co co чб. чў же 
545 24% 54% 54% 24% 44% 545 54% 54% 99% 90% 505 %0% 29% 99% 54% 54% 54% 24% 44% 24% 24% 24% за“ зе“ 905 HHS SFO TES Sad 248 245 24% 549 59% 98% 54% 447 24% 54% 24% 54% 24% 74% 94% 94% 50% 24% 24% 24% 54% 94% 54% 94% %9% 2)% 94% за 94% 94% 94% 49% 49% 945 24% +É + 


~ 
го 
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OUTER PERIME LE Reh ЕЕ 


РЕВ ТЕТЕ ВЕ 
DO 150 1-2, Х-1 
к= 
К=Р 


ТХЕУСІІК)- ТОЛТКУу- БІЗҒОҢ ТІЛІК + Тат. 
& 2*T(1J,K)) + 2*FO*( T(1+1,J,K) - T(LJOE) ) + 
& FRITO Te Po N Ó 


PCOS TECE 


> 
A 


ЭКЕЛЕ - 


DO 135 J= 2. `N-1 
[= М 
K=P 


ТЕЖИК) = ТИК) БОО је 
& 3*T(LEK) ) * 23FO*CTCI TERRE S TOI S 
K RTO TIIKI S PRK 


рабо ее 


x 


“PERIMETER AR? 
DO 140 J=2, N-I 
i 
k=! 


95 


БЕСТІК - TIJEK t RIFO T(OJ-LK) + T(I,J + 1,K) - 
pee ТО Ку) ғ 2ЕОЗСТ(І-ІЛК)-ТОІЛЕу)- 
& 2*R2*FO*( T(L.E.K * 1) - T(LJ.K) ) 


140 CONTINUE 


~ 


ж 


ШФЕКТМЕТЕК ХӘ. 4 


DO 145 J=2, `N-1 
В] 
Ке | 


ESTER = ТОК) + RI*FO'( I(IJ-L.K) + T(I.J + 1.K) - 
EE Ico FOS TONO - ТОК) + 
EN E FO (CT(J.K 10)- T(IJ.RO) 


145 CONTINUE 


E ERIMETERNR. 5 


DO 150 K=2, P-l 
I 
nS 


ШЕПТІ) - IJK) t REON TOJ, KFI) + ТаЈК-1) - 

& 2*1(1,J,K)) + 2*R1*FO* ТАЈА К) - ТаАЈК) ) + 

& 2*FO*( T(I-1J,K) - TUJ,K)) + 2*FO*BIO*(R1**0.5)*( TAMB - 
ИК) 


БОСОХТІХСЕ 


ж 


29 


* PERIMETER Ы 


DO У и = ul 


[=] 
J=N 
TNEW(LJ,K) = T(LJ,K) + R2*FO*( TUJ,K+1) + T(1J,K-1) - 
& 2*T(LLK) ) * 2*RI*FO*( T(LJ-,K) - T(LJ,EK) ) 4 
& 2*FO*( T(I- LK) - T(l,K) ) 4 2*FO*BIO*(R1**0.5)*( TAMB 
Gag SUR NK Y ) 


ІЗ SE Е 


СЕВЕ Е 


DO 160 K = 2. P-1 
| sj 
ps 


TNEW(IJ.K) 2 T(LJ.K) ^ FO*R2*( T(LIK c 1) SDT YI IE 
& 2*T(LJ.K) ) * 2*FO*( T(I-LJ,K) - TK) ) + 
& 2*FO?RI*( T(LJ- LK) - T(LILEK) ) + 27FO*RI*DELY*FLUX.K1 


СОСО УМЕЕ 


* PERIMETER NR. 8 


ж 
DO 165 K= 2, P-1 
[=] 
J=] 


ТХЕ И Б) = ТОЛК) Е КӨЗЕ (TJ O I IT I O E 


100 


ВЕ ВОЧ ТОТЮ ІК) ) + 
O RITTER) ALJIR) F O PO” RIF DELY* FLUX ‘KI 


tes CONTINUE 


ж 


x 


ЕВ: МЕТЕК ХА, 9 
DO 170 122, M-1 
J= N 
K=I 


| НЕК) = T(LALK) - 2*FO*RI*( T(1J-1.K) - TU1J,K)) + 
EMECEO*RO'( I(I-EK 1) - (LK) )- * FO*( T(I+ 1,J.K) + 

КИ КО) - 2*T(I VR) ) * 2*FO*BIO*(RI**0.3)*( TAMB - 
EXEUEULJ.K) ) 


КОК ОМГ МЕ 


x 


M ERIMETER NR. 10 


DONS I= 2, M-I 
је 
Кер 


TNEW(J,K) = T(LI,K) * 2*FO*RI*( T(LJ-LK) - T(LLK) ) * 
& 2*FO*R2*( T(JK-D - TLLK) ) * FOSC TOI LK) + 

& T(-LJ,K) - 2*T(1J,K) ) * 2*FO*BIO*RI1**0.5)*( TAMB - 
& T(I,J,K) ) 


We CONTINUE 


~ 
Du 


* 
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” PERPE TERT KS 


РО 150 [=2, М-1 
J I 
K=P 


TNEW(LJ K) = LIR) F FO TOIA ИК ОЕ 
& 2*T(LLK)) + ЗЗВОЗВИ ЕВИ И И вок: 
& 2*FO*R2*( T(LJ,K-1) - T(LJ,K) ) - 22FO*RI*DELY*FLUX/KI 


СЕ 


“РЕК МЕТЕ puls 


DO 185 I2 2, M-1 
J=] 
K=] 


INEW(LEK) — ТОКУ ВР РО T(1+1.J,K) + T(1-1,J,K) - 
& 2°) ) + FO RA etd Е 
& 2*FO*R2*( T(LJ.K - D - T(LJ,K) ) ^ 2*FO*RI*DELY*FLUX KI 


Iso CON VINEE 


ete Me sts NE che Me sks фе мо Јо te she te ate ote Ма фе пее he до а de ge ake se e ge e 
See Ja Ja" Ta cam e^ fae > She ЗА МУ За ЗА Pam ca Da OS C 1 > ТЕ > МА МАЗ = > 


ote ok oe Me ake ok fe he a oh ООО С kO R Rk 3 oe o fe oe a ok DER ДОМА ДЕ Ме Пу oh oe oo Oo SE 
Ж z eo МЕ Ae oe oR ak ae ae Ree Ae oh eK oe Me ole eM oh ee Me aie ate Д Me a ae ae he Me a Ха ea ake ae a ate кок Ж Ж Ж Me ok Me He I eae He HE EE R R 


ж 


LPDATING THE OLD NODE TEMPERATURES ТО СИНЕ ЕУ 
М NODE TEMPERATCRES 


DO 200 1= 1, M 
рога ] N 


BO 210 K=1. P 
TUJ.K)= TNEW(I,J,K) 


210 CONTINUE 
K EC ONTINLUE 
200 CONTINUE 


ж 
Me ale ate as fe ats ake ake sth ok tes ake ate abe oh abe ale te ake a ole ok sk ok oe ke ok oh ole ake Se te ote se ae she oe ak oe als ake oe oe ole ok oi oe ob ae Ж ok ak sh Ж Ж Ж de ak oh oe ae ak ok de ok she ate oh ok 


4, е, -+ ~, е LO г, ~ - ытта ~ 4, m е, ~, 2 » ~, +, EA -4 2, + СЯ е, + c - 9, as 
2% 26 Ха А OK МА За а МА МА Ха sh ok ok sh oh oh oe 3h 2h oe 2k a me ake ok oe fe ok he ae ie oe ee oR eh OK Ж Җ Ж: Ж Me Sho sk ok oh ae ЗРЕО ОР ЗЕРО f SE. YR Ха АЗ 2 


ПИ РАВКТ ОЕ РКОСКАМ DETERMINES IF SUBROUTINE VALID OR 

ИШ ОШКО |] NE VALID IS TO BE CALLED. IF THE USER DESIRES ТО 

ПИ“ ТБАТЕ ТНЕ РКОСКАМ IN CONJUNCTION WITH PROGRAM BRIAN 
ES OEROGRAMNALID. THEN THE SUBROUTINE OLTPUT WILL BE 
BA PASSED. 


IF(VAL .NE. 1) GO TO 250 
IF(TIME .EQ. DELT) CALL VALID 
IF(COUNT .NE. FREQ) GO TO 300 
CALL VALID 

| GO TO 300 


ВИИ ИМЕ ЕО. DELI) CALL OLIPUT 
eee OL NT NE. FREQ) GO ТО 500 
CALL OUTPLT 


* 


ШЕ ІЕТТІМЕ .LE. 3600) GO TO 1 


* NOTE: THE TIME PERIOD WAS ARBITRARILY TAKEN 
е AS 3600 SECONDS (1 HOUR). 


ШӘЛІ” 
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PRINTS 

CALCULATING AND PRINTING Chl umm IE 
TIME] = CTIME() 

CALL UMACH(2,NOUT) 


УКІТЕСХОГТ,% CPU TIME (SECONDS) = 7 IME TIMES 


STOP 
END 


Ys ate ate 4, ole > ....,.,......... У: ate at «мо ъм ote ale 2: wie ate ale ale EE > ye » at st x“ de de we she te ~ ale a's ум ale als als ale ote 
5 4,% 2,5 %,% 2,% 2% 2,5 2,% Spn sen Se Sue tu Jue Sym 44% %,% одо %4% 74% 4% %)% 49% 44% 49% Hee Дъ Дет mys Mee Dye Jes др ди ти Зи age oye aye ge тұз др те 2,5 ғұ% 2,1 tgs 49% 49% 29% 


Ye ate st ate ote > 2“ ote o's a's а, ate ale oe n's s's обл абл абл лл ale абл абл әб бл обл ote ats ate ob абл абл пін обл бо ум а ste vie ate ale a's ole ale абе ум ole ate 0% ab abe x. abe als ale ale ale ух ate ale ale ale a's ale 
$,* Mee 49% дуо 2,9 SaS re> SeS reS Ser met rgo Soo mu Hye Bee 94% %4% 4,4 24% %,% %,% 2,% 4,4 7,% %,% 49% 44% 44% 24% 94% 44% Ber 347 Bye 79% 44% 505 24% 44% 2% 24% 24% 24% 54% 2% 4,% 29% 4% 2): ғұз ғұ: 29% 49% 408 54“ за Сул зоб че“ 5% 


at ab ot ale a 
` + И ита 


ТН ЗС ВКО ЩЕ МТЕЕ ВЕТ ЗЕО ПОКОРІ С ТББ ODE 
TENPERATECRES: 


СЕКОЕ ОТР 

DIMENSION T(23,25.25) TEW 23 2 To 

REAL IA MB-TII-RHO. ALEAVCEO DEC S ISP 
KEAL DEEN, DEL. DELZ. РЕГ МБА МЕ КЕ Н) 
КЕ аа 

ЕБ А 

REAL CTIME; TIMEO F TIMEI 


INTEGER I,J, Kk, M, X, P, COUNT, FREOWV ARES el 


COMMON TNEW, TIME, DELT, M, N, P, TDIFEI, TDIFE2, TDIFE3, TINIT 
COMMON COUNT, FREQ 


WRITE TIME = TLIE 
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PRINT* 


Ж 


350 DO 370 K=1, P 
WRITE(*,*) TEMPERATURE DISTRIBUTION ON PLANE' 
WRITE(6,360) 'K 7 ',K 

360 FORMAT(A3,12) 
PRINT^ 
PRINT 380, ((TNEW(,J,K), J=1,N), 1=1,M) 

380 FORMAT(IX,F6.1.1 X, F6. 1,1 X,F6.1,1 X;F6.1, I X,F6.1,1 X, F6.1,1 X,F6.1, 
& IX,F6.1,1X.F6.1,IX,F6.], IX,F6.1:7) 


PRINT* 
PRINT* 


EUSECONTINUE 
EE LI NITIALIZENG THE COUNTER TO ZERO 


POUNT = 0 


е 
2 
RETURN 
whe we e Qe ate vto “}» <° Yo «је «Ро 2” чйо $ whe мб» ate ate ate whe at ake Фо «а o qu Mo Mo co do uo Mo so qu qu so no qo «Po Se ste до че че че е Јо че «је So «ја «је o vto Jo Jo «је te «је «до «Во «је че So te Ws So ate 
DaS $qV 744 7% ч;4 2,% 0,% Фуч 46% 24» 2ач 7;- 2,% 4)% 29% 04% 04% 2;% 29% 2,- 24% 2,» 2)" "qe 24" 24% 2;% гол одл оа од 20% 24° 24" 9)% 29% 04% 244 94% 2)% 244 2у% 2ұ% ғұ% 2,4 2,5 245 2,5 597 44% 20" Og" д, % 24% 2% oe 
уа 2% X “ EIS 2 Ж 2 x ГА 
= 
SUBROUTINE VALID 
he 
oq 


WAR emo LON 1 25,25,25); T NESV(25,25.25) 





Е ТЕЛЕ 1 H RHO, ALFA, FO, FLUX, R. R1, R2, K1 
Ба РОА p j INA PELZ. DELT, TIME, TINIT, TNEW, BIO 
ПАТИ РА О РЕЈ, ТРТВЕЗ 


REAL Wie 
REAL СТВ СОИ Е 


INTEGER LJ, NM; NB COUNT ВЕСЕ ЕЕ 


% 
COMMON TNEW, TIME, DELT, M, N, P, TDIFE!, TDIFE2, [DI FEs ИЦ 
COMMON COUNT, FREQ 
ТӨБЕ = TREN ST 
TDIFE2 = TNEW( 5, 5, 5)- TINIT 
TDIFE3 E I EW il 


IF(TIME .NE. DE; T) GO TO 450 
PRINT? 

PRINT“ 

WRITE(*.*) TIME(SEC) TDIFE] TDIFE2 TDIFES 
о | 


450 PRIN T60, ТИМЕ ОТРТРЕЕФИ ФТЕЕО ИТЕР 
ОТРОК МА с 


? REISIILALIZISC PIL C O C А АКО ТРАЈЕ О 


COUNT = 0 


RETURNS 
END 


ae же зу зу зу гк з «Y th Ж Ж Ж fe e e se ok ok ЕЕС СЕЕ ЕСЕ СЕЕ ЗЕВС Ж 4 Ж Ж Ж eo oc ox xo I x ex xo ok ok OK Ж 


Be oe te И ТИ МЕ Me ok oe ЖЕ МЕ ЖЕ МЕ Пр ТИ ТИ ТИ Roo oe Mn ТИ ТИ ТЕ he По aek iot zok aok Мр Ир ИЕ ДЕ Е k Е ЛЕО НОСОК ИЕ А ПА МА МО ДЕЛА НЕ Но Но ноя 
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APPENDIX D 
PROGRAM BRIAN 


X 3 X: XZ S 2 m уй E Z I Z $ I Xe $ Ç 2 Z еже £ Ik б фи: си Жоор TR к б К< Ж Жу 2: Ж Жол ж:ї EK X S Ec S $ 7 NK Qk а б ЖЖЖ f f Z $ 


* 
* PROGRAM BRIAN 
Ж 


киоска коже е 


ж 


* BRIANS METHOD FOR TRANSIENT HEAT CONDUCTION IN THREE 
* DIMENSIONAL CARTESIAN COORDINATES 


яя куки X. fe Se 3 Y Fr X: YK YK Y E ж Ж: ги Ж fe Ya fF 5 e S: к Ж Ж ЖЖ Ж: Ж: Ж ш ККК КОКЕ 
fom de n Зей: 2 г ЕЕ QU ONCE ETE ONU GR NU E C m m ege x e GR E m mou o c RO o р и р эб ти E 
с 

*  DESCRIPIION OF VARIABLES USED IN PROGRAM BRIAN 

* 

fee. of UNITS ARE USED IN THIS PROGRAM UNLESS 

5 OTHERWISE SPECIFIED 


BEEN oe SPECIFIED VARIABLES 

EE GEOMETRIC DIMENSIONS: 

У DELX-DISTANCE INCREMENT IN THE X DIRECTION 
* BELY-DISTASCE INCREMENT IN THE ¥ DIRECTION 
Ë DELZ-DISTANCE INCREMENT IN THE Z DIRECTION 
Ë I----X COORDINATE OF NODE 


° J----Y COORDINATE OF NODE 
. K-Z COORDINATE OFRNODE 
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Х--- СЕМЕТО ЕИ ОЮЕ Е ОРЕ ШЕН 
ГҮ-“ГЕХСТИ OPF NiODP PP IN Y IF T O jO 
LZ---LENGTH OF МОБЕЕ 7 ee eo 


М----\Е МВЕК ОЕ \ОБЕЗ 1 ТЕ ВЕС RS 
N-A L MBER OF NODES IN THE ЕИО 
P---- NUMBER OP NODES IN THE Z DIRECTION 


B. MATERIAL PROPERTIES: 


CP--SPECIFIC HEAT 
КА ТНЕКА AL 'CONDUCHIVI1 YX 
RHO--DENSITY 


C. INITIAL DATA: 


SL BSCIPTS HAVE THE FOLTON IAG MEAD: 
I LEFT FACE 
2-RIGHT FACE 
3-BOTTOM FACE 
4-ТОР ҒАСЕ 
a RON] PACE 
6-BACK FACE 


1 


BIO1,2,3,4,5,6-----BIOT NUMBER 
DELT------------- -- JIMEINCREMPE | 
EPS1,2,3,4,5,6-----EMISSIVITY 
FLUXI,2,5,4,5,6----SURFACE FLUX 


FREQ--------------- NUMBER OF TIME STEPS BEI WEEN SUCCES. Ree 
TEMPERATURE PRINTOUTS 
СЕХ---------------- ENERGY GENERATION TERM 


Н1,2,3,4,5,6------- HEAT IRANSFER БЕРЕРІ 
S1G1,2.3.4,5.6-----STEFAN--BOLTZMANN CONSTANT (USE 0 IF NO 
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Ke DON iE РЕК СЕ 5.67Е-5 ІР 
УМРИ У ОУ ИВАН ЕР О РЕКПРКЕЗЕХТ ) 
TAMBI.2,3,4,5,6----AMIBIENT TEMPERATURE 
TINIT-------------- ГУР ТЕМРЕЕ ЖСБ ОГАГІГХОРЕ5 
а УАГ----------------- РОО SE ВЕ ИК DETERMINE IF VALIDATION 
Е JST ISTO BERUA 


ПИО ADDITIONAL YARIABLES FOR SATELLITE APPLICATIONS: 


ABS1,2.3.4,5,6-----ABSORPTIVITY 


E ALBCO-------------- ALBEDO COEFFICIENT 
j \LT---------------- ALTITUDE оре В ЛЕГЪШЕЛ КИ О МЕТЕКУ 
EARTIT-------------- Е ТИЕ ИЕ Jili ADSORPIIMITY FOR 


БЕКЕЕВА ИЕ ЖЕУДІ Е А5 ІНЕ 
EMISSIVITY FOR RADIATION HEAT TRANSFER 
D E BOTH ARET THEINFRA RED REGION) 


ECLIPS------------- pP I ODOI LINE J ee ESE SATELLITE 
ISIN THE SHADOWOF THE EARTH 

Г.А----------------- Ci ONS СЕЛО СОКА ЕВЕВО ЕЕЕ 

FE ----------------- ЕО СОВОК ВОВЕ FLULNX 

SAC1.2.3,4.5,6-----SOLAR ASPECT COEFFICIENTS 

SOLAR-------------- SOLAR- FLER 


fem NON-LSER VARIABLES 


* A. GENERAL: 


* B C Ееее COEFFICIENT VECTORS FOR SOLVING THE 
х TEMPERATURE EQUATIONS IN SUBROUTINE TRIDAG 
* PETA... VECTOR OF INTERMEDIATE COEFFICIENTS SOLVED 
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ЗОРА ОДРИ РРА Пау 
ВТОТ1 2,3, 5.6--ВЛОЛ МЕОЗМЕВЕК 


COUNT------------ COUNTER FOR DETERMINING THE NUMBER ЕЕ 
ELAPSED TIME STEPS 

ЕО--------------- РОБ ВЕС MBER 

GAMMA------------ VECTOR OF МГЕ ДЕР ЛЕ ШЕЛФЕРЕ ЛЕГЕ ГЕ 
SOLVE TN SUBROUTINES eG 

IF ,L------------- NUMBERS CORRESPONDING TO IHE FIRST ANDES I 
EQUATIONS TO BE SOLVED IN SUBROUTINE TRIDAS 

QDOT------------- MATRIX OF NODES DEALING ҮҮГІН ЕХЕБС 
GENERATION 

І---------------- MATRIX OF TEMPERATURES AFTER HALF-TIME Sm 
INTHE Z DIREC IMOS 

ТЕМР------------- TEMPORARY STORAGE VECTOR FOR TEMPERATURE. 
CALCULATED TN SUBROUTINE TRIDAG 

TIME------------- DURATION OF CALCULATIONS 

TSTARI----------- MATRIA OF TEMIPERATURES AFTER НАТГЕ- ТП ЧЕ Ле 
IN THE DIREC TIOS 

TSTAR2----------- MATRIX OF TEMPERATURES AFTER HAXLF-TDEVIESSHESS 


INTHE 3 DIRECTION 


. ADDITIONAL YARIABTLTES FOR SATELEE PTICO T 


ALBED1I,2.34.5,6---ALBEDO FLUX RECPIYEDIO EUR CE 


DIST--------------- SUM OF EARTH RADIUS AND SATELLITE ALTIUS 
EARTH1,2,3,4.5,6---EARTH ЕБЕКЕСЕРЕР АЕ S 

M U----------------- GRAVITATIONAL CONSTANT 

PERIOD.------------- THE TIME FOR ONE SATEDEIIL ORDbII 
КЕ----------------- RADIUS OF THE EARTH 

SOLAR1,2,3,4,5,6---SOLAR FLUX RECEIVED OS A SURFACE 
SUN---------------- THE TIME IN MISNCIFS THES ATELLITE ЕВО Е 


SUNLIGHT DURING OXE ORBIT 
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| 


meee VARIABLES FROM EXTERNAL PROGRAMS: 

ЕВ ME-----CSED IX EXTERNAL IMSL ROUTINE FOR DETERMINING 
í AMOUNT OF CPL TIME USED 
- NOUT------ BSED IN EXTERNAL INIST ROL TINE FOR DETERMINING 
E АМО Б NT'OP CPU IIME USED 
ч TIMEO-----CPU TIME AT BEGINNING OF CALCULATIONS 
j TIMEI-----CPU TIME AT END OF CALCULATIONS 
i UMACH-----USED IN EXTERNAL IMSL ROUTINE FOR DETERMINING 
E PON Or CPL TIME USED 


хы мыл лла СЫ е зи Ха За УМ Ул УК ХА ste 04 ote aie ats ait is ote ts ote ots ate ok ste ote ole ake oh За та та 19 


whe ato Se ate <° ate эф» ate ote ote ate we ate ots ale Јо ме бе ко ето чо ate ots whe ale ae ate ate ate «је сре da ato Jo de ate de эф» Фе оо Фе че Фе Фе Фе Фе Фе 1; ate ale $: ado Jo Jo x“ s <” < whe de ale ale обе “o Jo whe ote we ato oa o'e 
обе 7.7 7.7 7,7 27,6 41% 2,7 74% 7,6 7,7 7,7 27 9ұ7 7“ CGF 7,7 247 94% 9,7 2,6 7,7 24% 7,7 4,6 46 8,7 24% 2,6 7ұ6 4,6 7ұ6 7,6 7,6 4% 446 4ұ6 446 4,0 4,6 74% 49% оро оце одо одо оце 4,6 9) 9% 24% ChE 24% 44% Оде еде ефт оце еце ор бро 917 44% 7,9 9% 49% 997 9,0 97 997 


DIMENSION A( 30). B( 30). C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAMMA( 30). QDOT(30, 30, 30) 
IIBSENSION T( 30, 30, 30) TSTARI( 30, 30, 30) 
DIMENSION TSTAR?2 30. 30, 30) 


REAL TINIT. TSTARI, TSTAR2, T. TEMP 

REAL CP, K1, RHO. R. R1, R2. TIME, DELT 

AL LX. LY. LZ DEL DELY, DELZ 

KA p C D. BETA. GAMMA 

LIP SJUTLUXI, FLUX2, FLUX3, FLWX4, FPWX5, FLUX6 
REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMB5, TAMB6 
REAL HI, H2, H3, H4, H5, H6 

REAL BIOI, BIO2, BIO3, BIO4, BIOS, BIO6 

REAL 5101, 5102, 5165, 5164, 5165, 5166 
ЗЕМЛСІГСИ ЖӘКЕН 

REAL V1, V2, V3, V4, V5, V6 

REAL W1, W2, W3, W4, W5, W6 

КЕЛ Хх X23. 889, X5, X6 
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REAL EPS1, EPS2, ЕР53, ЕР54, ЕР55, EPS6 

REAL ODOT CE MH 

REAL СТІМЕ, ТІМЕО, ТІМЕ! 

REAL SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLARI, SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
REAL EARTHI, EARTH2, EARTH3, EARTHA4, EARTHS, EARTH6 
REAL ALBEDI, ALBED2, ALBED3, ALBED4, ALBEDS5, ALBED6 
REAL ABSI, ABS2, ABS3, ABS4, ABSS, ABS6 

REAL SACI, $4€209 30909308459 46509 086 

REAL SUN, ECLIPS, QI, Q2, Q3, PERIOD, PI 

REAL P1, P2, P3, P4, P5, P6 


INTEGER IJ, IS ЕБ РЛАР 
ІМТЕСЕК СОСМТГЕКЕО БОБО ЗИ ЕИ ЕЕ 


ЕЛЕБИ ИЕЛЕ ав 


COMMON TINT, а КУКЕ ЕО ТЕ NID 

COMMON CP. K1, RHO R RI R ИЛЕ ЕРЕ 
COMMON LX; LY, LZ. DELY DEEI TDR 

COSMOS LIL RENS E EEE SIN 

CON NOS A.B C. D. pETA GS I [A GO 

COW МТО ЕТС ХЕ ОЩЕ ЕХО ЛЕВЕ Ве 
СОММОХ ТАМВІ. ТАМВ2, ТАМВЗ, ТАМВ4, ТАМВ5, ТАМВ6 
СОМ МОХ НЕН>, НЗ. НЕЗ, О 

СОММОХВІСІ, ВІО2, ВІО5, ВІО, ВІО5, ВІО6 

COMMON SIGI, 5102, 5103, 5104, SIGS, SIG6 

COM IOS LIT Abo DO 

COMLVIOUN V1, V2, VOV ao 

COMMON WIL W2, WS) Ws Wows 

COMMON AL A2, A3 NIAS A6 

COMMON COUNT, FREQ, ANS1, ANS2, ANS3, VAL, NOL T 
COMMO ЕРБІ, ЕР52ЕР55 ЕРЗ ЕР ЕО 

СОХ ОМ ОБОГО ss 


МОХ ИМЕО, ИМЕ! 

ИНО SOLAR EARTH. RE, DIST, ALT, FE, FA, MU 

COMMON SOLAR], SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
TORINON EAXARPHT, EARTH2, EARTHS, EARTH4, EARTHS, EARTHG 
COMMON ALBED1, ALBED2, ALBED3, ALBED4, ALBED5S, ALBED6 
COMMON ABSI, АВ52, АВ55. АВ54, АВ55, АВ56 

COMMON SACI, SAC2, SAC3, SAC4, SAC5, SAC6 

meri ytON SUN, ECLIPS, QI, Q2, Q5, PERIOD, PI 

SOSTON PI, P2. P3, P4, P5, P6 


X 
Xx e xe Me Me YE sh ok sh ck sh sh ake oh Yk Не ex xe cR o МОМ ЕН НИ КИНЕ МЕНИ ИЕ К КОК УК ЭБ 0 06 СЕ ОО Ж 


h ate ste або st Je 3 We ate вр ве абр 42 4% sto ste n's n'o fo ato ate ote ale x“ ate ate nt, y 4....... У ate al yz ate абр уб 2% 2% «2 «2 у ate ate ate wh ato a". ats ate ate ate s. «> оа atr до g's ate об 
ee mgt 29% Fee Dye 2.9 gs Oye Set 2522 1,М Mgt aye aye Nye Der Seo Peo 55 25 2,1 2,5 25 2,9 35 27 2% 2,5 2,4 2,2,5 2,5 Men Ian ES AS uS IA MS NLIS TALIS IS. SCC age tye See Д 2,5 2,2 Фе tee cee 


Pe IALIZING CPL TIME 


МЕЧКО - СТ! МЕС) 


X X5 siz ste siz Ys Yd ox nob за 36 де Ус и Но Ма 22 За е мо ие зе зе 1: 


КОШО 5 РАКА МЕГЕВК5 TO BE LSLD THROLGHOLT THE PROGRAM 
ВЕ NOW BE READ INTO THE PROGRAM. 


WRITE(*.") ‘ENTER TINIT, THE INITIAL TEMPERATURE.’ 
READ(*,*) TINIT 

WRITE(*,*) (ENTER K1, THE THERMAL CONDUCTIVITY.’ 
READ(*,*) K1 

WRITE(*,*) ENTER CP, THE SPECIFIC HEAT OF THE MATERIAL.’ 


WRITE(*,*) (ENTER RHO, THE MATERIAL DENSITY.’ 
READ(*,*) RHO 

WRITE(*,*) ‘ENTER GEN, THE INITIAL GENERATION TERM.’ 
READ(*.*) GEN 

WRITE(*,*) “ENTER LX, THE LENGTH IN THE X DIRECTION.’ 
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READ(*,? 

WRITE(*, i ‘ENTER LY, THE LENGTH IN THE Y DIRECTION. 

READ(*.*) LY 

WRITE(*,* ња ‘ENTER LZ, THE LENGTH IN THE Z DIRECTION.’ 

READ(*,*) L 

eee ‘ENTER M, THE NUMBER OF NODES IN THE X DIRECTION.’ 
READ(*,*) M 

WRITE(*,*) ‘ENTER NX, THE NUMBER OF NODES IN THE Y DIRECTION. 
READ(*,*) N 

WRITE(*,*) 'ENTER P, THE NUMBER OF NODES IN THE Z DIRECTION” 
READ(*.*) P 

WRITE(*,*) “ENTER DELT, THE TIME INCREMENT IN SECONDS.’ 
READ(*.* | DELT 

WRITE(*,*) “ENTER FREQ, THE NUMBER OF TIME STEPS BETWEEN’ 
WRITE(#.") ‘SUCCESSIVE PRINTOUTS. 

READ(*,*) FREQ 

WRITE(#,*) ‘ENTER ALT, THE SATELLITE ALTITIUDE IN KILOMETERS.’ 
W RITE(*. *) (NOTE: FOR AN EARTH-BASED MODEL, ENTER 0 FOR ALT.’ 
READ(*.* ALT 

WRITE(*.*) "ENTER ECLIPS, THE TIME IN MINUTES THE SATELLITE WILL’ 
W x ‘BE IN THE SHADOW OF THE EARTH DURING AN ORBIT.’ 
WRITE(*,*) ‘NOTE: FOR EARTH-BASED MODELS, ENTER 0 FOR ECLIPS.’ 
READ(*.*) ECLIPS 

WRITE(*,*) 'ENTER FE. THE GEOMETRIC FACTOR FOR EARTH FLUX.’ 
W Mi 'NOTE: FOR GEOSYNCHRONOUS SATELLITES AND' 
WRITE(*,*) EARTH-BASED MODELS, ENTER 0 FOR FE.’ 

READ(*,*) FE 

WRITE(*,*) ‘ENTER FA, THE GEOMETRIC FACTOR FOR ALBEDO FLUX.’ 
WRITE(*,*) ‘NOTE: FOR GEOSYNCHRONOUS SATELLITES AND’ 
WRITE(*,*) 'EARTH-BASED MODELS, ENTER 0 FOR FA.’ 

READ(*,*) FA 

WRITE(*,*) ENTER SOLAR, THE SOLAR FLUX.' 

WRITE(*,*) NOTE: AVERAGE ANNUAL VALUE IS 1353.’ 

WRITE(*,*) 'NOTE: FOR EARTH-BASED MODELS, ENTER 0 FOR SOLAR.’ 
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READ(*.*) SOLAR 
WRITE(*.*) ‘ENTER EARTH, THE THERMAL RADIATIVE FLUX OF EARTH.’ 
WRITE(*,*) NOTE: AVERAGE ANNUAL VALUE IS 237.’ 
DESETE(* S) NOTE: FOR GEOSYNCHRONOUS SATELLITES AND’ 
WRITE(*,*) 'EARTH-BASED MODELS, ENTER 0 FOR EARTH.’ 
рг) клктн 


WRITE(*,*) “NOTE: AVERAGE ANNUAL VALUE IS 0.3.’ 

WRITE(*.*) ‘NOTE: FOR GEOSYNCHRONOUS SATELLITES AND’ 
WRITE(*,*) 'EARTH-BASED MODELS, ENTER 0 FOR ALBCO.’ 
READ(*.*) ALBCO 

WRITE(*,*) ‘ENTER VAL: THIS WILL DETERMINE IF YOU WANT THE’ 
WRITE(*.* У ОЕ SEDATION SURROUIINE CARDED IF YOU WANT: 
WRITE(*,*) ‘THE VALIDATION SUBROUTINE CALLED, EXTER 1 FOR’ 
W RITEC. ПО IF OLR ANSWER TIS O ENTIERO 

READ(*,*) VAL 


Peels FOR LEFT FACE 


WRITE(*.*) ‘ENTER TAMBI, THE AMBIENT TEMP. ON THE LEFT FACE.’ 
READ(*.* *) TAMBI 

WRITE(*,*) ‘ENTER HI. THE LEFT FACE HEAT TRANSFER COEFFICIENT.’ 
READ(*,*) H1 

WRITE(*.*) “ENTER FLUX1, THE FLUX ON THE LEFT FACE.’ 

READ(*.*) FLUXI 

WRITE(*,*) ‘ENTER SIG1, THE LEFT FACE STEFAN-BOLTZMANN’ 
WRITE(*,*) ‘CONSTANT. NOTE: IF RADIATION NOT INVOLVED, 
WRITE(*,*) ‘ENTER 0 FOR SIGI.’ 


READ(*,*) SIGI 
WRITE(*.*) ‘ENTER EPS1, THE LEFT FACE EMISSIVITY.’ 
READ(*,*) EPSI 


ARNE) ESTER ABSI, THE LEFI FACE ABSORPTIVITY. 
WRITE(*,*) NOTE: FOR EARHT-BASED MODELS. ENTER 0 FOR ABSI.’ 
READ(*,*) ABSI 


READ(*,*) SACI 


INPUTS FOR RIGHI FACE 


WRITE(*,*) ‘ENTER TAMB2, THE AMBIENT TEMP. ON THE ’ 
WRITE(*,*) ‘RIGHT FACE.’ 

READ(*,*) TAMB2 

WRITE(*.*) 'ENTER H2, THE RIGHT FACE HEAT TRANSFER' 
WRITE(?,*) ‘COEFFICIENT.’ 

READ(*,*) H2 


УКЛЕС,У) ENTER SIG2°5THE RIGHT PACE STEFAN-BOEIZA L m 
МИТЕ АЗ ГААТ 

WRITE(*.*) “NOTE: IF RADIATION NOT INVOLVED, ENTER 0 БОҚ em 
RES D) SIC? 

WRIITEG) ENIER EPS2 IHE RIGIII II МЕ ЕЛЕК ИК 


WRITE(*™.*) ENTER ABS2. THE RIGHT FACE ABSORP Ia НА 
WRITE(*,*) “NOTE: FOR EARTH-BASED MODELS, ENTER 0 rOR A L m 


WRITE(*.*) ‘EXTER SAC2, THE RIGHT FACE SOLAR ASPECT’ 
WRITE(*,*) ‘COEFFICIENT NOTE: FOR EARTH-BASED MODELS.’ 
WRITE(*,*) ‘ENTER 0 FOR SAC2.’ 

READ(*,*) SAC2 


INPUTS FOR BO} TOM FACE 


WRITE(*,*) ENTER H3, [THE BO] TOMI EAGE EAT TRANS eee 
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WRITE(*,*) ‘COEFFICIENT.’ 

READ(*,*) НЗ 

WRITE(*,*) ‘ENTER FLUX3, THE FLUX OX THE BOTTOM FACE.’ 
READ(*,*) FLUX3 

WRITE(*,*) ‘ENTER SIG3, THE BOTTOM FACE STEFAN-BOLTZMANN’ 
WRITE(*,*) ‘CONSTANT.’ 

WRITE(*,*) NOTE: IF RADIATION NOT INVOLVED, ENTER 0 FOR SIG3.’ 
READ(*,*) SIG3 

WRITE(*,*) 'ENTER EPS3, THE BOTTOM FACE EMISSIVITY.’ 
READ(*,*) EPS3 

WRITE(*,*) ‘ENTER ABS3. THE BOTTOM FACE ABSORPTIVITY.’ 
WRITE(*.*) 'NOTE: FOR EARTH-BASED MODELS, ENTER 0 FOR ABS3.' 
READ(*.*) ABS3 


WRITE(*.*) COEFFICIENT. NOTE: FOR EARTH-BASED MODELS’ 
РКІТІ(” 2) ENTER О ГОК SAC3.’ 
REA D(*.*) SACS 


meee TS FOR THE FOP FACE 


WRITE(*,*) ‘(ENTER TAMB4, THE AMBIENT TEMP. ON THE TOP FACE.’ 
READ(*.*) TAMB4 
WRITE(*,*) 'ENTER H4, THE TOP FACE HEAT TRANSFER COEFFICIENT.’ 


WRITE(*,*) 'ENTER 0 FOR SIG4.' 

READ(*,*) SIG4 

WRITE(*,*) ENTER EPS4, THE TOP FACE EMISSIVITY.’ 
READ(*,*) EPS4 


ТШ 


READ(*.*) ABSJ 
WRITE(*,") ENTER SACA. THE TOP FACE SOLAR ASPECGTIGOEPEIGIESSM 


КЕАШ С S 2€4 


INPUTS FOR THE RO ИН ИТ 


WRITE(*,*) “ENTER TAMBS, THE AMBIENT TEMP. ON THE FRONT FACE.’ 
READ(*,*) TAMBS 


WRITE(™,*) ENTER FEUAS, TRE FEC S O Hi Во ЕЕЕ 
READ ЕЕ 

WRITE(*,*) ENTER SIG5, THE FRONT FACE STEFAN BOLCTATSSS 

OR PECs SG ONS VAN. 

WRITE(*,*) “NOTE: IF RADIATION NOT INVOLVED, ENTER O FOT ILE S 
READ GG] о то 

WRITE(*.*) ENTER ЕР55, ЕНЕ FRONT RAGS ЕЕ ae 


WRITE(*,7) ENTER ABSS, THE ЕКОХЕ РАСЕ АББОК РАКИЋ 
WRITE(*,*) “NOTE: FOR EARTH-BASED MODELS, ENTERS Oke == 


WRITE(™,”) ENTER SAC5 THI FRONI Е АСЕ ОЕ S RP Lt G 
WRITE(’,*) СОЕРЫСТЕХТ. ХОЗЕЗЕФВЕАВ ЕВЕ ОЕ В 
WRES ESTEK OPFOR о 

Е 


ISPLUISFORTHE BACK FACE 


WRITE(*,*) “ENTER TAMB6, THE AMBIENT TEME. ON THE BACHI ЕВ 
READ(*,*) TAMB6 


READ H6 
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ша ТШЕ E T ENTER FCC THE FLUX ON THE BACK FACE.’ 
ЈЕЛО“) FLUNG 


WRITE(*,*) “NOTE: IF RADIATION NOT INVOLVED, ENTER 0 FOR SIG6.' 
READ(*,*) SIG6 

WRITE(*,*) ‘ENTER EPS6, THE BACK FACE EMISSIVITY.’ 

READ(*,*) EPS6 

WRITE(*,*) “ENTER ABS6, THE BACK FACE ABSORPTIVITY.’ 

WRITE(*,*) “NOTE: FOR EARTH-BASED MODELS, ENTER 0 FOR ABS6.’ 
READ(*,*) АВ56 

WRITE(*.*) 'ENTER SAC6, THE BACK FACE SOLAR ASPECT COEFFICIENT.’ 
WRITE(*.*) 'NOTE: FOR EARTH-BASED MODELS, ENTER 0 FOR SAC6.' 
READ(*.*) SACO 


. " e*s 3°, об. чб абр обр чб абр абр эъ, ut, x бо аб абр абр абр абр абр або абр 2% ate ~ На ats абр абр або абр абр об > ss w^ sl oh ate 2% abl ab хх ah м» x“ +, +, У ~ >” E ~ whe ol у «Ӧн або nto sx 0% st 
7 9,% 6° ... wa? LII ~" "4" M id „47 2,” ve" ги 2... 6" ter 2” АГЫ ar ми Ми а” ` .. .. 4. .. "4 wae ~" 5." ~" + ГЫ wy? ~." ... 2” .. и .. ми `. war ч“ ~" ЫГЫ 4. „4. а” wae 4” .. 047 2”. ..” .." oe Tar dL we" 


о ahe e da ale ale ole ale ate ъде we ate або ale We ж/е ее de te Se «је је de he de ae Ye Ye do Se ate е ole we ъбе we <° af n о ао һб ђе w eo «To оле ate ate ate ци о йо мо ake ale we «fe а/е ste da as 
x 7 4 de а, За ate we ate ale ale We we Se «је ate Se «је ale је he de ake we we «је е ае ае «е «йо ale we ahs ale аҹ ъфе + wJe We ate Se ate ahe ale ale абе aho So Yo Аве ate ado we sP. We ake ws 0° 
% 94% 44% %%% 7;% 44% 7,% #49 74% 949 949 #65 99 969 905 94% 995 949 ngn yy FEF дэ гыч 64% 6,6 одо еде ода еде еј 74% 9)67 7;% 79% 74% 969 949 Pam до одо еде едо 94% 909 PES 6,5 #49 949 969 969 249 гот Pyn Sgr So? 945 95 949 649 одо дъ оце еде 


meee CO LATIONS OF CONSTANTS TO BE ESED THROLGHOLT 
THE PROGRAM 


IDEL NX 
DELY 
DELZ 


LX (M-I) 
р 
LZ (P-1) 


R = DELX (DELY*DELZ) 
Ri = (DELX**2) (DELY**2) 


ALFA = KL(RHO*CP) 
FO = ALFA*(DELT,2)(DELX**2) 


BIOI = HI KI*DELX 
ВНК "МЕТА 
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BIOS 
BIOS 
BIOS 
BIOG 


EARTHI 


= H3 KI*DELX 
= H4 KI*DELX 
= НТО ЮЕ 
- H6 KI*DELX 


RI*FO 
R2*FO 


ate 
*.* 


sta 
2. 


to r3 tJ 


ә ко м 


tJ 


ty 


RI*FO 
R2*FO 


*FO 


= EPS] 


EARTH? = EPS? 


EARTH3 = EPS3* 
EARTH4 = EPS4* 
EARTHS = ЕР55* 
EARTH6 = EPS6* 


QI = 2*FO* 


*FE 
*FE 
FE 
FE 
FE 
FE 


*FO*BIOI*SQRT(RI) 
'FO*BIO2*SQRT(R1) 
TO PI SOD] po) 
*FO*BIO4*SQRT(R2) 
*FO*BIOS 
= 2*FO*BI106 


БӘСІРЕ 
ЕР 
“EARTH 
"EARI 
ШӘРІП 
ЕАСИ 


RI*DELY КІ 


Оо РОВЕР В 
Q3 2 2*FO*DELX'KI 


P4 = 
1 


СО? 
ФЕ 


(Ti а 
Евы 


EARTHI) 
EARTH2) 
EARTH3) 
EARTH) 
EARTHS} 


P6 = Q3“(FLUX6 + EARTHó) 


КУ -2°FO*SIG]*EPS1*R1*DELY;K1 
КО = -22FO”SIG2*“EPS2*RI*“DELY'K1 
fee -2"FO*SIGO*EPSS"R2*DELZ'K1 
X4 = -2*FO*SIGA*EPSA*R2*DELZ/K1 
Es -2*FO*SIGS5S*EPS5*DELX:' KI 
Хб = -2*FO*SIG6*EPS6*DELX'K] 


КЕ = FO"(DELX**2) K] 


ALBEDI = ABSI*FA*"SOLAR*ALBCO 
ЕВЕ. = ABS2^rFA"SOLAR*ALBCO 
Ж ЕРЕПШЗ - АВЗЗ РА"ЗОГАК"АГВСО 
ALBED4 = ABS4*FA*SOLAR*ALBCO 
Pepe DS = ABSS*FA*SOLAR*ALBCO 
ALBED6 = ABS6*FA*SOLAR*ALBCO 


SOLARI » SOLAR*SACI"*ABSI 
ОГАН) - SOLAR"SAC2*ABS?2 
ПСЕ = SOLAR"SACS"ABS3 
SOLAR-4 = SOLAR *SAC4”*ABS4 
КОШ АКУ = SOLAR*SACS*ABS3S 
КОШ АКО = SOLAR*SACO"ABSO 


WI = Pl + QI*(SOLARI + ALBED!) 
W2 = P2 * QI*(SOLAR2 + ALBED2) 
W3 = P3 + Q2*(SOLAR3 + ALBED3) 
W4 = P4 + Q2*(SOLAR4 + ALBED4) 
W5 = P5 + Q3*°(SOLARS + ALBEDS) 
W6 = P6 + Q3*(SOLAR6 + ALBED6) 


Seer An  ih-bAseD MODELS, CALCULATIONS FOR PERIOD AND SUN 


- sea RE NOt REO REI: 


IF (ALI IE C O TOS 


PI = ATAN(1.0)*4.0 
МС = 398603.2 

ВЕ = 6378.165 

DIST = RE + ALT 
ECLIPS = ECLIPS*60 


ate ot whe ate x <Ó EET У ate ам ам аа ам ан АМ ам ам аи аи «мо «и Ми «и «и Ми аби аб и обо 409 
тул за За С, Са 264 $55 245 29% 144 4,55, 2,1 1,5 1,1 1,4 2,4 1,42, 2,% 1,5 2,4,4 2,5 1,1 541 544 2,2 Mgt tee Set Mgt et Зее Mgt Sgt Sgt ngs age nye Mgt Sgt sgt gt Set Ee Sgt ТАЗИ ТА Та ТА Та ТА ТА МИНИ ЛЕ 21 


в, * е, 4 е, е, %, е, 9, е. U °. 4, . 4, Q + %, +, . 9, е. 4. % е, е. е (ЗА U 9. е е +, t ` «t t, . 4, U е $, + е . $, 4, + . 4 +, +, Ф, + + + + +, е е + 4 + е. 
ооо 2 ааа с заета с а қақ уд Уа ақа саса А а ағы 


INITIALIZING THE HIME AN DIGG IEE 


> TIME = 0.0 
COUNT = 0 


OPENING FILES FORPRINT SITATEN TSI. 
СРЕ ОС =- FELES DATAT 
OPEN(UNIT = 2, FILE ='’DATA?') 
OPEN(UNIT = 3, FILE =’DATA3’) 
OPEN(UNTT =4, FILE=‘DATA4) 
OPEN(UAIT=5 FLES DATA) 


аа Заа. 


жасы ТИ ТИ ПУТА МИ ТИ ТА РА ТА ЖЕ ТЕ А МА сд си te ts ha ta os he ds oe ots he а Ме Ма Ма За ИА Ма a a oe os oe oe ae eo oe te fe he eo ae Ин gO SK ЗА Мо ПАОЛА 


~“ 
` ry ares 444 54% ера ча а 44. 





s TTIALIZE ALL TEMPERATURES TO TINIT, THE INITIAL 
Е IPERATURE. 


DO I01=1, M 
DO 20 J=1, A 
POSOK=I PF 


НИК) = ПАГ, 


30 GON TINCE 
рио ОК ПАСЕ 
ШСО-ТІХІЕ 


$, 4, 4 4 4 $, е, 4 4 г) * + $. J е $. е. я е + ^, $. . 4 е $, % . а $, . “4 е p 4, т 4. е $, 4 $, $ е в 4, 
75 2% 22 922 е 14 Ма tte Зи Зе и За ЗМ ОЗИ Зи Sede mx re cR E Mer e ge xe xe xx ys ok yo 0 З Мо за та За Ус Јо Ја Ке зи за Зи 204 Уа зи зк Уд зи. За Жл эл Ул Уз Уз су а Ја џез 


ate ote ale «Mo a". s) w° ate ote ate оба sta ate ate als abe ats ate whe wile ate 4%. «Уз «Әә 47 «би 47% ale ate ate ats ә)» ats ale ale als als ato whe 479 ale ate atv ats ate ats oe we ats ale we ole whe ale Јо Sw ote ate abs ale ale ote abe at де ate ote ate o she abe 
565 545 949 9,9 9,9 545 949 9,9 949 994 %4% %)% %)% 9/% 949 Bye %;% 76% %)% %;% %4% 505 505 545 949 949 9,9 9,9 949 949 949 949 9,9 949 905 SEF 9,9 949 9)» 949 949 94% %/% %/% 249 949 %/% %;% 9,% %)% %)% 9% %;% %/% %)% %)% Og одо одо рде объ одо одо одо одо 54% 905 545 349 54. 


ЖІШІСТІЛЕ АІТ СЕ-ЕКАТТОХ ТЕКМӘ ТО GEN, THE INITIAL 
= GENERATION. 


DO 40 I= L.M 
DO 50J2 LN 
DO 60 K- LP 
ODOT(IJ.K) » GEN 


60 SON TINUE 
ВО CONTINUE 
40 CONTINUE 


whe 

"~ 

Ж: :& эд :& Ж Ж: Ж: с ch Sh ah ok oe Se УМ ЗА Me oe oe ok dod c X X Y Y dede x Hh Me fe He Me fc fc Ж Жижић 
` Gd + е 7 ~ ` ` ` ` ` ` ` ` * ws + + ` ` ` - + ` ` ` ` Ф Ф + + + + - ` 

mx А ден о Seu Sos нижим See Set wat St MS He Me Be Se ee dle sh he ode Me okt oe oe ks ese ch SK oe Se cis te He ch ee as hs ee ke EP HE So Ee Ke KE HE 


x 


EN CIE VENTENG THE TIME AND COL NTER 


70 TIME = Ti i Eso ee 
COLAT = COL 


ete мо мо обр абр об оба ойр ols ob ate ate ate of aio st до мо пва мо ча s. об. мо ово абе Ми Фо «а а o's oh 492 Y ate als ale ale whe Wh st, o uu ate ale So ve Ve oe Se 3% чо ге Jis z's ya ghe Ae s'e a'a ale xk Де зе за Хе +. уа мо she oale 
145525552211 2,1,2 22 5225 22 222,2 2,2 2,5 2,9 2,2 2,2 2,2 5,2 2,2 2% 222,525 2,5 2,2 2,2 age aye Ife aye age age age 269 Фра го Зе Tao tae tuo Ще Сао ща aye Fae Tee age age Dao 2,2 2,5 2,5 2,5 Hye 208 206 S Ta Dao 22222,2 2,2246 242 245 


We Ae Me He Ae Ye she she Ne ae he Но а по а Ке ТИ Re he he oe Ne oe oe ee Me ake ЗА Ще Зи Me ae ae Re ae ae ae fe Re гм К Pk 2k 25 55 Ree Me ae Re Me Ke HEME ee ee SE Ж he ee oh ee ME He 22 


~ 
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* THISPART OF THE PROGRAM CALCULATES THE ЗХОГАК ЕГСХЕЗЕ ИЕ 
" ALBEDO FLUXES WHICH WILL CHANGE WITH TIME AS THE SATELLITE 
* “REVOLVES AROUND THE EARTH. DURING ECLIPSE PERIODS 
* THE SOLAR FLUX WILL BE ZERO. FOR EARTH-BASED MODELS, 
* THIS PART OF THE PROGRAM AND THE SCBROCUTINESPERI A p 
“ PER? WILL BE BYPASSED USING AN Ir HI v SI A PP NI IR 
ТНЕ ЗЕКТАРГТЕОРОКА ТЕ БАК ПАЈА РАВНИМ ОЮ 
THEN THIS PART OF THE PROGRAM AND THE Tw SUBGROUP. 
МУДАР О ПРИ а зе 


ЕЕС ЕЕЕ ОТО ИВ 


ПЕСМЕ Т ООРОО > 
CATO LES 

ES 
CALL PER? 

EADIF 


ALBEDI = ABSI*FA#*SOLAR*ALBCO 
ALBED? = ABS2“FA“SOLAR*ALBCO 
ALBED3 = ABS3*FA*SOLAR*ALBCO 
ALBEDS = ABSS*FA*SOLAR*ALBCO 
ALBED5 = ABSS*FA“SOLAR*ALBCO 
ALBEDO = ABS6*FA*SOLAR*ALBCO 


SOLARI 2 SOLAR*SACI*ABSI 
SOLAR2 = SOLAR SAC 52 


SOLARS = SOLAR*SAC3*ABS3 
SOLAR4 = SOLAR*SAC4*ABS4 
BOARS = SOLAR*SACS*ABSS 
SOLARG6 = SOLAR*SACG6*ABS6 


WI 
W2 
W3 
W4 
Ws 
WO 


PI 4 QI*(SOLARI 
P2 + QI*(SOLAR2 
P3 + Q2*(SOLAR3 
РА + Q2*(SOLAR4 
P5 + Q3*(SOLARS 
P6 + Q3*(SOLARG 


АБВЕР) 
+ ALBED2) 
AT BEDS) 
+ ALBED4) 
+ ALBED5) 
+ ALBED6) 


4........, У ate ote x° x als ale ats «ба s > лаззат ste во е мо УУ 72 h ай, бо >> КА 
z⁄ `; 3⁄2: 245 Sus ye See Teh то 2⁄2 2⁄2 2⁄5 2⁄2 2⁄4 $ тии те те ди ЗА ЗА ЗА ТА За За За 2⁄2 + 


> .,........,...-.,........ MEREVEE VEYE E Ey , 46, 0 
в а о воть 


ESNHUNG COEFFICIENT ARRAYS A, B, AND € IN THE X DIRECTION 
75 DO $0121, M 

B(D = (1 + US) 

IF(1.EQ. 1) THEN 


А(1) = 0.0 
C(I) = -U5 








ELSE IF((I .GT. 1) .AND. (1 .LT. M)) THEN 
А() = -FO 
C(I) = -FO 





ELSE 
A(I) = -US 
C(I) = 0.0 





ENDIF 
80° CON TENGE 


SETTING COEFFICIENT ARRAY DISA THE X DIRECTION 


DO 90 K=1,P 
DO 100 J=1,N 
DO 110 I= 1, M 


О ДУО АЦЕ s 
CALL NRAYI 
ELSE JTF (KG is) ND (I ies 


CA ee 
ELSE 

ое 
END li 


110 CONTINUE 


IF=] 
pon 


CALL TRIDAG 
DO 120 I=IF,L 
TSTARI(I,J,K) = TEMP(I) 
120 CONTINUE 


100 CONTINUE 
IOTC OSTIS CE 


ж 


ұлы Кин of he oe ote oe ok ake ah 2% os oh ok УИ ok ok да Ум си зи Уго За а 9: Ж: ok ok ok ok зи та 








Band ft bata aS DOS ESQ UM EEG XO E ЗД л л ЕГО ЛИЛ К ЫК ИЧ: ОЛУУ А? 


SETTING COEFFICIENT ARRAYS A, B, AND C IN THE Y DIRECTION 
DO 130 J=1,N 
BU) - (1 + U3) 
IF(J .EQ. 1) THEN 
A(J) = 0.0 
С(Ј) = -U3 
ВЕ ЕЕ СТ. 1) .AND. (J .LT. N) THEN 


A(J) = -UI 
СО) + -L1 


ELSE 
A(I) = -U3 
С(Ј) = 0.0 


ESVIDTE 


EUSCONIINUE 


pet NG COEFFICIENT ARRAY D IN THE Y DIRECTION 


DO 140 K=1,P 
DO 150 I2 1,M 
DO 160 J= LN 


IF(K.EQ.1) THEN 
БЕНЕН 
ЕЕ ТЕ (KOIL D. AND (KLT P) ТНЕХ 


CALL YANK2 
ELSE 

CALL YANK3 
Бхр 


160 GONTINGE 
IF=1 
L = N 


CALL TRIDAG 

DO 170 J= IF, L 

TSTARXIJ.K) = TEMP(J) 
170 CO JI UE 


150 CONTINTE 
140 CONTINUE 


* ate до ае «е зе ato ale Yo де або Vo «ве Yo We Аве «е 4» ще «а «е «ве «е «ве аа «е «е ма «а Jo «ве «ве «ве whe ha Ја ме ме d 
= 49% 909 947 Hae бұз 44% 44% 74% 7)% за“ Рап ca оде бре HES FES FPS TEs бро бро поч бро без mq бич ода ој FES бро 66% 74% еда 74% 60% едо едо е 


SETTING COEFFICIENT ARRAYS A, B. AND C IN THE Z DIRECTION 
DO 180 K=1,P 
B(K) = (1 + U4) 
IF(K .EQ. 1) THEN 
A(K) = 0.0 


UG) = Ue 


ЕСЕ ТЕКСТ.) АОК ИКО ИЕ 


КИК) = -[2 


eon) = -U2 
ELSE 

A(K) = -U4 

С(К) = 0.0 
ENDIF 


lO CONTINUE 


x 


BESUU SG COEFFICIENT ARRAY D EIN THE Z DIRECTION 


ЕШ! |90 1-1,Х 
[O 000 [= ].M 
ШО О О К = Ј|.Р 


IF(J.EQ.1) THEN 
ШЕТІ ZULU] 

ELSE IF((J.GT.1) AND. (J.LT.N)) THEN 
ШТІ 711122 

ELSE 
ЕСІП 71123 

ENDIF 


ate 
"въ 


210 CONTINUE 


ote 
2% 


NH] 
Ег 


CALL TRIDAG 
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x 


220 
200 
190 


ж 


DO 220 K = IF, L 
T(LJ,K) » TEMP(K) 
CONTINUE 


COE ТЕ 
CONTI CE 


и и 8 Ы ИА 09 КОРА КО 0 Д а Жо то За ake ak ak ake a ok ai ok Ал ТЕ SE ии: 


Me Ne she Me oe Me sh sie ke oe sh aks ote Но Не Не «o6 xxx Кс зк СХ К R Ik XZ Zk Se oh sh Me oie Se fk she 5 Ж he oh oR oe ak oe ok ak ok ok У oe e si oe % she oe oe 76 26 0020 


3e 
7“ 


22 


240 


THIS PART OF PROGRAM IS CUSED TO DETER ШЕ ИЕ ИШЕ 
SUBROUTINE УАЗ ВЕ САТЕЕО В ИСЕ 

IS LSED TO VALIDATE THE PROGRAM BRIAN AND COMPARE 
TEMPERATURE DIFFERENCES AT SELECTED NODES ИШЕ 
DIFFERENCES FOUND BY THE PROGRAMIS EXPLICIT ANE ТИ 


IF SUBROUTINE VALID IS USED THE SUBROUTINE OUIPUT ША є 
By PASSED. IT S!IOULD BE NOTED THAT THE SUBROU HINES wee 
WILL NOT NORMALE ВЕСЕ 


IF(VAL .NE. 1) GO TO 250 
IF(TIME .EQ. DELT) CALL VALID 
IF(COUNT .NE. FREQ) GO TO 240 
САБЫ МУАШ 


IF(TIME .LE. 3600)СО ТО 70 


GO TO 310 


THIS PART OF PROGRAM WILL DETERMINE IF THE TIME AND/OR 
COUNT CRITERIA FOR CALLING THE SUBROUTINE OUTPY Tiss ies 
THIS SUBROUTINE IS USED РОК ДО ЩО ОЕ 


130 





ии И МЕ ЕО. ОБЪТ) САЩ, ОСТРСТ 
ШСОГХТ.ХЕ. ҒКЕО) СО ТО 500 


Sr OLTPUT 


* NOTE: IF THE SUBROUTINE OUTPUT IS CALLED, THE COUNT WILL BE 


% Mme TTAEIZED IO ZERO TS IHE SUBROLTINE OUTPLT 


EN IE: THE TIME PERIOD SELECTED FOR THIS PROBLEM IS 86400 


5 pee DS ORDA ANOTHER TIME LIMIT CAN BE USED 
Í ПОЛОВИ РА Se RODESIRES A RERIOD DIFFERENT PHAN ] DAY. 


300 IF(TIME.LE. $6400) GO TO 70 


м» 49, ат, 49, <° 470 27 4% abo <J, <" «де А, «бо Абе ае Ае АТ. s”, А ale ake x" a's ate ate 4) te ае аи «е Wo a's <Ó, «Ло «до «Че бе +" ote whe ale ate Jo ate ОУ ИЗ ote we wie ate а е ме ае Айо мо че е а“ пле ale a's oe 
94% 03% 7% 74% 59% Пра 905 249 25 55 245 75 зр5 75 зуз 75 205 905 905 зуб Рдз Руб Руз 945 7,9 %9% 7;% 7% 89" 7)% 0;0 09% %9;% %;% 04% *q* 07,4 04% 00% пре Фр 04% ПА еда Фра еда Ре 74% 79% 60% 74% 7,% 74% 79% 79% 64% 644 74% 7ұ% 7,% 69% 64% 64% 74% Fy 44% 64% 6% ғұз 


meee LATING AND PRINTING CPL TIME USED DURING THE 
SEERATING OF THIS PROGRAM 


Poe PRINT* 
БАТ“ 


ШЕГЕСІ = СПМЕО 
ПРУ LMACH(2.NOLT) 
Smee SNOUT”) CPU TIME (SECONDS) =, TIMEI-TIMEO 


STOP 
END 


` и ~. ~. ~ КАС; 
ИИО ЖНИВ ЕЕС ЕЕЕ ЕЕЕ ЕСИ 
ж 


SCEROLULTISE PERI 
IDEN SION AC 30). BC 30),.GG30), D( 50), TENP( 30) 


[51 


DIMENSION BETA( 30), GAMMA( 30), QDOT(30, 30, 50) 
DIMENSION I( 30, 350730) TSTAR ITON OS ци 
DIMENSION TS TAR вил вет бу 


REAL TINIT, TSTARI, TSTAR2, T, TEMP 

REAL CP, K1, RHO, R, R1, R2, TIME, DELT 

REAL LX, LY, LZ, DELX. DELY, DELZ 

REAL A, B, C, D. BETA, GAMMA 

REAL FLUXI. FLUX2, FLUX3, FLUXA, FLUX5, FLUX6 

REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMB5, TAMB6 
REAL HI, H2, H3, H4, H5, H6 

REAL BIOI, BIO2, BIO3, BIOJ. BIOS, BIO6 

REAL SIGI, SIG2, SIG3, 5104, 5105, 5106 
РЕ. 

Ял 6 

REATO NS VUE 

REAL Nl, Noe eo 

REAL EPS}. EPS2. EPS3, EPS4, EPS5, EPS6 

REAL QDOT. GEN, YI 

REAL CTIME. TIMEO. TIMEI 

REAL SOLAR, EARTH, ALBCO, RE. DIST, ALT. FE. FA, MU 
REAL SOLARI. SOLAR?. SOLAR3. SOLAR4, SOLARS, SOLARS 
REAL EARTHI, EARTH2, EARTH3, EARTH4. EARTHS, EARTHG 
REAL ALBEDI. ALBED2. ALBED3, ALBED4, ALBEDS, ALBED6 
REAL АВ51. АВ52, ABS3, АВ54, АВ55, АВ56 

REAL SACI, SAC? SAGs SAG nS AGNES 

REAL SUN, ECLIPS, О1, Q2, Q3, PERIOD, PI 

REAL PI, P2, P3, РА, Р5, Рб 


INTEGER I. J K NL РС ПЕ oD oli Tl ео 
ІХТЕСЕК СОГАТУЕКЕО, АЗОТ АЗСА ЕЕ 


COM MON TIN T. ISESRI TSIAR? T TENE 
CON MON CP KI PHO PF pi Si E ЕТ 
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COMMON LX, LY, LZ. DELX, DELY, DELZ 

MESON L], K. M, N, P, IF. L. IFPI, LAST 

COMMON A. B, C, D. BETA, GAMMA. G, Q 

@OMMON FLUNI, FLUX2, FLUX3, FLUX4, FLUXS, FLUX6 
COMMON TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMB6 
COMMON HI, H2, H3, H4, H5, H6 

COMMON BIOI, BIO2, BIO3, BIO4, BIOS, BIO6 

СОММОХ 5101, 8102, 5103, 5104, 8105, $166 

ШОММОХ UI, U2, U3, U4, U5 

COMMON VI, V2, V3, V4, V5, V6 

COMMON W1, W2, W3. W4, W5, W6 

BONDMON XI, X2, X2. X4 X5, X6 

ВАО COLNT, FREQ, ANSI. ANS2, ANS3, VAL, NOUT 
COMMON EPSI, ЕР$2. EPS3, EPSJ, EPS5, ЕР56 

most ON ODOT, GEN, ҮІ 

COMMON TIMEO. TIME! 

COMMON SOLAR, EARTH. RE, DIST, ALT, FE. FA, MU 
COMMON SOLARI, SOLAR?. SOLAR3, SOLAR4, SOLARS, SOLAR6 
COMMON EARTHI, EARTH?2, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI, ALBED2, ALBED3. ALBED4, ALBEDS, ALBED6 
COMMON ABSI. ABS2, ABS3. ABS4, АВ55, АВ56 

On! ON SACI. SAC2, SAC3. SAC4, SACS, SAC6 

COMMON SUN, ECLIPS, Q1, Q2, Q3. PERIOD, PI 

COMMON PI, P2, P3, P, P5, P6 


ETIME. LT. SUN) THEN 
GO TO 10 


ELSE IF((TIME .GE. SUN) .AND. (TIME .LT. (SUN+ECLIPS))) THEN 
GO TO 20 


БЕ РА РИМ СЕ. (50 Х ЕСІІР5)) А ХЮ. 
БИ ПАШЕ (SL N+ ECLIPS Y THESN 
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СОТО10 


ELSE IF((TIME .GE. (2?SUN + ECLIPS)) .AND. 
& (TIME .LT. Q*(SUN+ECLIPS)))) THEN 
GO TO 20 


ELSE IF((TIME .GE. (2*(SUN 4 ECLIPS))) .AND. 
& (TIME .LT. (3*SUN + 2*ECLIPS))) THEN 
GO TO 10 


ELSE IF((TIME-:GE-( SUN РОЗЕОИ ТЕС 
& (TIME LI БЕСТЕ D hi 
GO TO 20 


РЕЗЕ ПРОСТИМ „СЕ. (385 РО КАЕ 
& (TIME LT. @*SUN + 3°ECLIPS))) THEN 
GO TO 10 


LLSE IF TIME GE. (F SUN F S ECEE D 
& (TIME LTP ЕР ЕТ ЕАС 
GO TO 20 


ELSE IP(( TIME СЕ. 455 БЕРЕТ SIM 
& (TIME Li (SUN Б Б ИУ: 
GO TOTO 


ELSE IF((TIME .GE. (5*SUN +4*ECLIPS)) .AND. 
& (TIME.LT. (5*(SUN+ECLIPS)))) THEN 
COTO 


ELSE IF(TIME .GE. (5(SUN-- ECLIPS))) .AND. 


& (TIME ЕТО КО ЗОВУ АВСТРИИ) АИ ` 
GO TO 10 
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ШЕ ЛЕ (ТЕ СЕ. (6“5( + ЭЗЕСГІР85)). АХ. 
РОМЕ LT. (6% SUN +ECLIPS)))} THEN 
Eo ТО 20 


ELSE IF((TIME .GE. (6*(SUN 4 ECLIPS))) .AND. 
& (TIME .LT. (7?SUN * 6*ECLIPS)) THEN 
GO TO 10 


ELSE IF((TIME .GE. (7#SUN + 6*ECLIPS)) .AND. 
& (TIME .LT. (7*(SUN+ECLIPS)))) THEN 
GO TO 20 


ВЕ ТЕМЕ «ОЕ. (“ОС + ECLIPS))) . AND. 
ПРИ ЛР LI (SSSU`N + 7“ECLIPS)) THEN 
GO TO 10 


ELSE IF((TIME .GE. (8®SUN + 7#ECLIPS)) .AND. 
& (TIME .LT. (S*(SUN+ECLIPS)))) THEN 
GO TO 20 


ELSE IF((TIME .GE. ($*(SUN + ECLIPS))) .AND. 
& (TIME .LT. (9*SUN + 8*ECLIPS))) THEN 
GO TO 10 


mee IEQ TIME GE. (9°SUN + 8*ECLIPS))-AND. 
Е Е ГГ. (905 СХ Е ЕСЕ1Р5}))}) THEN 
СО ТО 20 


ELSE IF((TIME .GE. (9%(SUN + ECLIPS))) .AND. 
& (TIME.LT. (10*SUN + 9*ECLIPS))) THEN 


GO TO 10 


Е live GE.(l0eSUN + 9*ECLIPS)) -AND. 
cm UE PICO (SUNSECLIPS))) THEN 
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GO TO 20 


ELSE IF(TIME .GE. (10*(SUN+ ECLIPS))) AND. 
& (TIME.LT. (11*SUN + 10*ECLIPS)) THEN 
COTONE 


ELSE IF(TIME .GE. (11*SUN + 10*ECLIPS)) .AND. 
& (TIME.LT. (11°(SUN+ECLIPS)))) THEN 
GO TO 20 


ELSE IF( TIME ово ПРУ ВБ ВИ АИ ЈАНО 
б (ТЕМЕ. ЕТО ОМО C pir ance 
GO TO 10 


ELSE 
GO TO 20 


Е ЮЕ 
10 ЗОГАК - 1353 


WRITE(6,") ‘SOLAR 
GO TO 30 


SOLAR 


|| 


ДБО С) 
WRITE(O-*) SOLAR = ӘГАК 


КЕЛИН 
EXD 
ини ДЖ f ts te t t te ot fe ots ot ol te ok ae ok ak ake > Ж эй ak ak gai ak ot ts c doc ote s Ж: Ж Ж ak d d aoo oco dei Xo оу 


ЖОАО МОРДА МОКОО МАМОН МАИ МАМА АДОИ О ДОО МАМИНОЙ ИИМИ ЛА ЖА ЖООК КК А А А 


SUBROLTINE PER? 
DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
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ШГТЕХЗІОХ BETA( 30), GAM MAC 30), QDOT(30. 30, 30) 
I ALP NSION TC 30, 30, 30), TSTARI( 30, 30, 30) 
DIMENSION TSTAR2( 30, 30. 30) 


REAL TINIT, TSTARI, TSTAR2, T, TEMP 

REAL CP, K1, RHO, R. RI, R2, TIME, DELT 

REAL LX, LY, LZ, DELX, DELY, DELZ 

REAL A, B, C, D. BETA, GAMMA 

REAL FLUX1I, FLUX2, FLUX3, FLUX4, FLUX5, FLUX6 

REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMB6 
REAL HI, H2. 13, H4, H5, H6 

REAL BIOI, ВТО“, BIO3, BIO4, BIOS, BIO6 

REAL SIGI. SIG2, SIG3, SIG4, SIGS, S1G6 

о 1 12 13.14, 115 

REAL Vi, V2, V3, V4. V5, V6 

REAL WI, W2, W3, W4, W5, W6 

ERE XI, X27, X3, X4, X5. X6 

REAL EPSI, EPS?2. EPS3, EPS:, EPSS, EPS6 

REAL QDOT, GEN, Y1 

REAL CTIME, TIMEO, TIME] 

REAL SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLARI, SOLAR?. SOLAR3, SOLAR4, SOLARS, SOLAR6 
REAL EARTHI, EARTH2, EARTH3, EARTH4, EARTHS, EARTH6 
REAL ALBED!, ALBED?. ALBED3, ALBED4, ALBEDS, ALBED6 
REAL ABSI, ABS2, ABS3, ABS4, АВ$5, ABS6 

MEI) SACI, SAC2, SAC3, SAC4, SACS, SAC6 

REAL SUN, ECLIPS, Q1, Q2, Q3, PERIOD, PI 

REAL P1, P2, P3, P4, P5, P6 


УЕ ЕКОЛ КОМ, ХВИ ИВА Т, б, О 
ИБРМЕКСОЛХТ, ЕКЕО, АХӘРАХЭШТАХОЗУАР, ХОСТ 


e. ЛО МИША ТЕ 15 ТАКТ,. TSTAR2, T, TEMP 
ИША 5IONCP PF, RHO, R. Ri. R2, TIME, DELT 
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COMMON LX. LY LZ DELS DELI DA 

COMMON I, J. K. M, APIF ТЕРІ ГАЗТ 

COMMON A, B, C, D, BETS, GCAM MARMO 

COMMON FLUXI, FLUX2, FLUX3, FLUX4, FLUX5, FLUX6 
COMMON TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMB6 
COMMON HI, H2, H3, H4, H5, H6 

COMMON BIO], BIO2, BIO3, BIO4, BIOS, BIO6 

СОММОХ 6ІСІ. 5168; 55613, 5164 5105 9156 

соммо хр то еее 

СОММОУ ХІ, V2, V3, V4, V5, V6 

COMMON W1, W2. W3, W4, W5, W6 

COMMON XL ADN. ALA 5 

COMMON COUNT. FREQ. AnS S ал в 
COMMON EPS1. EPS2. ЕР53, ЕР54, ЕР55, ЕР56 

COMMON QDOT, GEN. Y1 

COMMON TIMEO, TIMEI 

COMMON SOLAR. EARTH, RE, DIST, ALT, FE, FA, MU 
COMMON SOLARI, SOLAR2. SOLAR3, SOLAR4, SOLARS, SOLAR6 
COMMON EARTH}, EARTH?. EARTH3, EARTH4, EARTHS, EARTHÓ 
COMMON ALBED1, ALBED2. ALBED3, ALBED4, ALBEDS5, ALBED6 
COMMON АВ51. 4852. 4853. АВ54, АВ55, АВ56 

COMMON SACI SAC SACO ЕС 

COMMON SUN, ECLIPS, Q1, Q2, Q3, PERIOD. PI 

COMMON PI, P2. P3, P4, P5, P6 


IF( TIME .GE. (12*(SUN - ECLIPS))) .AND. 
& (TIME.LT. (13*SUN + 12*ECLIPS))) THEN 
GO TO 10 
ELSE IF((TIME .GE. (13*SUN - I2*ECLIPS)) .AND. 
& (TIME.LT. (13(SUN- ECLIPS))) THEN 


GO TO 20 


ELSE IF((IIME ЖЕ. (155(5 БЕ РЕСЕ РА 
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mE — = 


EEUU TEL (IJ SEN T 13° ECLIPS))) THEA 
GO TO TO 


ПИО IPO(TIME.GE. (14°SU`X + 13“ECLIPS)) .A`ND. 
c (TIME .LT. (14*(SUN+ ECLIPS)))) THEN 
GO TO 20 


ELSE IF((TIME .GE. (14*(SUN + ECLIPS))) .AND. 
& (TIME.LT. (15*SUN 4 14*ECLIPS))) THEN 
GO TO 10 


BEELIF(LIME.GE.(15*SUN + 14“ECLIPS)) .A`ND. 
ПОАТЈЕ ыо" МЕС) THEN 
GO TO 20 


Pier ite ОЕ, (1555 Х-ҒЕСІІР5))) AND. 
EI [PE LI (I6 SUN T ЕС) НЕХ 
GO IO 10 


ELSE IF((TIME .GE. (16*SUN + 15%ЕСІЛРЅ)) .АХР”. 
& (TIME .LT. (16 (SUN- ECLIPS))) THEN 
S STO 20 


ЕЕЕ МЕ GE. (167*(SUN + ECLIPS).AND. 
КИ II A IPL PT (17“SUN` + 16°ECLIPS))) THEN 
GO IO 10 


ELSE IF((TIME .GE. (17*SUN + 16*ECLIPS)) .AND. 
& (TIME .LT. (17*(SUN+ECLIPS)))) THEN 
GO TO 20 
ELSE 
Gono 10 
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ENDIF 
10 SOLAR = D3: 
WRITE(6,?) SOLAR 2 ,SOLAR 
GO TO 30 


20 SOLAR = 0 
WRITE(6,*) ‘SOLAR =’,SOLAR 
y 
30 RETURA 
END 


5УСВКОСТХЕЗ ХКАЛ!. ХКА 2, АМО ХКАТЗ СОМРЕТЕ СОВЕ I R 
ARRAY DIN THE X DIREC TI? 


ВОДЕ ИЕСИ 


SUBROUTINE XRAYI 

DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30). GAMMA( 30), QDOT( 30, 30, 30) 
DIMENSION T( 30, 30, 30), TSTARI( 30. 30, 30) 
DIMENSION TSTAR2 30, 30, 30) 


REAL TISSIT, TSTARI, ISTAR? T TEME 

REAL Cr AN RHO РОК АЦА Ова hi 

REAL LCA; LCY LZ. DELS DER DEEZ 

REAL A, B, C, D, BETA, GAMMA 

REAL FLUXI, FLUX? FLUX) FLUXIT М ЕРІ ЖАО 
REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMB5, TAMB6 
КРАВ ВА АРАТА H HS Hio 

REAL BIO1, BIOL, BIO3. BIO4, BIO5, BIO6 
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ЕСІ 5101, 5167, 5165. 5163, 5165, 5166 

ШЕ CI, U2. U3, U4, U5 

REAL V1, V2, V3, V4, V5, V6 

REAL W1, W2, W3, W4. W5, W6 

ШЕСІ ХІ, Х2, Х3, Х4, Х5, Х6 

REAL EPSI, ЕР52, ЕР53, ЕР54, ЕР55, ЕР56 

REAL QDOT, GEN, У] 

REAL CTIME, TIMEO, TIMEI 

REAL SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLARI, SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
REAL EARTHI, EARTH2, EARTH3, EARTH4, EARTH5, EARTH6 
REAL ALBED1I, ALBED2, ALBED3, ALBED4, ALBEDS. ALBED6 
REAL ABSI, ABS2, ABS3, ABS4, ABSS, АВ56 

REAL SACI. SAC2. SAC3, SAC4. SACS, 5АС6 

REAL SUN, ECLIPS, QI. Q2. Q3, PERIOD, PI 

REAL P1, P2, P3, P-1, P5, P6 


КО ОЕК ГЛ К.М. NN PLITESESIFPLI, EAST, 6, Q 
ПИО СОВА FREQ ANSI ANS522ANS5 МАГ, ХОСТ 


NO UUON IUNII, TSTARI, ISTAR2, T, TEMP 

more iON. CP, Kl. RHO, R, Ri, R2. TIME, DELT 

meet viON LX, LY. LZ, DELA., DELY, DELZ 

ПАО ТОЈ NK. M, N P. IF, L, IFPL, LAST 

mer isiON A,B,C, D, BETA. GAXIMA, G, Q 

ПАОЛА РЈ ОУ PLUS) FLUNS5SSFLUSmI I FLUX5, FLUX6 
COMMON TAMBI, TAMB2, TAMB3, TAMB4, TAMB5, TAMB6 
COMMON H], H2, H3, H4, H5, H6 

COMMON BIOI, BIO2. BIO3, BIO4, BIOS, BIO6 

COMMON SIG1, SIG2, S1G3, SIG4, SIG5, SIG6 

COMMON LI L2, U3, L4, L5 

COMMON VI Y2, Y3. V4, V5. V6 

COMMON WE W2 W3. W4, W5, W6 

ШОХ М ОМА ХМ? Хэ. ХӘ, Х5‚ хб 
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COMMON COUNT. FREQ 

COMMON EPSI. EPSV EPS) EPS4, PPF i 

COMMON ODOT, GESTY] 

СОМОНИ МЕО ТИЕ 

COMMON SOLAR, EARTH, RE, DIST, ALT, FE, FA, MU 

COMMON SOLARI, SOLAR2, SOLAR3, SOLARA, SOLARS, SOLAR6 
COMMON EARTHI, EARTH2, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBED1, ALBED2, ALBED3, ALBED4, ALBED5, ALBED6 
COMMON ABSI, ABS?, ABS3, ABS4, ABS5, ABS6 

COMMON SACI, SAGZ, SACs SACI Ger (O 

СОМ МОМОСУ О ООЗ ИЕТ 

СОММОХ РР ЕРГЕ ЕЕ 


ТЕ. РО. АР. Бок 


& 
& 
& 
& 
& 


D(I)- I(LJ,R) + ОЗЕ Ва 


ІЗТОЛІК АІ) ИІЛЕ = СТА е ES 
ҰЗҰ(ТАМВЗ-ТІГІК! (ТАМ ое Е 


EESE IFQ(I.GT.I) ANDER Py aN ТЕРГЕ 


DT) = IJR) LRF TE 


а ЛА КЕ то Де Qh A pi p I EE 
V3*(TAMB3 - T(1.J.K)) + W1 + W3 + 
XI*(T(LJ,K)**4 - TAMBI**4) 

X3*(T(1J,K)**4 - TAMB3**4) + YI*QDOT(LJ,K) 


ELSE IF((I.EQ.M) .AND. (J.EQ.1)) THEN 


Э) = ТКТ ГӘЗИТЕ ЕТЕР, 


U4*(T(1,J.K +1) - T(1.J,K)) + VI%TAMBI - № 
V3*(TAMB3 - T(LJ,K)) + V5*(TAMBS - T(I,J,K)) + 


МАТКУ АВ ии 
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x 
x 





& ОИ КОА ЕВРА ale ODOT(1.J-K) 


ШЕР IF(1.EO.1) AND, (3.GT.1) AND. U.LT.N))) THEN 

ШШ ey Cli del.) + 10J+1,k)- 2°101,J,K)) + 
U4*(T(1J,.K +1) - TJ,K)) + V3*(TAMB3 - T(I,J,K)) + 
V6*(TAMB6 - T(I,J,K)) + W3 + W6 + 
X3*(T(LJ,K)**4- TAMB3**4) + 
X6*(T(1J,K)**4 - TAMB6**4) + YI*QDOT(I,J,K) 


р р р Р 


Peo TIF((IL.GT.1) AND. (LLT. M) AND. (J.GT.1) AND. 
ЕТЕТ М))) THEN 

EO S= LILK) t LU1S(T(LJ-1.K) + T(I,I + 1,K)- 2°T(I.J,K)) + 
& ЛЬ =. = У АВ. КО тт 
& Е В) QODOTI.J.K) 


ШЕР ІРЕ(ІЕО./М).АХхХІ. ((1.6Т.1).АХР. (7.1.7.Х) »1 THEN 


Е ТОК) Е току а ТТК) отл) 
& UAJ*(T(LEK 8 1) - T(LI,K)) + V3°(TAMB3 - ӨӨ а 
& V5°(TAMBS - T(1J.K)) + W3 + WS5 + 
& N3*(T(LJ.K)?*4 - TAMB3**4) + 
& X5°(T(1.J.K)**4 - TAMB5**4) + YI®QDOT(I,J,K) 


ELSE IF((I.EQ.1) .AND. (J.EQ.N)) THEN 
D(l) = T(LLK) * U3*T(IJ-1.K) - T(LJ,R)) * 


& ОТИ) ТИЈ К) + У2=(ТАМВ2 - Та.) + 
& 39 (ТАМВЗ - ТП/.К)) + Ууб тАМВБ - ТК + С 
& Мо + W3 + W6 + X2*(T(I,J,K)**4 - TAMB2**4) + 
& X3*(T(1.J,K)**4 - TAMB3**4) + 

& X6*(T(1J.K)**4 - TAMB6**4) -- YI*QDOT(1,J, K) 


ELSE IF(((I.GT.1) .AND. (1.LT.M)) AND. (J.EQ.X)) THEN 
О ЮЕ ЗТ К) Так) 

& U4“ T(LI.K +1) - TUJ.K)) + V2*“(TAMB2 - T(I,J,K)) + 

& VSSCPANIBS - T(1.J.K)) + W2 МО 
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& Хо Ноу 
& X3*(T(LJ.K)**4 - TAMB3**4) - YI*QDOT(1J.K) 


ELSE 
Dil) = T(1,1,K) + U3*(T(1,J-1,K) - T(1,J,K)) + 


& Са “ТОК 1) - ТК) В У ИТТЕ ИСА 
& V3*(TAMB3 - T(1,J,K)) + V5*(TAMBS - T(1,J,K)) + 

& W2 + W3 +4 W5 Е ОПОЛ) A MB2**4) + 
& X3*(T(LJ,K)**4 - TAMB3**4) + 

& X5*(T(IJ,K)**4 - TAMB5**4) + YI*QDOT(I,J,K) 
ENDIF 

RETURN 

END 


SUBROUTINECSISAYS 


SUBROU DESI ATP 

DIMENSION A( 30). B( 30). CC 30) D ОЗО ЕМІН ЕСІ 
DIMENSION BETA( 30), GAMMA( 30 ODOT Оо) 
DIMENSION 1( 30:250 050] рол Оо) 
DIMENSION DSWAR2(30, 30989 


КЕЛТА TSTARI ISI AP s T P oj 

REAL CP. RI RHO R. RIP PTI Ni phi h 

REAL LX LY LZ DELS DI px pi pz 

REALA, B,C, D, BETA GCS MMA 

КЕЛІ, ЕБЕХІСКІЛХ РЕГХФ ЕБ ЕРЕ ЕЕ О 
REAL TAMBI, ТАМВ2, ТАМВ3, ТАМВ4, ТАМВ5, ТАМВ6 
REAL HI H2, H3 H4 H3 HC 

REAL BIOI, BIO?, BIO3, BIO4, BIOS, BIO6 

REAL SIG1, SIG2 SIG3 SIG ЛЕС 


БЕСІТІ (215 14,125 

REAL V1, V2, V3. V4, V5. V6 

МЕИ ОЈ. АКО. W3, W4. W5, W6 

FISSE X1, X2, X3, X4, N5. X6 

REAL EPS1, EPS2, EPS3, ЕР54, ЕР55, ЕР56 

REAL QDOT, GEN, YI 

REAL CTIME, TIMEO, TIMEI 

REAL SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLARI, SOLAR?, SOLAR3, SOLAR4, SOLARS, SOLAR6 
REAL EARTH], EARTH?, EARTH3, EARTH4, EARTHS, EARTH6 
REAL ALBEDI, ALBED2. ALBED3, ALBED4, ALBED$, ALBED6 
REAL ABSI. ABS2, ABS3. ABS. ABS5, ABS 

ПЕ 5401. SAC?. SAC3. SACA, SACS, SACÓ 

REAL SUN. ECLIPS. Q1. Q2, Q3, PERIOD. PI 

MEAL Pl. P2, P3, P4. PS. P6 


A CIEI J K.N.SN.P. IF. L.IFP1 LAST, G, Q 
C IP COL ОЕ РКЕО А А\52. А\55, УАГ, ХОСТ 


BEND UONTINIT,TSTARIH STARS, T, TEMP 

О СВ. БВ КТО, В, КВ. TIME. DELT 
ОХ Гл. | (и, ПЕРА DEL Y. DELZ 

ОХ ВОВА, РАЕ, Е ТЕРЕЕАЗТ 

ОХ A. B. C, D. BETA, GAMMA, G, Q 

I IONS TITI ad, PLUN2, FLUNS, PLUASS PEL AS, FLUXG 
COMMON TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMB6 
GO VLVION HI, H2, H5, H4, H5, H6 

COMMON BIO], BIO2, BIO3, BIO4, BIOS, BIO6 

COMMON SIG1, SIG2, 5163, 5164,5165, 5166 

SONIVION Ul, U2, 030425 

C O NI MOSS Vl V2, V3, V4eV5, V6 

COMMON W1, W2, W3, W4, W5, W6 

COMMON XI, 2, X3, X4, X5, X6 

COT VION COUNT. FREQ 


СОММОС-ЕРБІГЕС52 ЕР ЕР5ЗЗЕЕЗЗ ЕНСЕ 

CON MO N QDOT, GEN. Y1 

СОММОХ ПМЕО,ТІМЕІ 

COMMON SOLAR, EARTH, RE, DIST, «Ее Е КЕ 

COMMON SOLARI, SOLAR2, SOLAR3, SOLARA, SOLARS, SOLAR6 
COMMON EARTHI, EARTH2, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI, ALBED2, ALBED3, ALBED4, ALBEDS, ALBED6 
COMMON ABSI, ABS2, ABS3, АВ54, АВ55, АВ56 

COMMON SACI, SAC2, SAC3, SACA, SACS, SACÓ 

COMMON SUN ЕВ 0 ОВ ом 

СОММОСУРЕЕЛ РЗ БӘ 


ІЕ( ЕО. Ар ТЕО Ех 
DT) = ТАЈ,К) 95 152 ШЕ Е О: 


& Зета К) О) r VI Е 
& V6*(TAMB6 - T(IJ,K)) + W1 + W6 + 

& X1°(T(J,K)""4 - TAMBI**4) + 

& X6*(T(IJ,K)**4 - TAMB6**4) + Y1*QDOT(I,J.K) 


ELSE IF(((1.GT.1) -AND. (1.LT.M)) .AND. (J.EQ.1)) THEN 
>) + ТК) + U2%(T(1.J.K-1) + T(1J,K+ 1) - 22T(1J.K)) + 

& U3*(T(LI- LK) - TL,K)) -- VITAMBI - TÁK) + 

& WI + XIXT(LJ.K)^*4- TAMBI**4) - YI*QDOT(LJ,K) 


ELSE IF(LEQ.M) .AND. (J.EQ.1)) THEN 
О) = так) + U2*(KIJ,K-D + a = 
U3*(T(IJ- LK) - T(LJ,K)) * VI*(TAMBI - T(LLK)) + 
V5*(TAMBS - T(1.J,K)) + WI + W5 + 
X1*#(T(1,J,K)**4 - TAMBI#*4) + 
X5*(T(1,J,K)**4 - TAMB5**4) + Y1*QDOT(I,J,K) 


рр р р 


ELSE ІЕ(П.БОЛ) АХ, (БӨЛІ) АЗЫ ЕЕ 
рЫ(І)- ТОЛҚУ КІЛТІНЕ ШЕГЕ О j S 
е U2*(T(LJ.K-1) * T(IJ,K £1) - 2*T(LJ,K)) + 
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6: V6"(TAMBG6 - T(I.J,K)) + W6 + 
& X6*(T(1J.,K)**4 - TAMB6**4) - YI1*QDOT(IJ,K) 


ELSE IF(((I.GT.1) .AND. (I.LT.M)) .AND. ((J.GT.1) .AND. 
& (J.LT.N))) THEN 
Ра) LJK) + UI*(T(,-LEK) * T(J - LK) - 2*T(0],K)) 
& ОТОК) TO,K 0) - 22T(J K)) + 
& Y1*QDOT(I.J,K) 


ELSE IF((I.EQ.M) -AND. ((J.GT.1) .AND. (J.LT.N))) THEN 
D(I) - T(1.J,K) + U1*(10,J-1,K) + T(LJ+1,K) - 2*T(1,J.K)) + 


& Ë “(I РАИ А РУАН АЕН РРА ЈА ју) + 
& НЕАП NIB = EIN NS 
& Peden) a= PA IDS 23) T Y12QODOT(LJ,K) 


ELSE IF((I.EQ.1) AND. (J.EQ.N)) THEN 

DI) = TUJ,K) + U2*(T(J.K-1) + TULJ.K+1)- 2*T(LJ.K)) + 
ШООК КОЛГЕ ЛЕ GA Ip (IJ K) + 
V62(TAMB6 - T(1.J,K)) + W2 + W6 + 
X2*(T(LJ.K)**4 - ТАМВ2%%) + 


р р р р 


ELSE IF(((1.GT.1) .AND. (LLT.M)) .AND. (J.EQ.N)) THEN 
D(I) 2» T(LJ.K) 4 U2*(T(LJ.K-1) * T(OLLK o 1) - 2*T(LL K)) + 


& U3*(T(LJ-lL.K) - T(LJ.K)) *. V2(TAMB2 - T(OLK)) + 
& мо + N2*(T(IJ,K)**4 - TAMB2**4) + Y1*QDOT(I,J,K) 
ELSE 


We DER не де Е ТООКЕ - 2500 P R)) 9 
U3*(T(1,J-1,K) - T(1,J,K)) + V2*(TAMB2 - T(I,J,K)) + 
V59(TAMBS - T(1,J,K)) + W2 + W5 + 
X2°(T(1,J,K)**4 - TAMB2*#4) + 
X5*(T(J.K)**4 - TAMBS**4) + Y1*QDOT(I,J,K) 


po P» P Ро 
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ENDIF 


КЕТЕК 
ESD 


SUBROUTINE XRAY; 


SUBROUTINE Sans 

DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAM MA( 30), QDOT( 30, 30, 30) 
DIMENSION T( 30. 30, 30), TSTARI( 30, 30, 30) 
DIMENSION TSTAR2 30, 30, 30) 


REAL TINIT, TSTARY, АЕ ТЕ S 

REAL CP, Rig RHO? Re Rie Reve ie ее 

REAL EX, L3 ЕХ DEL БИЕ МЕТІ 

REAL A.B.C, D BETA GAIA 

REAL FLUXI, FLUA? FLUNSSA FECIT EC ee го 

REAL TAMB1, TAMB?. TAMB3. TAMB4. TAMB5, TAMB6 
КЕАТ НЫ В. В мо 

REAL BIO], BIO2, BIOS, BIO4, BIOS, BIO6 

REAL SIG], SIG2ZISIG3s, S| GaesiG oes Ge 

КЕЗЕ СЕТУЕТ 

REAL. \V зен лека 

REAL W1, W2, W3, W4, WS. W6 

ВАЈАТ АРА за А О 

КЕДЛІСЕББІСЕР52 ЕР ХЕБ ERSS PESO 

REAL ODOT, GES, YI 

КЕ ОЛИ РО ВУ РО УП и 

REAL SOLAR, EARTH, ALBCO, RE DIST ALTHER TAE 
REAL SOLAR1, SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
КЕЛІ, ЕАКТНІ, ЕАКТН2, EARTH3, EARTH4, EARTHS, EARTH6 
REAL ALBEDLI ALBED2, ALBEDS, ATBEDASSASPDEDS SAT DISIDO 
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REAL ABSI, ABS2, ABS3, АВ54, АВ55, АВ56 

МЕИ SACI. SACI- SACI, SAC4, SACS, SAC6 
REAL SUN. ECLIPS, Q1, Q?, Q3, PERIOD, PI 
REAL PI, P2, P3, P4, P5, P6 


СЕВ ИР Г. [ЕР]. LAST, G, Q 
BSIEGER COUNT, FREQ. ANSI, ANS2, ANS3, VAL, NOUT 


ШОУ МОХ ТГМЇТ, Т5ТАК1, Т5ТАК2, Т, ТЕМР 

БИОС СР, КІ, КНО, R, RI, R2, TIME, DELT 

O I NMON LN. LY. LZ, DELX. DELY. DELZ 

ша MON i, J, kK, MN. PIF. L, IFPI. LAST 

Ps I [ON AD COD БЕТА GAMMA, G. Q 

ПЕРА Е пор Вера ЕР ХО O FLUN3S Евка ЕО Хо, ЕС ХО 
COMMON TAMBI. TAMB?. TAMB3. TAMB-4. ТАМВ5, ТАМВо 
ШОМ МОМ HI. F2, H3. H4. HS. H6 

COMMON BIOI, BIO2, BIO3, BIOJ. BIOS, BIO6 

ШОЛ МОХ 51651. 51652, SIG3. 516-4. 51655. 51656 

ШОММОХМ 11. 62. 23, 24205 

E MON VI, V2. V5. V4, V5. VO 

ШОО МОМ М. М2. М5. 4 \5, Мо 

ШОМ ОХ NI. AZ, AS, Ad. AS, AG 

SOXIMON COUNT, FREQ 

MONON EPSI EPS2.CPS3, EPS4. EPS5., EPSG 

SONINION ODOL GEN. Y1 

СОММОХ ПІМЕО, TIMEI 

COMMON SOLAR, EARTH, RE, DIST, ALT, FE, FA, MU 
COMMON SOLARI, SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
СОММОХ ЕАКТНІ, EARTH2, EARTH3, EARTH4, ЕАКТН5, ЕАКТН6 
COMMON ALBEDI, ALBED?2, ALBED3, ALBED4, ALBED5, ALBED6 
COMMON ABSI, ABS2, ABS3, ABS4, ABS5, ABS6 

COMMON SACI, SAC?, SAC3-SAC4, SAC5, SAC6 

CONOS SUS ECLIPS. QL O2, Q3, PERIOD, PI 
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COMMON PIP? Е 9 


[F(LEO.1) AND. (.EO.1)) THEN 
РО сс ТОК A а T IJ PK ЕТИ 


& LA4*(T(1,J.K-1) - T(LJ,R)) - VI*(TAMBI - T(LJ,K)) + 

& V4*(TAMB4 - TU,J,K)} + V6*(TAMB6 - T(I,J,K)) + W1 + 
& W4 + W6 + X1*(T(U,J,K)*"4 - TAMBI**4) + 

& AJA*(T(L,J,K)**4- TAMB4**4) + 

& Хе ЦК) 74 - ТАМВО 4) а РО РОМА) 


ЕТЗЕ ТЕ ОТ1) АМОС (ЕЕ ТЗ ода ТЕЛЕ 
рају = ЕЛЬ) = О ОИ ЕК а EE 


& САТ К р TANET ен 
& VA*(TAMBA - T(LJ,K)) + WI + W4 + 

& ХАТТЫ 

& XA*(T(LK)**4 - TAMB4*#4) + Y1"QDOT(IJ,K) 


ELSE IF((I.EQ.M) .AND. (J.EQ.1)) THEN 
DU) = TUK) + (T py ea ee 


& САТИ К ТООК) ЕТА ee Tee 
& VITAMBI MIIR) + V5*(TAMBS - T(1.J,K)) + W1 + 
& ХА + WS + X1*(TULILK)**4 - TAMBI*#4) + 

& XA*(TOLJ.R)*74- TAMBA**4) 

& XS*(T(LJ.K)*^4 - TAMB5223) + YI*QDOT(LJ,K) 


ELSE IFGLEQ I) .AND (J.GI 1) AND (JDI Е 
D(D = ТК) + UI*CT-LEK) *- T(J4 LE) - 2?T(,A K)) + 
U4*(T(1,J,K-1) - TUJ,K)) + V4*(TAMB4 - T(IJ.K)) + 
V6*(TA MB6 - T(1,J,K)) + W4 + W6 + 
X4*(T(1,J,K)**4 - TAMB4**4) + 
X6*(T(IJ,K)**4 - TAMB6**4) + Y1*QDOT(I,J,K) 


Po фо Po Р 


ELSE 1F(((1.GT.1) AND. (LLT.M)) .AND- (UJ.GT.D AND. 
SNO DIES SEDENS 
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EIL MIU coi Ns ILEHR)-2*I(L Jk) + 


& 
& 


UA*(T(LJ.K-1) - T(IL.L,R)) -- VA*(TAMBA - T(LJ,K)) 9 


BESETE(QEEQO. MV) ACND-(U.GT.D) -.AEND. CLELT.N)) THEN 
Em —ULNUTUOUPCOJLEK&)- c I(LJT LR) -2*T(IJ,K)) + 


& 
& 
& 
& 


UA*(T(LJ.K-1) - T(LJ,K)) - VÁ(TAMBA - T(LJ,K)) * 
V5*(TAMBS - T(1,J,K)) + W4 + W5 + 
XA*(T(IJ.K)**4 - TAMBa**4) + 

Х5“Т(1,1.К)4%4- ТАМВ5%%4) -- YI*QDOT(L,K) 


ELSE IF((I.EQ.1) .AND. (J.EQ.N)) THEN 
may = IJK) + U3ST(J-EK) - TCR) + 


Л ЕК Е О ИВА T LJ E) + 
VS (TANIB3 - TO(LJ.E)) + V6*(TANMNIB6- T(LJ K) + W2 + 
И о ОК 04 <- TAMB2**4). T 


о БЕ: YI ODOT(I,J K) 


i | IP(((LGTLI) AND. (ILT.ND) . AND. J.EQ.N) THEN 
may I(LJ.K)-* LS*(TI(IJ-I. KR) ТК) + 


ELSE 


Aik lel Ree TANE E TIEN 
V4*(TAMB4 - TULJ.K)) + W2 + W4 + 
XO*(T(L,K)**4- TAMB2**4) + 

XJ*(T(LJ,K)**4- TAMBA*94) + YI*QDOT(LA,K) 


ID SOBRE T USTCIQLI-Li&)- I(LJ; R0) + 


рр р р р 


UA*(T(LJ,K-1) - T(LJ,K)) * V2S(TAMB2 - ТА ЈАК)) + 
V4*(TAMB4 - T(I,J,K)) + V5*(TAMBS - T(1,J,K)) + W2 + 
W4 + W5 + X2*(T(I,J,K)**4 - TAMB2**4) + 
X4*(T(1J,K)**4 - TAMB4**4) + 

ОШИО МВО А) p YISODOT0 JE) 


liil 


ES DIEF 


БЕТ М 
ESD 


ои о ие Ук е Хе И е И sh 32 За ти 212 22 ПТ Ни и За Ха 55:5 5 ус Ж: 55:5 SE У У F> 3ç 55 X; X; X 392 Z: X; Yz Z 352% E 9:25: Ж 


~ аб 3% ate ФУ «и ад ай, $34 2594-5 те уа Уа wie ate ole ае «е то ада уби зи cle ole ob 
тт зо до го Ser Bye sgt set ge oye е Ме ѓа е оо ге ок е 


% 
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SUBROUTINES YANKI YANK2, AND NANK5 COMPUE TEGOEFFIC I Sl 
ARRAY DIN TFE DIRECTION 


УРА О И и и hi 


SUBROUTINE М 

DINME`NSION A( 350). ВО ОИ о ВО 
DIMENSION BETA( 30) GA MATA НА А ТА РОБНЕ) 
DIMENSION T( ЭР РТ ЕГЕ ЕТ 
DIMESSION TSTAR ПТ! 


REAL TINIT. TSTARI, TSTAR2. T. TEMP 

REAL CP. K1. RHO. R, R1. R2. TIME. DELT 

REAL EN. LY, LZ. DELN, DELY DELZ 

REAL A, B. C. D, BETA. GAMMA 

РЕА ЕСУІ РИО РОХ ВИВЕСТИ С 
REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMBG6 
REAL Hl, H2, H3, H4, H5, H6 

REAL BIO}, BIO2, BIO3, BIO4, BIOS, BIO6 

REAL SIGI, SIG2, SIG3, 5104, 5105, 5106 

REAL CI PONE: 

REAL о аа 

REAL W1, W2, W3, W4. W5, W6 

ВЕРА оо 

REAL EPSI. EPS2, EPS3, ЕР54, ЕР55, EPS6 


— —  — —— —a s s ƏswcuS1yH!.x Fn — 


REAL QDOT. GEN, Y1 

REAL CTIME, TIMEO, TIME] 

REAL SOLAR, EARTH. ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLAR}, SOLAR?. SOLAR3, SOLAR4. SOLARS, SOLAR6 
REAL EARTHI, EARTH2, EARTH3, EARTH4, EARTHS, EARTH6 
REAL ALBEDI, ALBED2, ALBED3, ALBED4, ALBEDS, ALBED6 
REAL ABSI, ABS2, АВ53, АВ54, ABS5, ABS 

REAL SACI, SAC2, SAC3, SAC4, SACS, SAC6 

REAL SUN, ECLIPS, Q1, Q2, Q3, PERIOD. PI 

REAL P1, P2, P3, P4, P5, P6 


ПИ CGO TLJ. K, M; A, P, IF; L, IFPI- LAST, G, Q 
СЕК COLNT FREQ ANSI ANS2..ANS3, VAL, NOLT 


ШОО МОМ ТГМЇТЇ, УЛАК TSTAR?. T. TEMP 

Шик МО СР. КТ, RHO, R. RI, R2, TIME. DELT 

ИШЕ ON LN TY EZ DELX DELY. DELZ 

ша ON], J, hh, E N. P IF, L IFPI, LAST 

О ^. В. С, Б. ВЕТА, САММА, 6. Q 

О РЕ FEU N22 PLUXS, FLUMS. FLUXS, FLUXG 
БӘСІМО TANIBI, TAMB2, TAMBS. TAMPBA, TAMBS, TAMB6 
МАЧЕТО)“ ТТ. НО. НИ. HS5. НӨ 

EO DMOSN BIO!, BIO2, BIOS, BIOJ, BIOS, BIOG 

er ON SIG, SIG2, SIGS. SIGGRSIGS4S1GG 

ШОМ Ох 1,12, [ 5, [м5 

ПАЛА IOS VL v2, V3, V4, VS, V6 

COMMON W1, W2, W3, W4, W5, W6 

ЕО МО XI, N2, X5, X4, X5, X6 

COSIMOASA COLAT. FREQ 

SS TVION EPSI, EPS2, EPS3, EPSA, EPSS, EPSO 

COM MOS ODOT, GEN YI 

COMMON TIMEO, TIME] 

CON MON SOPAR EARTH RE, DIST, ALT, FEFA, MLE 
COMMON SOLARI, SOLAR2, SOLAR3, SOLARA4, SOLARS, SOLAR6 
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COMMON EARTHI, EARTH2, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI. ALBED?2, ALBED3, ALBED4, ALBEDS5, ALBED6 
COMMON ABSI. АВ$>, АВ$3, АВ54, АВ55, ABS6 

COMMON SACI. SACI RAC ЗАСТ аас 

COMMON SUN, ECLIPS, QI, Q2, Q3, PERIOD, PI 

COMMON PI, P2, P3, P4, P3, P6 


IF((J.EQ.1) AND. (1.EQ.1)) THEN 


DJ) = ТОК) + СВЕН Е 


U5*(TSTARI(I+1,J,K) - TSTAR1(1,J,K)) + 
VITAMBI- TOLKI + V3*(TAMB3 - T(I,J.K)) + 
V6"(TAMBG - T(LJ.K)) + WI + W3 + W6 + 
Х1*(Т(1..К)**4 - ТАМВ1**4) + 


ELSE IF(OU.GTED) 2A NDIGDISISBSS ОИЕ чх 


& 
& 
& 
& 


D(J) = ОК) + 14 ЕО КВ I 0 Т 


Lo (TSISRI I IJ r T TI I | 
УЗТАМВЗ - ТАЈК)) + УбЧТАМВе - ТАЈ КУ + 


Хо (ПК) 74 - TASB 6" 23) 9 УЛЕС БО чи 


ELSE IF((J.EQ.N) AND. (LEQ.1)) THEN 


Po р PP Po Pp 


DJ) = ИЛЬ) + ЕКО + Ш. ЦЕ 


US*TSTARI(I+1J.K) - ТТАВЛИ Ј,К)) + 
V2*(TAMB2 - T(LJ,K)) - V3*(TAMB3 - T(LJ,K)) * 
V6*(TAMB6 - T(I,J,K)) + W2 + W3 + W6 + 
X2*(T(1,J,K)**4 - TAMB2**4) + 


ELSE IF((J.EQ.1) AND. ((I.GT.1) .AND. (LLT.M))) THEN 


D(J) ТИЈ Ку САНИ ЈУ CSS OI ЕР 
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& Fo on Ae ee ers DARK) = 2°7STARI(1,J,K)) 


& ee elec the WN yet evo (PAMB3- TU,J,K)) + WI + 
& W3 * XI*(T(LJ.K)**4 - TAMBI**4) 
& DENM ИУ ОИ Орои, К) 


ШЕЕ ТЕСЕ) АЮ. О Г.А) .ANND. (L.GT.I) .AND. 
Е. Т.М))) ТНЕХ 
DD S= ШЕК) САТИЛ КЕТ) ОК) + 


& FO*(TSTARK(I-1,J.K) + TSTARI(I+1,J,K) - 2*TSTARI(L,J,K)) 
& ИЕЗИ ТК) + W3 
& X3*(T(LJ.K)**4 - TAMB3**4) + YI*QDOT(LJ.K) 


ЖЕБЕ ТЕС (ОЕ АЛЬ. ТЕТ. М))) THEN 
р ИЛ - T(1J.K)) + 


& FO*(TSTARI(I-1..K) ^ TSTARI(I LJE.K) - 2*TSTARI(I,J.K)) 
& В К» + М + 
& УЗ + Х2*(Т(1.1.К)* 24 - ТАМВ2**4) + 

& X3*(T(LJ.K)**4 - TAMB3**4) + Y1*QDOT(IJ,K) 


ВЕ ЕТЕ( ЕО) АХР. (ШЕО. МП) THEN 
ШО К j I OCT Ida 1 CPO DE) + 


& ОКИ Е ТУЛЕ ВИИ Е) t 

< MISSIS ТАМ ТЕЛЬ + 
& ИВ ПЕР» NICE NS NT 

SS АГ (АЕК) TANEI cr 

& ть)” А ТА МЕР А) + 

& AS*(TUJ,K)**4 - TAMBS**4) + YI*QDOT(I,J,K) 


Been Giles Neel aN) ekND. (1.EO.M)) THEN 
DJ) » T(LJ,K) + U4*(TUJ,K +1) - TUJ.K)) + 


& U5*(TSTARI(I-1,J,K) - TSTARI(I,J,K)) + 
& V3*(TAMB3 - T(I,J.K)) + V5%(TAMBS - T(I.J,K)) + W3 + 
с W5 + X34(T(LJ.K)**#4 - TAMB3**4) + 
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& X5*(T(1J.K)**4 - TAMBS**4) + Y1*QDOT(IJ.K) 


ELSE 
ро) - T(LI.K) + U4*(T(LJ.K +1) - T(LJ,K)) + 
& L5*(TSTAR ICAL IK) = ISmee Gn ere eer 
& V2*(TAMB2 - T(1.J,K)) + V3*(TAMB3 - T(I,J,K)) + 
& VS*(TAMBS5 - T(L,K)) + W2 + W3 + W5 + 
& X2*(T(IJ,K)**4 - TAMB2**4) + 
& X3G*(T(IJ.K)**4 - TAMB3**4) + 
& XS*(T(IJ,K)**4- TAMBS**4) + Y1*QDOT(I,J,K) 
ENDIF 
RETURN 
END 


ОЗО ОСИ УЕ 


SUBROUTINE YANK? 

DIMENSION A( 30), BC 30), C( 30). D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAMMA( 30), QDOT( 30, 30, 30) 
DIMENSION Tt 30, 30, 30) TSTARIN ТТ 
DIMENSION TSTARX 30, 30, 30) 


КЕЗЕ ИО ТУЕ ИИСИ И У 

REAL CP, RI, RIO R RI PO HINDI I | 

REAL LA Li LA DEESIDE DEE 

REALA, B,C. D BETA. GAMMA 

REAL FLUX), FLUX ЕЕЗ ци ла Ри а вна 
REAL TAMBI, TAMB2, ТАМВ5, ТАМВ4, ТАМВ5, ТАМВ6 
КРАЈИНА А HS B3 EIS ERG 

REAL BIO1. BIO?2, BIO3. BIO4, BIOS, BIOG 

КЕЛДІ; 8ІСІ, 5162, 81027516 У ас бо 


| 





REAL UT, U2, U3, U4, US 

REAL V1, V2. V3, V4. V5, V6 

REAL WI, W2, W3, W4, W5, W6 

ВЕК X1 XO, X3, X4, X5, X6 

REAL EPSI, EPS2, EPS3, EPS4, ЕР55, ЕР56 

REAL QDOT, GEN, Y1 

REAL CTIME, TIMEO, TIME] 

REAL SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLARI, SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
REAL EARTHI, EARTH?, EARTH3, EARTH4, EARTHS, EARTH6 
REAL ALBEDI. ALBED2, ALBED3. ALBED4, ALBEDS5, ALBED6 
REAL ABSI, ABS2, ABS3. ABSJ. ABS, ABS 

REAL SACI, SAC2. SAC3, SAC4, SACS, SAC6 

REAL SUN. ECLIPS, Q1. Q2, Q3. PERIOD, PI 

REAL P1. P2, P3. Pl, P5, P6 


EN CER LE K. NM. N, P. IF, E; IFPILEAST, 0, Q 
СЕК СОС ЕКЕО АМУ, А52 А5), VAL, NOUT 


COMMON TINIT, TSTARI, TSTAR?, T, TEMP 

COMMON CP, K1, RHO. В. RI. R?, TIME, DELT 
ISONDMON EX, LY. LZ, DEEX"DEEY; DELZ 

ШООК Ох | J. K. N.N. P. IE. L. IFP1L LAST 

COMMON A. B. C, D. BETA, GAMMA, G, О 

Moron FILUX1, FLUN?, FLUNG, FLUNG, ELUXS, FLUX6 
COMMON TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMB6 
COMMON HI, H2, H3, H4, H5, H6 

COMMON BIOI, BIO2, BIO3, BIO4, BIOS, BIO6 

COMMON SIGI, SIG2, SIG3, SIG4, SIG5, SIG6 

IBDNIMON UL E2, СЗ, U4, US 

COMMON VI, V2, V3, V4, V5, V6 

COMMON WI, W2, W3, Wd, W5, W6 

(ОК ОХ е X2, 3, 8, 885. Хо 

COMMON COUNT, FREQ 


ен 


COMMON EPSI, EPS2, ЕР53. ЕР54, ЕР55, ЕР56 

COMMON ODOi Gia 

COMMON TIMEO, TIME! 

COMMON SOLAR, EARTH. RE, DIST, ALT, FE, FA, MU 
COMMON SOLARI, SOLAR?, SOLAR3, SOLAR4, SOLARS, SOLAR6 
COMMON EARTHI, EARTH?, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI, ALBED?, ALBED3, ALBED4, ALBED5, ALBED6 
COMMON ABSI, ABS2, АВ55, АВ54, АВ55, АВ56 

COMMON SACI, SAC2, SAC3, SACA, SACS, SACÓ 

COMMON SUN, ECLIPS, Q1, Q2, Q3, PERIOD, PI 

COMMON PI, P2, P3, P4, P5, P6 


ЕЕ АО ПЕ ТТ а 
D(J) = KLIK) +A CANLI E РАИ АРА 0) У ИЕ ПИГ х 


& US*(TSTARI(I+1J.K) - TSTARI(LAJ,K)) 

& - VIP(TAMBI - T(LL,K) t Vé*(TAMBG - T(J,K) + WI 9 
& W6 + XI*(T(LJ.K)**4 - TAMBI?*4) + 

& X6*(T(LJ.K)*^4 - TAMB6**4) + Y1*QDOT(IJ.K) 


ELSE 1F(((.GT.1) AND. U0. Nj) ANDMORE Oi ERES 
DJ) = T(LLK) * UAIR D - TOEK * 1) - 2?T(01J.K)) 9 


& US*(TSTARIG  1J,K) - TSTARI(I.K)) 
& + V6°*(TAMBG - T(1J.K)) + W6 + 
& Х6 ТИЈ К) А - TAMB6#*4) + Y1*QDOT(L.J,K) 


ELSE IF(J.EQ.N) .AND. (LEQ.D) THEN 
О(р е ЦК) * U2^(T(LLK-D) * T(LAK- D) - 2?T(LJK)) + 


& US*(TSTARI(I * 1J,K) - TSTARI(LJ,K)) 

& + V2*(TAMB2 - T(LJ.K)) + V6%(TAMB6 - T(I,J,K)) + W2 + 
& W6 + X2*(T(1,J,K)**4 - TAMB2**4) + 

& X6*(T(IJ.K)?*4 - TAMBG**4) 4- YI*QDOT(LJ,K) 


ELSE IF((J.EQ.1) .AND. ((I.GT.1) .AND. (LLT.M)) THEN 
D(J) = T(Q.LK) * U2NT(LJ.K-1) * T(OLALEK + 1) - 2*T(LJ,K)) 9 





& В в РЕ) 221STARKIJVINKIÇ I 
& РАВ (Е 
& DM cJ RS - DNNIBIS d) 5- YPEODOT(IJ,K) 


BISETE(G.GT.D .AND- (.LT.N)) АХО. (@.СТ.1) .АХО. 
& (LLT.M)) THEN 
ВИ К) (ТК) + T(LJK+ I) - 2*T(1J,K)) + 
& FO*(TSTARI(I-1,J,K) + TSTARI(I+1,J,K) -2*TSTARI(I,J,K)) 
& + Y1°QDOT(I,J,K) 


mae (()-EO.N) AND. ((1.GT.1) AND. (.LT.M))) THEN 
EX ЛОК У АВА ОЈ 1) = ТАЈА РО А ТА ЈАКУ + 


& ОИЕ ИТЕ Е ИБИ ЕО ЈО 2 TSTARNILIK)) 
& ee ANI = UL ISK) ае F 
& ПАЛА =- TAMB?" IP VI ODOT TK) 


ELSE IF((J.EQ.1) .AND. (LEQ.M)) THEN 
D) 9 T(LLK)- U2*«(T(LLK-1) + TULJ,K +1) - 2*TQ,J,K)) + 


& ü (IST APHIIIE) TSTARI INDO 

& БІРІ ТІК) AS TA MBS] ЛЕ» + Wd + 
& WS + NI“(T(LJK)**4 - TAMBI“*4) + 

& XS*(T(LJ.K)**4 - TAMBS**4) + YI*QDOT(IJ,K) 


ELSE IF((J.GT.1) .AND. (J.LT.N)) . AND. (LEQ.M)) THEN 
ШО = ТК) + U2*(T(LEK-0)) 9 TOJ,K 40) - ТОК) + 


& ВЕТ К, ISIARI IRE) 
& + V5*(TAMBS - T(I,J,K)) + W5 + 
& X5*(T(1J,K)**4 - TAMB5**4) + Y1*QDOT(I,J,K) 
ELSE 
Ра) = T(LLK)- U2*T(IJ,K-D) + NLK D 2*T(1,J,K)) + 
& ОЕШ ТК) ТВИТ Же) 
& + V2*(TAMB2? - T(L.J,K)) + V5%(TAMBS5 - T(IJ,K)) + W2 + 
& \'5 + X2*(T(LI,K)**4 - TAMB2**4) + 
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& ХОСЕ УНА TAN BS T Yr ODOTTI] 


ENDIF 


КЕТЕ 
ESD 


SCBROULUTISE TAn 


SUBROUTINE YAN 

DIMENSION AC 30), B( 50. C ВОИ D ( oO RE [I O, 
DIMENSION BET 4( 30). GAMMA( 30), QDOT( 30, 30, 30) 
DIMENSION T(. 305 30. ЗО ви OEC UD 
DIMENSION ISTAR2( 30, 30. 50) 


REAL FINIT, TSTARI TSTAR? TRENIE 

REAL CP. KR RHO R. RI R ISPP DI p ñ 

KEAL EN У LZ DELN DEROA 

КЕСА АИВ, С ОВОГА Сая: 

REAL ЕОСХІ РЕСХОСЕСЕХЕТ ЕВЕ ЕЕ ШЕКЕ: 
REAL TAMBI, ТАМВ2, ТАМВ3, ТАМВ4, ТАМВ5, ТАМВ6 
БЦП ИТТЕ ИИ ЕГО 

REAL BIOI BIO? BIOS IBIOTBIO ОТЕ 

REAL SIGI,SIOS SIGS. SIGSMSIGONSIGOÓ 

REAL UITE ГЕ 

REAL Vie sawn aN SN 

REAL W1, W2, W3, W4, W5, W6 

REAL AIFA2 AD AAI AD О 

REAL EPS EPS? TEPS ЕРОЧРБРБЕНЕС 

REAL ODOT GEN iI 

REAL CTIME, TIMEO: ТРЕ 

REAL SOLAR, EARTH: ALBCO, RE, Dis ЕВ ВЕ o I 
REAL SOLAR],  SOLAR2 SOPAPR3 ЗОВ ORSR 
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ШЕЛГЕАКТИЕП, ЕАКТН?” ЕАКТИ5, БАКТН4, ЕАҚТИ5, ЕАКІН6 
DESI ALBEDI. ALBED?2, ALBED3, ALBED4, ALBEDS, ALBED6 
REAL ABSI, ABS2, ABS3. ABS4, ABS5, ABS6 

ШЕЛІ ЗАС 5АС?2, 5АС3, SACA, SACS, SACO 

ШЕСІ SUN, ECLIPS, QI]. QO?, Q3, PERIOD, РІ 

ESL P1, P2, P3, P4, P35, P6 


A EGERI J, K: M,N, P, IF, L, IFP1, LAST, G, Q 
E EGER COLAT, FREQ, ANS51, ANS52, ANS53, VAL, NOUT 


FOM MOAN TINIT;, TSTAR1I, TSTAR2, T, TEMP 

BO INON CP. RE RHO R, R1, R2 TIME, DELT 

ше O [SY LZ DELXNX. DELY DEL Z 

ЕО L J. F M, A.P, IF L. IFPI LAST 

movi NON A. B,C, D. BETA. GAMMIA. G, QO 

ИШЛЕ O [I I Sl АРОНА А А УУ О А PEUX, FLUXNO 
СОММОХ ТАМВІ. ТАМВ2, ТАМВ5, ТАМВ4, ТАМВ5, ТАМВ6 

G [SON Hl. HI? Hš. H4. H5S, H6 

COMMON BIO1, 3102, BIO3, BIOJ. BIOS, BIOG 

EOS ION SIGI. SIG2. SIG5, SIG, SIGS, SIG6 

ПАО МАЈА ОУ 65. 63. 65 

ОТОМ УТ, хоу. У4, Ка, VO 

CONIMON WE W2 W3, WI W5, Wo 

EOALMON AL ADNKS A1 ASA 

АМОС COUNT, FREO 

NND IONEPSI, EPS?2, ЕР55, ЕР54, ЕР55, ЕР56 

COMMON QDOT, GEN, Y1 

COMMON TIMEO, TIMEI 

Meri VION SOLAR, EARTH, RE, DIST, ALT, FE, FA, MU 
COMMON SOLARI, SOLAR2, SOLAR3, SOLAR4, SOLAR5, SOLARó6 
SOM VON EARTH! EARTH? EARTHS, EARTHS, EARTHS, EARTHG 
COMMON ALBEDI, ALBED2, ALBED3, ALBED4, ALBED5, ALBED6 
COMMON ABS], АВ52, АВ53, АВ54, АВ55, АВ56 

(C OO Ms ON SACI, SAC2. SAC3 SAC4. SAC5, SACó 
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COMMON SUN, EGF S, Of, O27 Oconee 
COM NPO PSPP np P Sp ПО 


IFJ EQD AND (LEO ІШЕ 

DJ) = TUJ,K) + U4*(T(J,K-1) TE + 
L5*(TSTARI(I+1,3,K) - TSTARIG I,K) + 

& V1*(TAMBI - T(1,J,K)) + V4*(TAMB4 - T(IJ,K)) + 

V'6"(TAMB6 - T(1,J.K)) + WI + W4 + W6 + 

XI*(T(LLK)**4 - TAMBI**4) + 

XA*(T(1J,K)**4 - TAMB4**4) + 


бо 


ою бє 


БІСҢЕЛЕ(ТОТІ АХР ПАН БОЕ. 
D= TOLK) = CANARI- eee 


& LO*( TS РАК ТИ РА А БОНИ ЛАА E 

& V4"(TAMB4 - T(ILJ.K)) - VO (TAMBOG - TUJ,K)) + ХА + 
& WO * NJ*(I(LJ,R)**4 - ТАМВЧ““4) + 

& ХОЛИЈК) А ПВО КЕИ ОЕ Оа) 


ELSE IF((J.EQ.N) .AND. (LEQ.I)) THEN 
DJ) 7 T(LIK) + U4*¢(T,J.K-1 - TUE) + 


& US*(TSTARH(I- LJ.K) - TSTARI(LJ.K)) + 
& VOSCEA MB? - T(1J_.K)) + козе тава а SD ES) 
& V6*(TAMB6 - ТАЈК)) + W2 + W4 + W6 + 

& N2*(T(LJ.X)**4 - TAMB2**4) + 

& XJ*(T(1J,K)**4 - TAMBA**4) + 

& X6*(T(lJ.K)**4 - TAMB6**4) + YI*QDOT(I,J,K) 


ЕЅЕ ТЕ(О.ЕО КАРЕЛ A NREL Poy nies 
DU) = T(1J.K) ГАЗ К SUUS E 


& FO*(TSTARK(I-1,J,K) + TSTARK1+1,J,K) - 2*TSTARI(1,J.K)) 
& + VI*(TAMBI - T(IJ.K)) + V4*(TAMB4 - T(1.J,K)) + WI + 
& ма  XI*CT(LJK)**4 - TAMBI**4) + 
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& AJA*(T(IJ.K)**4 - TAMB4**4) - YI*QDOT(1J,K) 


ВЕ ЕО СТ.) АБ. TILTA) AND. (LGT. AND. 
ESI TM) THEN 
D(J) = T(1.J,K) + U4*(T(1,J,K-1) - T(1,J,K)) + 


& ВН + ISTARI I IR) - 2 ТЪТАКТ ИТ К)) 
& + V4"(TAMB4 - T(1,J,K)) + Wd + 
& AA"(T(LJ,RK)**4 - TAMB4**4) - YI*QDOT(I,J,K) 


ШЕЕЕ ТЕС ШЕОО) АА (ТЕСТІ) АХӘ С (ТІЛІМУ) THEN 
Е ТИК) = АТК) - Т) + 


& FO*(TSTARI(I-1.J.K) + TSTARI(I+1,J.K) - 2*TSTARI(I.J.K)) 
& ЕС ПІР ТРИ ЕПА у Ене + 
& W4 cO XO*(T(LJ.KO*74 - TAMB2**4) + 

& XA*(T(LJK)5*4 - TAMB4224) + Y12QDOT(LJ,K) 


ELSE IF((J.EQ.D .AND. (LEQ.M)) THEN 
DJ) T(LALK) * UJ*T(IJ,K-1) - TCLJ.K)) + 


& US*(TSTARI(I-1J.K) - TSTARI(LJ,K)) + 
& VI*C(TAMBI - Т КЕ) *- VJ(TAMBA - T(L,K)) + 
& V5*({TAMB5 - T(LJ.K)) + W1 + W4 + WS 

& ХІЗЧТА ІК). ТАМВІ%94) + 

& XA*(T(LJ.K)*54 - TAMB4*#4) + 


ОО (ГОА) ОИ iB S**4) ale ODOT(1.J,K) 


ELSE IF(((J.GT.1) AND. (J.LT.N)) AND. (1.EQ.M)) THEN 
D(J) 2 T(LJ,K) * UJ*(T(LJ,K-1) - T(L,K)) 9 


& US*(TSTARI(1-1,, K) - TSTARI(LJ,K)) + 

& V4*(TAMB4 - T(1.J.K)) + V5*(TAMBS - T(I,J,K)) + W4 + 
& W5 + XA*(T(L,K)**4- ТАМВА““А) + 

& X5*(T(1J,K)**4 - TAMBS5**4) + Y1*QDOT(I,J,K) 

ЕЕСЕ 


РО) = T(1,J,K) + C4*(T(LJ,K-1) - T(1,J,K)) + 
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& СО (ТЪТАКЦЕ ЦИР - ПЪТ КГ ЕТЕ; 
& W2"(TAMB2= TR) VS" A Dp D] P K ја 
& МО ТАМЕ И И y T 


& XA*(T(LJ.K)**4 - TAMB4*#4) + 
& XS*(T(IJ,K)**4 - TAMB$5**4) + YI*QDOT(IJ,K) 


ENDIF 
RETURN 


вр 


ate ale ate ate ate ale ato абл we ale ам ole uu ае ufo абе аби абе x< бо або фи «и фр е аи 
5 2,4 27 2,1 2,5 4% 2,4 7ұ% 29% 2% 245 7% 44% 24% 2)% 4% 4% 5,5 5% 29% 2,1 age sga зул aye age 


4% 47, ate ate we ate ate ate ats x". ate ate ate ate ale whe aks ate ate ale alo ale ale ale whe ale ole ate ate ate wale ate ate ato ја ale ate ale ate ate ale ate ate айде ГА we alo ale we Ме Ма ah ale we ale ~, 
+ 14055 2 2,4 241 2,1 2,2 2,1 2,5 5,2,4 2,5 25 2,4 24% 245 2,4 24% 2% 2,5% 54% 29% 24% 291 1,5 24% 79% 2,5 HGS 74% 7% 2,4 2,% 2,% 2,% 2)% 29% 2,5 2,% 2,4 


` v to ate we ate wie we ate ale whe we we ate we we ale we we sto te 
59555 55 2,% 2,5 645 2,5 9% 2): 2/1 2,% 64% 44% 24% 24% 54% 64% 29% 2)% 29% 2)% Hye Age тук бул буз зет Жуз суз 


SUBROUTINES ZULU]. ZULL2, AND ZULUs COMPUIB Ger IO OS 
ARRAY DTS THE Z DIREC HO: 


SUBROUTINE Z LECI 


SUBROUTINE ZULU] 

DIMENSION A( 30). B( 30), C( 30). D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAMMA( 30), QDOT( 30, 30, 30) 
DIMENSION T( 30, 30, 30), TSTARI( 30, 30, 30) 
DIMENSION TSTARX 30, 30, 30) 


REAL TINIT, ТУРАК УКО ГЕЕВ 

REAL CP KI RIO PYRI FO JI IPE pri i 

REAL EX LY DZ DEDXPSDERSSPEIMZ 

REAL- A,B, C. D, BETA CANINA 

REAL FLUXI, FLUX2, FLUX3 FLUX4A TIES FEG 
REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMB5, TAMB6 
REAL Hd 2 НН 
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REAL BIO], BIO2. BIOS, BIO4, BIOS, BIOG 

ПРАГ 5101. 3102, 5103, 5104, 51905, 5106 

ЈАТ (1, 23. 23, СА, 25 

BENI V1, V2, V3, V4. ¥5, V6 

REAL W1, W2, W3, W4, W5, W6 

eat Al, X2, X3, M4. A5, ко 

ШЕСІБЕР5І; ЕР5?> ЕР525, EPS3, ЕР55, ЕР56 

MeL ODOT, GEN, Yl 

FALC CTIME, TIMEO;, ИМЕ! 

fee SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
РЕСІ SOLARI, SOLAR?, SOLAR3, SOLAR4. SOLARS, SOLARG 
REAL EARTHI1., EARIHS EARTES, EARTIA. EARTHS. EARTIH6 
ime ALEBEDI. ALBED2, ALBED3, ALBEDJ, ALBED3, ALBED6 
ВЕР ABSI, ABS?, ABS3, ABSJ, АВ55, АВ56 

ReAL SACI. SAC2. SACS. SACI, SACS, SAC6 

moist N. ECEIPS, Ql. Q2, Q5, PERIOD, P] 

fee vt Pl, P2. PS, PJ, P5, P6 


ESI GER T. № > P. IF. L, IFPI LASI G; Q 
meteor COLNT, FREQ, ANSI, ANS2, ANS3S, VAL, NOUT 


ШО ГМ ОМ TEINIT, TSTAR1, TSTAR2, T, TEMP 

mobi ION CP KI, RHO. R. RI, R2, TIME; DELT 

Gs SON [N . TY LZ DELAXA. DELY, DELZ 

ШОО МОМ, bA MN. S.P. IF, L. IFPI, LAST 

НОММОХМ A,B,C, D, BETA, GAMMA, G, Q 

Morin Vee Ia PLE K?, FLUX3, FLUNG, FLUUXS, FLUX6 
COMMON TAMBI, TAMB2, TAMB3, TAMB4, TAMB5, TAMB6 
ЖОЛИ МОХМ Н' НП РБ H4, H5, H6 

COMMON BIOI, BIO2, BIO3, BIO4, BIOS, BIO6 

COMMON SIG1, SIG2, SIG3, SIG4; SIG5, SIGO 

ВИО Е 2, U3, U4, US 

WOM VION V1. V2, V3, V4, V5, V6 

COMMON WL W2, W3. W4, W5, W6 
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COMMON XI, X2. X3, X4. X5, X6 

COMMON COUNT. FREQ 

COMMON EPSI., EPS2, EPS3, ЕР54, ЕР55, ЕР56 

COMMON ODOL CE S] 

COMMON TIMEO, TIME] 

COMMON SOLAR, EARTH, RE, DIST, ALT, FE, FA, MU 

COMMON SOLARI, SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
COMMON EARTHI, ЕАВТН?, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI, ALBED2. ALBED3, ALBED4, ALBEDS, ALBED6 
COMMON ABS], АВ52, ABS3, ABS4, АВ55, АВ56 

COMMON SACI, SAC2, SAC3, SAC4, SACS, SAC6 

COMMON SUN, ECLIPS. Q1, Q2, Q3, PERIOD, PI 

COMMON PI. P2, P3, P4, P5, P6 


[F((K.EQ.1) AND. (1.EQ.1)) THEN 
D(K) = TSTAR2(LJ.K) + U3*(TSTAR2(1.J+1.K) - TSTARXIJ.K)) + 


& US*(TSTARI(I LK) - TSTARI(LJ,K)) Ж 
& VITAMBI- TSTAR MII K) ВОЗИ 
& + V6*(TAMB6 - TSTARX1J.K)) + W1 + W3 + W6 + 

& X1*(TSTAR2.J.K)**4 - TAMBI#*4) + 

& ХА(ТУТАВИЈК) 4 - ТАМВ3=У4) + 

& X6*(TSTAR2(1.J.K)**4 - TAMB6**4) + YI*QDOT(L,J.K) 


ELSE IF((K.GT.1) .AND. (K.LT.P)) АХР. (1.ЕО.1)) THEN 

D(K) = TSTARAIJ.K) - U3(TSTAR2(IJ- LK) - TSTARZ(LJ.K)) + 
U5*(TSTARI(I+1,J,K) - TSTARI(I.J,K)) + 
VI*(TAMBI - TSTAR2(LJK)) + V6*(TAMB6 - TSTAR2(1,J,K)) 
+ WI + W6 + XI*(TSTAR2(I,J,K)**4 - TAMB1**4) + 
X6*(TSTAR2(LJ,K)**4 - TAMB6**4) + YI*QDOT(LJ,K) 


ELSE IF((K.EQ.P) .AND. (LEQ.1)) THEN 

О(К) » TSTARXLJ,K) + U3*(TSTAR2(1,J + 1,K) - TSTAR2(1,J,K)) + 
& US*(TSTARI(I L.,K) - TSTARI(IA,K)) + 
& VI*TAMBI - TSTARX(UJ,K)) *: VA*(TAMBA - TSTARX(U.K)) 
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& СОИ Поне По а реа рок и ке Др тъ ААА + W6 + 


& рот TO TANDP Po O T YISQDOI(LI,K) 


ELSE IF((K.EQ.1) .AND. ((I.GT.1) .AND. (IL.LT. M)) THEN 

БО) = ТӨТАВДТІЕ) 2 US*(TSTAR2(IJ LK) - TSTARZGJ.K) 4 
FO*(TSTARI(I-14,K) - TSTARI(I  LJ,K) - 2?TSTARI(J,K) 
+ VI*TAMBI - TSTARXLA,K)) + 
V3*(TAMB3- TSTARZ(LJ.K)) - W1 * W3 + 


> P P» P fe 


HESSE GTA) AND. (K.ET.P)) AND. ((1.GT.1) AND. 
ERU T. ND THEN 
D(K) » TSTARXIJ.K) * U3TSTAR2Z(IJ  L.K) - TSTAR2Z(LJ,K)) 9 


& Jer psl AP] Tl] ролевите и О БАКТА) 
& EDTA NBI TSTARX TLE Tt WI + 
& ОКА ЕВУ) YIQDOT(IJ.FR) 


Bese TE((K.EO.P) AND. ((1.GT.1) AND. (1.LT.M))) THEN 
D(k) = TSTARXIJ,K) + U3*(TSTARIU.J+ 1.K) - TSTARXIJ,K)) + 


& КО u T ОО P” TSTAPHISIJK)- 2šTSTAR(IJK) 
& - VISTAMBI - TSTARXI.J,K)) + 

& VA (TAMBA - TSTAR2UJ.K)) + WI + W4 + 

& XI*(TSTAR2UJ.K)55*4- TAMBI**4) 4 

& X4*(TSTAR2(1.J.K)**4 - TAMB4**4) + YI*QDOT(IJ,K) 


ELSE IF((K.EQ.1) .AND. (I.EQ.M)) THEN 
D(K) = TSTAR2(IJ,K) + С3“Т5ТАВ2 + 1,К) - ТЅТАВ2(1.7,К)) + 


& U5*(TSTARI(I-1.J,K) - TSTARI(IJ,K)) + 
& VI*(TAMBI - TSTAR2(,J,K)) + V3*(TAMB3 - TSTAR2(I,J,K)) 
& + V5*(TAMBS - TSTAR2IJ,K)) + WI + W3 + W5 + 

& XI*(TSTARXI.LK)**4 - TAMBI?*4) + 
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& X3*(TSTARXI.J.K)**4 - TAMB3**4) + 
& XS*(TSTARXLJ,K)**4 - TAMBS**4) -- YI*QDOT(I,K) 


ELSE IF(K.G1 1) ASD КАТ РАЗВИТ ЕО СЛ THES 
D(K) = TSTARXIJ.R) * US (ISLAR2ODIT LK) - БАЊАЦ а 


& US*(TSTARICG-1,,K) - TSTARI(I.K)) = 
& VI*(TAMBI - TSTARXIJ,K)) *- V$*(TAMBS - TSTARX(I,J.K)) 
& + №1 + \5 + XI*(TSTARXIJ,K)**4 - TAMBI**4) — 
& X5*(TSTAR21,J,K)**4 - TAMB5**4) + Y1*QDOT(I,J,K) 
ELSE 
рК) = TSTARXIS K) + U3Ss(1STAR2TJ IK ШЕТ ЕРЕ. 
& U5%(TSTARI(I-1J.K)- TSTAR1I(IJ.K)) + 
& \1*(ТАМВ! - Т5ТАВ2(1.1.К.)) + VITAMBI. TSTAR2(1.J,K)) 
& + V5*(TAMBS - TSTAR2(1.J.K)) + W1 + W4 + W5 + 
& XI*(TSTARO(LJ.K)**4 - TAMBI**4) + 
X (ТОТАЛ К AE TANE а) 
& X5*(TSTARXLJ.K)**4 - ТАМВ5%%4) + Y1*QDOT(,1,K) 
ENDIF 
RETURN 
END 


SC BISOU LISEZ ЕО 


SUBROUTINE ZULU2 

DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAMMA( 30), QDOT( 30, 30, 30) 
DIMENSION T( 30, 30, 30), TSTARI( 30, 30, 30) 
DIMENSION TSTAR2 30, 30, 30) 


БЕЕАГТІХІП ТАКІ ЛО ИЕ ЕЕ 
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REAL CP. K1, RHO, R. R1. R2, TIME, DELT 

REAL LX, LY, LZ, DELX. DELY, DELZ 

REAL A, B, C, D. BETA. GAMMA 

Pee FLUX! PWUN2 FLUNS, FLU X4, FLUXS. FLUNG 

REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMB6 
REAL HI, H2, H3, H4, H5, H6 

REAL BIOI, BIO?2, BIO3, BIO4, BIOS, BIO6 

REAL SIGI. SIG2, SIG3, S1G4, SIG$, 5166 

ШОО Ul U2, 3, U4, U5 

REAL VI, V2. V3, V4, V5. V6 

REAL WI, W2, W3, W.1, W5, W6 

ЫСА XI, N2. Жз. ХА, Х5,Х6 

REAL EPS1, EPS2, EPS3, EPS4, EPS5, EPS6 

РЕК Г ODOT, GEN. Y1 

REAL CTIME. TIMEO. TIME! 

REAL SOLAR. EARTH, ALBCO, RE, DIST, ALT, FE, FA. MU 
REAL SOLARI, SOLAR?, SOLAR3. SOLAR4, SOLARS, SOLARG 
REAL EARTHI. EARTH?, EARTH3, EARTH4, EARTHS, EARTH6 
REAL ALBED!. ALBED2. ALBED3, ALBED4. ALBEDS, ALBED6 
REAL ABSI. АВ52. ABS3, ABS4. АВ55, АВ56 

REAL SACI, SAC2, SAC3, SAC4, SAC5, SAC6 

REAL SUN, ECLIPS, QI, Q2. Q3. PERIOD, PI 

REAL PI. P2. P3, P3, P5, P6 


ЕСЕК ИО Коме В. ТЕ ТЕРІ САО 
BSDEGERCOUNLI FREQ, ANSL, ANSE, A.NSS, VAL, МОСТ 


COMMON TINIT, TSTARI], TSTAR2, T, TEMP 

НОРА ПАЛО СРОКТТКНО, КУ ЕКІ ЕОТТІМЕ, БЕТ 

SOON Ewer, LZ, DELX, DELY,DELZ 

COMMO ЕРИМБЕМ, NXP.- IF, EPIFPIPEAST 

COMMON A, B, C, D, BETA, GAMMA, G, Q 
EONINIONISIISSMCELUNZ. FLUNSSELUXISFLUXNS, FLUX6 
COMMON TAMB1. TAMB?, TAMB3, TAMB4, TAMB5, TANB6 
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COMMON HI. H2. H3. H4, H5. H6 

COMMON BIOL, BIO2, BIO3, BIO4, BIOS, BIOG 

COMMON 5101. 5102, 5103, 8164, $165. $166 

COMMON UI TARU UI 

COMMON VI, V2, V3, V4, V5, V6 

COMMON W1, W2, W3. W4, W5, W6 

COMMON XLI, X2, X39. X5. 6 

COMMON COUNT. FREQ 

COMMON EPSI, ЕР$2, ЕР53, ЕР54, ЕР55, ЕР56 

COMMON QDOT, GEN, YI 

COMMON TIMEO. TIMEI 

COMMON SOLAR, EARTH. RE, DIST, ALT, FE, FA, MU 

COMMON SOLARI, SOLAR2. SOLAR3, SOLAR4, SOLARS. SOLAR6 
COMMON EARTHI. EARTH?2, EARTH3. EARTH4, EARTHS, EARTH6 
COMMON ALBEDI, ALBED2, AEBED3, ALBED4, ALBED>, ALBEDG 
COMMON ABSI, ABS2. ABS3. ABS4. АВ55, АВ56 

COMMON SAGINSAG2NSACs ОАА К С 

COMMON SUN, ECLIPS, Q1, Q2, Q3. PERIOD, PI 

COMMON P1. P2, P3. P4, P5. P6 


IF((K.EQ.1) .AND. (1.EQ.1}) THEN 
D(K) = TSTARQ(LJ,K) + Ul*(TSTARA(1J-I,K) + TSTAR2(1,J5+1.K) - 


& УУТУТАВИЈ КУ) + 5“ Т5ТАЕЦІР ШИК). ТОПАН! 
& + V3*(TAMB3 - TSTAR2XI.J,K)) + 

& V6*(TAMBG - TSTARX1J.K)) + W3 + W6 + 

& XO*(TSTAR2(IJ.K)**4 - TAMB3**4) + 

& X6*(TSTAR2(IJ,K)**4 - TAMB6**4) + YI*QDOT(I,J,K) 


ELSE IF(((K.GT.1). .AND. (K.LT.P)) .AND. (1.EQ.1)) THEN 
D(K) = TSTAR2AI,J,K) + Ul*(TSTARAI,J-1,K) + TSTAR2(,J + 1,K) - 


& 2*TSTARXI,J,.K)) - USTSTARI(I- LK) - TSTARI(LJ,K)) 
& + V6*(TAMB6 - TSTAR2(1,J,K)) + W6 + 
& X6*(TSTAR2(1,J,K)**4 - TAMB6**4) + Y1*QDOT(I,J.K) 
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ELSE IF((K.EQ.P) .AND. (1.EQ.1)) THEN 
DUK) = ТЕТАВХІТК) ғ Г.ІХТ5ТАВ2(1,7-1.К) - TSTAR2(ILJ 4 LK) - 


& 2*TSTARXJ.K)) + US*(TSTARI(I4 1,J,K) - TSTARI(LJ.K)) 
& + V4*(TAMB4 - TSTAR2I,J,K)) + 

& V6*(TAMBG - TSTAR21,J,K)) + W4 + W6 + 

& XA*(TSTARZ(IJ.K)?*4 - TAMB4**4) + 

& X6*(TSTARXI,,K)**4 - TAMB6**4) + Y1*QDOT(I,J,K) 


ASE IF(R.EO.1) AND. ((1.GT.1) AND. (.LT.M))) THEN 
Р(К) = TSTARXIJ,K) * UI«TSTARAIJ-1,K) - TSTAR2(LJ- Ll,K) - 


& 2*TSTAR2IJ.K)) + FOX TSTARI(I-1.J.K) + TSTARI(I+1.J,K) 
& - 2*?TSTARI(I,J,K)) + V3*(TAMB3 - TSTAR2I.J,K)) + 

& W3 + X3 (TSTAR2(LJ.K)**4 - TAMB3**4) + 

& Y1*QDOT(.J.K) 


ШЕРІР(СЕСТІІ) АХр (КІТР) AND. (1.GT.1) . AND. 
& (LLT.M)) THEN 

D(K) 2 TSTARXIJ.K) ^ UI*(TSTARXIJ-l.K) * TSTAR2(LJ * LK) - 
& 2STSTAR2UIK) = FOSTSTARIG-1J,K) 3 TSTARA(I- 1J.K) 
& - 2*TSTARI(IJ.K)) + Y1*QDOT(I,J.K) 


ELSE IF((K.EQ.P) .AND. (1.GT.1) .AND. (1.LT.M))) THEN 
D(K) = TSTAR2(1,J,K) + U1L*(TSTAR2(,J-1,.K) + TSTAR2(1J + 1.K) - 


& 2*TSTARXIJ.K)) + FO Т5ТАВ1(1-1,К) + ТЗТАВЦ + 1.1,К) 
& - 2*TSTARI(I.J,K)) + V4*(TAMB4 - TSTAR2I.J,K)) + 
& W4 + X44(TSTAR2I,J,K)**4 - TAMB4**4) + Y1#QDOT(I,J,K) 


EESE IF((K EQ.) .AND. (LEQ.M) THEN 

DE ЕКОО) + СЛРОРТА КА ЈА КОЈЕ 1БТАК2АЈ +1,К) - 
Z TSTARX LI K) + 56 АВЕ РА КАЈ )) 
ove (| AMB - TS АКО IBS) 
Vor GhavVibs= TSTAR2Z Орион МЕ, + 
AS*"(TSTAR2X(IJ,K)**3- TAMB3**3) + 


P P P P 
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& 


NS"*(TSTAR2IG, JK)? "4 - ТАМО НА ЈЕ ТЕОРИЯ 


ELSE IFR GT PASP KLIP ph r o М) ШИЕ 


ЦК) = TSTAR2(IJ.K) + СОІАИ) БИСЕ ЕЕ 


& 2*TSTAR2(,J,K)) + U5*(TSTARI(I-1,J,K) - TSTARI(L,K)) 
& + V5*(TAMBS - TSTARQ(UJ,K)) + W5 + 

& X5*(TSTAR2(1,J,K)**4 - TAMBS**4) + YI*QDOT(I,J,K) 
ELSE 


D(K) = TSTARXLJ.K) + LI*(TSTAR2(1,J-1,K) + TSTARX(IJ FEK): 


& ОТЛАРИ) СТАРТ TERTE 
& + V4*(TAMB4 - TSTAR2(1.J,K)) + 

& V5*(TAMBS - TSTAR2(1,J.K)) + W4 + W5 + 

& XA*(TSTARXIJ.K)**4 - TAMB4A**4) + 

& NS*(TSTAR2X(LJ.K)?*4 - ТАМВЗ““А) + YI*QDOT(IJ.K) 
ENDIF 

RETUR 

END 


SUBROUDENE-ZISIS 


SUBROUTINE 21143 

DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAMMA( 30), QDOT( 30, 30, 30) 
DIMENSION T( 30, 30, 30), TSTARI( 30, 30, 30) 
DIMENSION TSTAR2( 30, 30, 30) 


REAL TINIT ISTARI, ISTARy I TEMP 
REAL CP KIL RHO, RRR? TENDERI 
REALLCA ЊОЈ КЕ ОЕ БЕР DELZ 
ЕВ Ь ЧУ 
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NE PERE FTCA FLUAS. FLEXIT FLUAS, FLULXG 

КЕЛІ. ТАМВІ. ТАМВ2?, TAMB3, TAMB4, TAMB5, TAMB6 
GSL Hi H. H5. H4, H5. H6 

REAL BIO], BIO2, BIOS. BIO. BIOS, BIO6 

REAL SIGI, SIG2, SIG3, SIG4, SIG5, SIG6 

DEALUI US, U5, UJ, СЗ 

REAL VI, V2, V5, VÀ, V5, V6 

REAL W1, W2, W3, W4, W5, W6 

ШЕЛІ ХІ. А2, А3 А4, А5, хб 

DIE SWNEDSI. EPSX.EPSS, ЕР54, ЕР55, ЕР56 

REAL QDOT, GEN, YI 

ЕБАГСТТМЕ, ТІМЕО, ТІМЕІ 

mew SOLAR, EARTHS ALBGO, RE, DIST, ALT, FE, FA, ML 
REAL SOLARI, SOLAR2. SOLAR3, SOLAR4, SOLARS, SOLARG 
GL ESKINI EARTH? EARTHS, EARTH EARIHS EARTH6 
ITAL ALBEDI, ALBED?, ALBED3, ALBED4, ALBEDS, ALBEDG 
МЕЗІ АПр51, А852, АВ55, АВ5-, ABS5, ABSO 

BESE SACI SAC? SACS, SACA, SACS, SACO 

ПА УРА У S ECIIPS Q1.O2. OS PERIOD, PI 

EAL PI P? Po P3. PS. P6 


о m. Р. ТЕ ЕТЕРІ МЫТ С, О 
BRE L COL NI БВЕО, А А м ЛОГ 


GO J ON TISPPATSrTARI], ISTIAR2 E T. FEMP 

NOXIVS CP RIT RHO R RI. R2 TIME, DELI 

SONM MOSTA LOLZ, DELA. DELY, DELZ 

COMMOANTI K, M,N. P.IF, L, IFPI LAST 

COMMON A, B, C, D. BETA, GAMMA, G, Q 

(ОИ ЧИТАО АРХ, БОЕХ 3, РРЕХА, БИСАО FEUXG 
СОММОХ ТАМВІ, ТАМВ2, ТАМВ5, ТАМВ4, ТАМВ5, ТАМВ6 
ООММО а Ин Та. На Ноене 

COMMON BIOI, BIO2, BIO3, BIOJ, BIOS, BIOG 

CONMNIOSN SI[GI. SIG2, SIG5, SIG4, SIG5, SIG6 


COMMON TILU UO 

COMMON VI, V2, V3, V4, V5, V6 

COMMON MEN: KNIN o 

COMMON XI AL AS NIN N 

COMMON COUNT, FREQ 

COMMON EPS1. EPS2, EPS3, EPS4, EPS5, ЕР56 

COMMON QDOT, GEN, Y1 

COMMON TIMEO, TIMEI 

COMMON SOLAR. EARTH, RE, DIST, ALT, FE, FA, MU 

COMMON SOLARI, SOLAR?, SOLAR3, SOLAR4, SOLARS, SOLAR6 
COMMON EARTHI, EARTH?, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI. ALBED2, ALBED3, ALBED4, ALBEDS, ALBED6 
COMMON ABSI, ABS2, ABS3, ABS4, ABSS, ABS6 

COMMON SACILSAG2 S3634550: 095 s Ge 

COMMON SUN, ECLIPS, Q1. Q2, Q3, PERIOD. PI 

COMMON PI. P2, P3, P-1, P5, Ре 


IF(K.EQ.1) .AND. (LEQ.1)) THEN 
D(K) = TSTARAXLLK) t US*(TSTAR2X(LJ-L,K) - TSEAROCE RYE 


& US*TSTARI(I T LJ,K) - TSTARI(LJ.K)) + 

& V2SC(TAMB2: TSTARMNIJ K) £ VSS TA DS S SJ SF Е 
& + V6"(TAMBG - TSTARZ(LJ,K)) + W2 + W3 + W6 + 

& NO*(TSTAR2X(LJ.K)**4 - TAMB2**4) + 

& NOSNTSTAR2X(LJ,K)**4- TAMB3**4) + 

& N6*(TSTARXIJ,K)**4 - TAMB6**4) + Y1*QDOT(I,J.K) 


ELSE IF((K.GT.1) .AND. (K.LT.P)) .AND. (LEQ.1)) THEN 

D(K) = TSTAR2(1,J,K) + U3*(TSTAR2(I,J-1,K) - TSTARZ(L,K)) + 
U5*(TSTARI(1 + 1,J,K) - TSTARI(,J,K)) + 
V2*(TAMB2 - TSTARXIJ.K)) + V6*(TAMB6 - TSTAR2(LJ,K)) 
+ W2 + W6 + X2#(TSTAR2(1,J,K)**4 - TAMB2**4) + 
X6*(TSTAR2XJ,K)**4- TAMB6**4) + YI*QDOT(I,J,K) 


р р р Р 


PLSE TF((K.EQ.P SA SD (TEOM THE 
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ЕКВ А О У А (ПУТА ОН ЈЕ К)- ТУТАКА ЈК)) + 


& ЕГО ро ие сипеи 

& SO TAMB? ISTARA ID i V4" 1A MB4 - TSTAR2(1,K)) 
& + V6"(TAMB6 - TSTAR2(1.J,K)) + W2 + W4 + W6 + 

& X2*(TSTAR2(1,J,K)**4 - TAMB2**4) + 

& X4*(TSTAR2(1,J,K)**4 - TAMB4**4) + 

& X6*(TSTAR2(1,J,K)#*4 - TAMB6**4) + Y1*QDOT(I,J,K) 


ӘРЛЕТЕПЕСЕОЛУ-АХТ». (БӨЛ ІСТІ ЕЛЕУІ)) THEN 
D(K) = TSTARZ(I,J,K) + U3*(TSTAR2(I,J-1,K) - TSTAR2(1,J.K)) + 


& FO*?(TSTARK(I+1.J.K) + TSTARK(I-1,J,K) - 2*TSTARI(I,J,K)) 
& + V2“™TAMB2 - TSTARA(I,J,K)) + 

& V3*°(TAMB3 - TSTAR2L.J,.K)) + W2 + W3 + 

& m ОАР КОЧА МВО) + 

& Х54Т5ТАВ2(1./.К)%%4- TAMB3**4) -- YI*QDOT(LJ.K) 


ЕЕ GT 1) AND (KI LT.P Ахр (ШӨЛІ) Ах. 
& (LLT.M))) THEN 
D(K) = TSTAR2IJ.K) + U3*(TSTAR2(1,J-1.K) - TSTARXIJ,K)) + 


& Ors УЗ А АИ hee ome AL ) А ОАЗА КАЈ К)) 
& TA (TAMB? ТУТАЕВ) H W2 + 
& да ДЪ ПАК TIR 1- TAMB 3 7 Y1*ODOT(IJ,K) 


ELSE IF((K.EO.P) AND. ((1.GT.1) AND. (LLT M) THEN 

КН АСТЕ К) less GnsmAR21J-1.K)- TSTAR2(1.J.K)) + 
FO*(TSTARI(I1+1,J.K) + TSTARM(1-1,J,K) - 2?TSTARI(IJ,K)) 
+ V2*(TAMB2 - TSTAR2(I.J,K)) + 
V4*(TA MB4 - TSTAR2I,J,K)) + W2 + W4 + 
X2*(TSTARXIJ.K)**4 - TAMB2**4) + 
XA*(TSTARXIJ,K)**4 - TAMBA**4) + Y1*QDOT(1,J,K) 


e» 


fo P mm 


ELSE IF((K.EQ.1) .AND. (I.EQ.M)) THEN 
В ват лозата ва К) - ТЗТАВ(,,К)) + 
& U5*(TSTARI(I-1,J,K) - TSTARI(I,J,K)) + 


Ig» 


& V2"(TAMB2 - ITSTAR2(I J K); V3(TANDp T p Po n 
& T VSSCLAMBS - ЛЕ ЕЕ алва најрани Ве ЗАНИ 

& ОГТАЈ ОИ ЗАЕЛА ИВЕТ И 

& AS" (TSTARXI,IK)**4J- TAMB3**3) + 

& XS*(ISTAR2(.J.K)**4 - TAMB5**4) + Y1*QDOT(I.J,K) 


ЕТЗЕЛЕЦК GT. 1). AND (K.L p) A SD (T O И ИЫ: 

D(K) = TSTAR2(I J.K) + СЗ РОГАТ TSTAF ЕЛЕС 
LS*(TSTAR I-LI- TSTAR KLIK a 
V2*(TAMB? - TSTARX1,J,K) + YS5S*(TAMBS5 - TSTAR2(IJ,K)) 


> 2 фр № 


ВЕЗЕ 
ОК) = ІГ5ТАКОПТКУ-ТОІ Т ТАКИЕ ТКИ Б 
& Б ТУРА ИЕР ТОЛА О Ы 
& VO*(TAMB2 - TSTAR2X(IJ.R)) + V4*(TAMB4 - TSTAR2(1.J.K)) 
& Et VS TAN BS = ISIS IP) Ww ВИА eee K. 
& A2A(TSTAR2A ЕИО: 
& AA (TSTARXIJ,K)**J - TAMBJ4**J) + 
& AS*(TSTARZC J Re "S = ТАМ СВО) 


ENDIE 


RETR 
ESD 


SSK NESS А ЖЫЛ ДЕЛИ ЗА СИ ИЗИ ЖЕ К СИЕ ДЕ ЛИ ДЕЕ Н ДЕТ ie aks fe oe oe oe ok ЗИ ЛЕ ы СЕ ККЕ ЖЕК Де МА Пи ЗА МЕ МЕ УБ МЕ ZE PE МЕ МЕ УБ ОК МЕ УБ ЖЖЖ 


2 25 ТЫЛ АЗА e she om oie a ¿K SX z Na si ур Ус ма ум Ус мо Ме шо лилии 


SUBROUTINE TRIDAG IS A SUBROUTINE FOR SOLVING A SYSTEM 
OF LINEAR SIMLLTANEOLSEQLATIO>S HAY GCA [IRIDI GO AL 
COEFFICIENT MATRIX. THE ЕОСАТТОХЗ АКЕ ХАМВЕКЕРІР Ili o R 
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L, AND THEIR SUBDIAGONAL, DIAGONAL, AND SUPER DIAGONAL 
ESEPEICIENISAARE STORED'TN THE ARRAYS A, B. AND C. THE 
ПЕРА РИМА SOLUTION VECTOR TEMPO). TEMP(L) 1S STORED 
ПАЛЕ АККАЛ ТЕМР. 


SEDBIRSOUTINE TRIDAG 

DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30). GAM MA( 30), QDOT( 30, 30, 30) 
DIMENSION T 30. 30, 30), TSTARI( 30, 30, 30) 
BEEN SION TSTAR2( 30, 30, 30) 


REAL TINIT. TSTARI. TSTAR?, T, TEMP 

REAL CP, K1, RHO. R, R1, R2. TIME. DELT 

RAL LN, LY. LZ, DELX, DELY. DELZ 

REAL A. B. C, D, BETA, GAMMA 

REAL FLUXI, FLUN2. FLUX3, FLUXA, FLUXS, FLUXG 

REAL TAMBI, TAMB2. TAMB3, TAMB4. TAMBS5, TAMBG 
REAL HI, H2. H3, H4, H3, H6 

REAL BIOI. BIO2, BIO3, BIOJ. BIOS, BIO6 

REAL SIGI. SIG2, SIG3. SIG4. SIG3, SIG6 

EDS). U3.U3, U5 

REAL V1, V2. V3, V4, V5, Уб 

DES SEO N23 03. W5, W6 

ШЕЛІГІ NO UNS N4, X5, X6 

REAL EPSI. ЕР52, EPS3, EPS4. ЕР55, ЕР56 

REAL QDOT, GEN. YI 

REAL CTIME, TIMEO, TIMEI 

REAL SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLARI, SOLAR?, SOLAR3, SOLAR4, SOLARS, SOLAR6 
REAL EARTHI, EARTH2, EARTH3, EARTH4, EARTHS, EARTH6 
REAL ALBEDI. ALBED2, ALBED3, ALBED4, ALBEDS, ALBED6 
REAL ABS1, ABS2, ABS3, ABS4, ABSS, ABS6 

REAL SACI, SAC2. SAC3, SAC4, SACS, SACG 

REAL SUN, ECLIPS, Q1, Q2. Q3. PERIOD. PI 
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REAL PIL РА О АВА EO 


ГАТЕСЕК 1.2. RKR. MA p Ir papi oe 
INTEGER COUNT, FREQ. ANSI ANS2ANS35 VAP ХОШ 


COMMON TINIT, TSTARI, TSTAR2, T, TEMP 

COMMON CP, KI, RHO, R, R1, R2, TIME, DELT 

COMMON LX, LY, LZ. DELX, DELY, DELZ 

COMMON L,J, E, M, N, P, IF, L, IFP1, LAST 

COMMON A, B, C, D, BETA, GAMMA, G, Q 

COMMON FLUXI, FLUX2, FLUX3, FLUX4, FLUX5, FLUX6 
COMMON TAMBI. TAMB2, TAMB3, TAMB4. TAMBS, TAMB6 
COMMON HI, H2. H3, H4, H5. H6 

COMMON BIOI. BIO2, BIO3, BIO4, BIOS, BIO6 

СОММОХ 51610 51627 5165. Чен STONES 

COMMON UINUDNIECNINENINS 

COMMON ню 

COMMON WI, W2. W3, W4, W5, W6 

COMMON ХА? М Xh n S Sç 

COMMON COUNT. FREQ 

COMMON EPS1, EPS2. EPS3, EPS4, EPSS, EPS6 

СОМА ТОМ ОКО GEN tl 

СОММОХ ПМЕО. ПМЕ 

COMMON SOLAR. EARTH, RE, DIST, ALT, FE, FA. MU 
COMMON SOLARI. SOLAR2, SOLAR3, SOLARJ, SOLARS5, SOLARÓ 
COMMON EARTHI, EARTH?, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI, ALBED?2, ALBED3, ALBED4, ALBEDS, ALBED6 
COMMON ABS1, ABS2, ABS3, ABS4, АВ55, АВ56 

COMMON SACI) SAG@2, SACs SA G45 Ne5 a eo 

COMMON SUN, ECLIPS, Q1, Q2, Q3, PERIOD, PI 

COMMON PI, P2, P3, P4, P5, P6 


COMPUTE INTERMEDIATE ARRAYS BETA AND GAMMA 
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BETA(IF) = B(IF) 

GAMMA(IF) 2 D(IF) BETA(IF) 

Pep? = IF + 1 

DO 400 Q- IFPI, L 

ЕО) = В(Оу- А(О)“С(О-Г)/ВЕТА(О-1) 

GAMMA(Q) = ( D(Q) - A(Q)*GAMMA(Q-1) )'BETA(Q) 
400 CONTINUE 


* — COMPUTE FINAL SOLUTION VECTOR TEMP 

TEMP(L) = GAMMA(L) 

LAST = L-IF 

O40 G=1, LAST 

О А 

TEMP(Q) = GAMMA(Q) - C(Q)*TEMP(Q + 1) BETA(Q) 
MO CONTINUE 

RETURN 

EXD 


$. е 4 3, $, U 3, 4 s 4 е 4, U U е, 4, е, 4, J е $, 4 
2 25 54 5 >< > SE sk е ме та Зи КИ 25 ae OR Oe ча НО За ЛА Ме НК ЗА а ЗЕ Ма За а ЗА За За са 


во то а мо чо мо що «во аа аа Мое мо аи 9, 22 «е мә Жә до ай; ate ate во ote Йо чо чо do Йо чо ate ate do до ма Ло аа De ме ate Jo Mo do uo co o eo eo uo o e Mo lo Ж Wa we We Se de te oe whe Ae le te we whe We ate ate we 
748 74% 7% 74% 7,% 2;% бо поч 7,% 7,% 44% 7)% 9,% 74% 04% 974% 9;% 3)% 7ұ% 99% 24% 24% 94% 7,% 94% 7)% 7,% 7,% %;% брз Ура 665 6,9 Fee Hed вач ој уч вот еде жез буз ја „Ја 745 714% 24% 7)% 7ұ% 74% 7% 74% 94% 74% 94% 74% 9)% 749 73% ја 9,% 7,4 7% 7,4 жуз жуз tq tq дч и 


flict BROUTINE WILL BE CSED FOR PRINTING THE NODE 
ШЕЛІРЕКАТСКЕ5В 


SUBROUTINE OUTPUT 

DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAMMA( 30), QDOT( 30, 30, 30) 
DIMENSION T( 30, 30, 30), TSTAR1( 30, 30, 30) 
DIMENSION TSTAR2 30. 30, 30) 


КБ ЕЕ ГК АЛБ ТАЕ А АВА Т, TEMP 
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REAL CP, KI, RHO. R. RI. R2, TIME. DELT 

REAL LX, LY, LZ, DELX, DELY, DELZ 

REAL A, B. C, D, BETA. GAMMA 

REAL FEUX. FLUX? FLUAS FLENI TLO ХС 

REAL TAMBI, TAMB2, TAMB3, TAMB4, TAMBS, TAMB6 
REAL HI. H2, НА, H4, H5, H6 

REAL BIOI, BIO?2, BIO3, BIO4, BIOS, BIO6 

REAL SIGI, SIG2, SIG3, SIG4, SIG5, SIG6 

БЕТИ ЗА 

REAL V1, V2, V3, V4. V5, V6 

REAL W1, W2. W3, W4, W5, W6 

REAL STA oe 

REAL EPSI, EPS2, EPS3, EPS4. EPS5. EPS6 

REAL QDOT. GEN. Y1 

REAL CTIME, TIMEO. TIME! 

REAL SOLAR, EARTH, ALBCO, RE, DIST, ALT, FE, FA, MU 
REAL SOLARI, SOLAR?, SOLAR3, SOLARA, SOLARS, SOLAR6 
REAL EARTHI, EARTH?, EARTH3. EARTH4, EARTHS, EARTH6 
REAL ALBEDI, ALBED2, ALBED3. ALBED4, ALBEDS, ALBED6 
REAL ABSI, ABS?, ABS3, ABS4, ABSS. ABS6 

RESL SACI. SAC? SAGs SAG SNe 55516. 

REAL SUN, ECLIPS, Q1. Q2, Q3, PERIOD, PI 

REAL PI. P2. P3, P1, P5, P6 


INTEGER IvJ, K. M. N. PË IP b s HI nipa sin s 
INIEGER COLUNT, FREQ ANSI ANS2 AF53. 81 ЛУЦИ 


COMMON TINIT, TSTARI, TSTAR2, T, TEMP 

COMMON. CP..K1, RHO F io ЕБЕТ 

COMNMO N LX. LY, LZ, DEP > mp Pp np ss 

COMMON I,J, K, M, N Pal И РАИ 

COMMON A, B, C, D, BETA, GAMMA, G, Q 

COMMON FLUA, FLUX2, FLUX АЕ О FEC SER no 
COMMON TAMBI, TAMB2, TAMB3, TAMBI TITANES ai 
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COMMON HI, H2. H3. H4. H5. H6 

COMMON BIOI. BIO2, BIO3, BIO4, BIOS, BIOG 

@ONINMON 51601. 5162, 5103 5124. 5105, 5106 

ВО ВО LI, U2, U3, U4. U5 

COMMON VI, V2, V3, V4, V5, V6 

COMMON WI, W2, W3, W4, W5, W6 

SONXIMON X1, X2, X3, X4, X5, X6 

COMMON COUNT, FREQ. ANS1, ANS2, ANS3 

COMMON EPS], EPS2, ЕР53. ЕР54, ЕР55, ЕР56 

COMMON QDOT, GEN, Y1 

COMMON TIMEO, TIME! 

SONDMON' SOLAR, EARTH, RE, DIST. ALT, FE, FA, MU 

COMMON SOLARI, SOLAR2, SOLAR3, SOLAR4, SOLAR$, SOLAR6 
COMMON EARTIII, EARTH2, EARTH3, EARTH4, EARTHS, EARTH6 
COMMON ALBEDI, ALBED2, ALBED3, ALBED4, ALBED$5, ALBED6 
COMMON ABSI. ABS2. ABS3. АВ54, АВ55. АВ56 

О | 10Х 5АСТ. ЗАС2, 5403. 5АСА, 5АС5, 5406 

COMMON SUN. ECLIPS, QI. Q2. Q3. PERIOD, PI 

COMMON P1, P2. P3, P4, PS. PG 


MERSER CAN PRISNI OLIT IHE TEMPERATURES AT SELECTED 
Otic ORIEI NODES ARE DESIRED THEA THE USER MAY 
G GE TIIE WRITE AND PRINT STATEMENTS TO THE 
DESIRED NODES. 


ИЕ IME NE. DELT) GO TO 500 


PRINT* 
WRITE(1,*)’ TIME(SEC) (6,111) (6,6,) (16,6) 
Due il E. c 07 


= 


500 WRITE (1,501) TIME, T( 6, 1. 11), T( 6, 6, 6), T( 1, 6, 6) 
501 FORMAT(IX.FS.0,4X,F$.3,4X,F8.3.4X,F8.3) 
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THE ULSER CAN ALSO HAVE ALL XODETENEPERS TECE TERE TEE 

* OUT. IF A PRINTOUT OF ALE PEANES AND OR Viewer or 
NOT REQUIRED, THEN THE DO VOOFS 15 BE Cl Sr p rO 
ACCOMODATE TRE ЕВ ЕВЕ 


WRITE(6,*) ‘DO YOU WANT K-PLANE PRINTOUT? (FOR YES, ENTER” 
WRITE(6,*) ‘ANSI AS 1; FOR NO ENTER ANSI AS 0.) 

READ(6,*) ANSI 

IF(ANSI .NE. 1) GO TO 699 


WRITE(2.600) *TIME(SEC) = TIME 
GOO TOP S AJ AS] 0 


DOGOS P ale 


WRITE(2,*) TEMPERATURE DISTRIBUTION ON PLANE :' 
WRITE(2.610) 'K 2 .K 

610 FORMAT(IX,A3,I2) 
WRITE(2.620) ((TUJ.K), J=1.N), P= 1M) 

620 FORMAT(IX.F6.1,1 X.F6.1.1 X, F6.1,L X,F6.1.1 X, F6. 1.1 X,F6.1,1 X, F6.1, 
& IX.F6.1.1X.F6.1.1X.FG. 1, ЕС) 


СОС СОГ МЕ 


699 WRITE(6,*) ‘DO YOU WANT J-PLANE PRINTOUT? (FOR YES, ENTER’ 
WRITE(6,*) ‘ANS2 AS 1; FOR NO ENTER ANS? AS 0.) 
READ(6,*) ANS2 
IF(ANS2.NE. 1) GO TO 799 


WRITE(3,700) 'TIME(SEC) = TIME 
700 FORMAT(A12,F12.3) 


ооо ЕТ, 


WRITE(3,*) TEMPERATURE DISTRIBUTION ON PLANE :' 
WRITEQ.710) J = J 
NY RITE(3,*) жж oh of oe oe ok af «fe xxx Ж А Аа” 

710 FORMAT(1X,A3,12) 
WRITE(3,720) ((T(1.J,K), 15 1,M), K 2 P,1,-1) 

720 FORMAT(IX,F6.1,1X,F6.1,1X.F6.1,1X,F6.1,1X,F6.1,1X,F6.1,1X,F6.1, 
6 1Х.Е6.1,1Х,Е6.1,1Х,Е6.1,1Х,Е6.1) 

ж 


S CONTINUE 


799 WRITE(6,*) ‘DO YOU WANT I-PLANE PRINTOUT? (FOR YES. ENTER’ 
WRITE(6.*) 'ANS3 AS I; FOR NO ENTER ANS3 AS 0.) 
READ(6.*) ANS3 
IF(ANS3 .NE. 1) GO TO 899 


ША ШЕ 500) TIME(SEC) = ИМЕ 
EXXEEOR MA TCA 12.F 12.3) 


ВО = М 


Pell Es?) TEMPERATURE DISTRIBUTION ON PLANE: 
S RITE(41.510) 1 =I 
K а сан атанады ыы ала ааа 
EDEDEOR MEA І(1Х,А3,12) 
WRITE(4,820) ((I(LJ,K), J 2 1,N), K2 P,1,-1) 
EEOSEOR VAT(CIS,F6.L,IX,F6.1, I X, F6. Ll, I X,F6.], IX; F6.1, IX, F6.1, 1 X,F6.1, 
ESI SEO. INEO. LE IN,F6. LI, EX, F0. 1) 
+ 


suse CONTINUE 


% 


EE CI NLBIAUDSPABSGSTHE:.COUNTER TO ZEROG 
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599 COUNT = 0 
z 
КЕШ PSS 
END 
Ж 
eh MH se we ok Me sh ae ok ЖЖ Ж: Ж зы эй з к су. Ж зу. Me E Tk Ya Mh th me me fe k Ye Kk ЖИ ЖЖ ИЖ ии 


8, ЗН de whe ар sb. sf * <. . О +, ~ ~ > "s 9%, <$ О * x EJ 
ЗИМИ ИСА Изи Зи Зи Зи ЗМ ЗМ ok sts th oh Из и Ме x xx x xe Me ok oh oh эн ЖЖ ЗОЖ 


SUBROU TINEW AnD 

DIMENSION A( 30), B( 30), C( 30), D( 30), TEMP( 30) 
DIMENSION BETA( 30), GAMMA( 30), QDOT(30, 30, 30) 
DIMENSION T30. 30, У ТЕТЕ ОИ 
БЕМЕХ ОХ Ти 


REAL TISII А ВЕ Е АА 

REAL CP. RL RHO R RCR ЕВЕ И 

REAL EX TZ DI INA PI ID ЗЕЛЕ ЛОРА 

УВА А АВА НОРА вг 
REALTLUALTLUNSILUNS PED SCISDI eee «C 

REAL TAMBI. TAMB2. ТАМВЗ. ТАМВ4, ТАМВ5. ТАМВ6 
REACTIILE. H TOS lo 

ISEXICBIOL BIOSADIOSSDIOS BIOS DIOO 

REAL SIGHITSIO2. SIS ЕС 

Reale Дл р. 

БАБ ҮК УО 6 

REAL WI УО АА ONS а ВИ ОИ ВА 

REAL ADTA А856 

ҚЕА!І- ЕРӘ5І; ЕР52: ЕРӚЗЗ“ЕРЗ ЕБ ЕЕ 50 

КЕЗЕ ОРОЛ ЕЕ YI 

ЕГА СИЕ ЕО ИЕ ТЕ 

REAL SOLAR, EARTH, ALBCO, RE, DIS, Abie ge F o sP. 
REAL SOLARI, SOLAR2, SOLAR3, SOLARA, SOLARS, SOLAR6 
REAL EARTHI, EARTH2. EARTH3, EARTHA EARTHS EARTHC 
REAL ALBEDI, ALBED?2, ALBED3, ALDED4, ALBEDS, ALBEDO 
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REAL ABS!, ABS2. ABS3, ABS4, ABSS, ABS6 

REAL SACI, SAC. SACS. SACA SACS SACO 
REAL SUN, ECLIPS, Q1, Q2, Q3, PERIOD, PI 
REAL P1, P2, P3. P4, P5. P6 


ENSE:GERI,J-K, M, N, P, IE, L, IFPLEAST, С, 
BSIEGER COUNT, EREQ, ANSI, ANS2, ANS, VAL, ХОСТ 


ВМО TINIT, ISTARI, ISTAR2, T, TEMP 

ge JViON CP, KT, RHO, R, RI, R2, TIME, DELT 

ES DIONLSLY.LZ, DELN, DELY, DELZ 

MERRION TLI KE ALN. PTF, L, IFPI; LAST 

OAMI OAN A B.C. D, BETA, GAMMA. G, Q 

ES no DLE NLLUNXSFLUNXS.PEUNA,ELUNXS,FLUXO 
ВАНО ГА ЗВ, ТАМВ2, ГАМВЬ. ТАМВА. ТА МВО, ТА МВО 
МЕСО“ И: НЗ. НЯ, Нэ, H6 

eo SION BIO!, BIO2, BIOS, BIO4, BIOS, BIOG 

ща ЦЕМОХ УТО]. 3102. $1G3, 516-4, $165. 53166 

ОХ Е. 62, (5. (4. (5 

ВОГ ОХ \1, \2.%5. V4. V5. VO 

ОЛОМ Ур W2, W5, W4. W5, W6 

BG IIO NONI, NO, X5, XJ. X5. X6 

BO DIONCOUNT.EREQ. ANSI, ANS2, А №53, VAL, ХОСТ 

BS IUOS EPSL EPS?. ЕР55. ЕР54, ЕР55, ЕР56 

go YON ODOT GEN. Y] 

COMMON TIMEO, TIME! 

FONIO SOLAR. EARTH, RE DIST -ALT FEFA MU 
BIOS SOLARIS SOLAR2, SOLAR3, SOLAR4, SOLARS, SOLAR6 
EON MOS TEARIHI,EARTH2, EARTH3, EARTHA, EARTH5, EARTHO 
BOSISIONSADPDBEDI P ALBEDO, ALBED3, AEBED4, ALBED5, ALBED6 
СОМ МОХ АВ51, АВ52, АВ53, АВ54, АВ55, АВ56 

COMMON SACI. SAC?, SAC3, SAC4, SACS, SAC6 

MeO NO wesw ECEIPS: OI, Q?, Q5, PERIOD, PI 


COMMON PIP? PS FARE В 


VARIOUS NODES WERE SELECTED TO SAN IFES Terie heat 
DIFFERENCES. THE USER CAN CHANGE THE SELE TED ORE 
IF DESIRED. IF THE USER GHANGES TE SELECTED ODES 
THE NODES SELECTED IN PROGRAM EXPLIC] I aD Dist ers 
CHOSEN IN PROGRAM VALID SHOULD BE ADJUSTED ACCORDINGLY. 


ТОВ! = О 0) - ТМТ 
TDIFB2 = T( 6, 6,11) - TINIT 
TDIFB3 = Tori 6) = TI AIT 


ТРОН МЕ NE ВЕЕТ ОО ТО 
ЕК 


WRITE(S.*) ‘TIME(SEC) TDIFB! TDIFB2 TDIFB3’ 
А зс eo ии 


920) WRITE (3.930), TIME, TDI Ee ise РЕКЕ 


г 


950 РОК МА А РОЗА ВА ВА ВИ ТАЕ 


REINITIALIZING THE COUNTER Were 


CO = 0 


RETURN 
END 


Ж ж:х ЖЖ Ж 2 о ЖЖ ужи 


x s x 


* w T 5 + 2, КА ~ ~ +J. КА г, 
ХАСО ИА ie he ie si и Ма МА Ла ЛИ МА ДОЗИ НОЕ Де 3 ТК x dk d До ПК ЖК ЛА ПК КОЕ МА СК ЗА oh oe oe OG ok oe ЗА ПК oR oe ЛЕ МК НК ПЕ ЛЕ Ne € 


186 


APPENDIX E 
PROGRAM VALID 


PROGRAM VALID 


ф жос Ж ЛЕ ЛЕ Ж ЛЕ Ес 2 ЛЕ сы ЛЕ ЖЕ @ Б ЖО ЕК 5 ЖО ЖАТ КАКВА АА КА КАА 


можа кажа кик ке ааа оф йй Ж Л иф k X X $ Жї їй йй $ $ 5 $ в кое каки: а 


* 


ж 


ж 


Ф 


ж 


VALIDATION PROGRAM 


ШИЕРІЕЕРОЗЕ ОҒ ПІ1І5 PROGRAM IS TO VALIDATE THE 
TEMPERATURES OBTAINED BY THE EXPLICIT AND BRIAN 
PROGRAMS. CONDITIONS USED IN THE VALIDATION PROCESS 
КИ CONSTANT SURFACE HEAT FLUX ON THE LEFT FACE 
ESSDUSDUINBATIC CONDITIONS ON-NEL OTHER SURFACES. 


INI S PROGRAM PROVIDES THE CLOSED-FORM SOLUTION OF 
ieee | HEAT CONDUCTION IN A SENII-INFINITE SOLID. 
BUSCLOSED-FORM SOLUTION FOR CONSTANI SURFACE HEAT 
eels FROM INTRODUCTION TO HEAT TRANSFER’ BY 
ROPEK AND DEWITT. THE RATIONAL APPROXIMATION FOR 
 БЕПККОКТГХСТТОХ IN THE CLOSED-FORM SOLLTION IS FROM 
ЕЕ DBOOR OF MATHEMATICAL FUNCTIONS” EDITED BY 
ABRAMOWITZ АХР Б5ТЕОСХ. 


ТЕТЕ. 


фес 


ж 


«кф Жс 7k k k k 22 2 кй» t: 5: & 5: $ k $ 2 $ G X $ $ $ k 2 E X E fk k 5 Yk 5 $ G 52 $ AE fx SE Ж: Ж: Җ-: k у у уй у 5 I 4: жож 5. k k 5: X Ë 


REAL Al, A2, A3, A4, A5, B, C, D, P, V1, V, W1, W, PI 
REAL TDIFV, FLUX, ALFA, K, RHO, CP, TIME, DELT, Y, DELY, ERRFC 
INTEGER COUNT, FREQ 
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THE USER WILL READIN Bak NIL pP | 


WRITE Е СТЕБ оа 


WRITE(*,*) ‘ENTER K, THE THERMAL CONDUCTIVITY.’ 
READ(*,*) K 
WRITE(*,*) 'ENTER RHO, THE MATERIAL DENSITY.' 


МЕРЕ?” ENERGIE SPEC EG RIEA 

eye Bile 

WRITE(*,*) ENTER ¥, THE DISTANCEA SE ob ih ee iG 
READ 

WRITEC,*) ENTER DELY, THE INCREMENT IN THE DIRE TLES 
КЕТЕРІ 


CONSTANTS VSED THROUGHOUT ie NO Gino. 


NOTE: VALUES ЕОК А АЛ АЗ, АЧ РО А У ОВО ЕНТОНИ 
“HANDBOOK OF NLATHEMATICAPI U NC TIO SS 


A] = 0595599 
A2 = -0.2S4496736 
АЗ = 1.421413741 
А4 - -1.455152027 
А5 - 1.061405429 
P = 0.3275911 

РІ = 3.141592654 
ПЕ =) 
COUNT = 0 


CALCULATION OF CONSTANTS LSCED THROUGHOUT TREPROCR Y 
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Severe OLATIONS FOR ALEA, Bb Gabe aN ОЕГТ АКЕ ҒЕКОМ 
I IRODLUCITIO TOTI Ji Siete EQOUATIONS FOR 
VILAND WI ARE FROM "HANDBOONR OF MATHEMATICAL 
Be SCTIONS’: 


ALFA = K'(RHO#CP) 

B = 2*FLUX/K*SORT(ALFA‘PI) 
C = FLUX*Y:K 

D = Y**2'(4*ALFA) 

DELT = 0.125%(DELY**2)'ALFA 
VI = P^Y (2*SQRT(ALFA)) 

WI = Y (2*SQRT(ALFA)) 


ANAT STATEMENTS FOR YARIOLS PARAMETERS USED 
IN PROGRAM 


PRINT? 
PRINT? 
NV R ] BEI m . E "ж Me chk sh Mee ch se sk oe аси сис СЕ у сизи су саси снн зр сис иод си сизи си су сиз с: эк сйс НЕ ЕНЕ 
WRITE(*.*) (FLUX = РІС 
WRITE(*.*) К = К 
WRITE(*,") ‘RHO=’.RHO 
БІБРІТЕ( 9) :СР- „СР 
WRITE(?,*) Y 2 ,Y 
WRITE(*,*) DELY 7 ;DELY 
WRITE(*,?) 'DELT 7 ,DELT 
ARITE *) ‚зи Та ти Se oe Sk ee Sh oh ie eh a ke He ee HH He eK Ж Ж: Ж ж Жек Ж: Ж: Ж: He He HK eH у Ж: з” 
PRINT* 
PRINT* 


— 


CALCULATION OF TDIFV 


КОШЕ FOCS TOSS FORY. W, AND ERRFC (COMPLEMENTARY 
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ж 


ж 


25 


ERROR FUNCTION ARE PROM НА НОВО ОВАКО 
MATHEMATICAE FUNCTIONS . THE EOUAIIONF epi. 
IS FROM “INTRODUCTION TO HEAT ТЕЗ РЕВ 


10 COUNT = COUNT + 1 
TIME = TIME + DELT 
V = 11 + ViSQRT(TIME)) 
W 2 WLSQRT(TIME) 


& (ASV S) EAPC) 
TDIFV = B*SQRIT(TIME)*EAPECD TIME) C ERRI 


PRINTING TINE МОГ ОГГЕ 


IF( TIME .NES DELI GO ee) 
УВПЕС S) TI (r mipi E 
WRITE(*.") --------- ----- | 

РАМА АИ АТА Ричи ОНА 


20 РОК и 


30 IT(COU NI NE-LRSEOQUOOSTOJ00 
PRINT 46, VINES DIES 


JO РОК I РО Ано) 


МИА ЕВА О THE COUNTERTOZER® 


50 СОСХТ - 0.0 
* NOTE: THE VALIDATION TIME PERIOD WAS ARBITRARILY TAKEN 
+ AS 3600 SECONDS (1 HOUR). 


we 
27“ 


100 IF(TIME Л.Е. 3600) GO TO 10 


~. 
° 


STOP 
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10. 
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